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ABSTRACT 


Methods  of  exterior  ballistics  applicable  to  utilization  in  modern 
airborne  fire  control  system  design  are  documented  herein.  The  fund¬ 
amentals  of  exterior  ballistics  are  included, along  with  a  description  of 
currently  used  ballistic  and  aerodynamic  notations,  a  discussion  of  the 
limitations  of  the  semi-empirical  aerodynamic  force  and  moment  system, 
and  methods  of  preparing  aerodynamic  data  for  use  in  trajectory  compu¬ 
tation.  Tutorial  material  is  provided  to  give  the  reader  an  understanding 
of  windage  jump  caused  by  the  complicated  angular  motion  of  a  spinning 
projectile.  The  Siacci  method  is  described  and  means  for  improving  its 
accuracy  are  developed.  Six-degree -of -freedom  equations  are  derived 
in  several  different  formulations  for  exploratory  studies  and  digital  com¬ 
puter  computations,  and  the  development  of  approximate  equations  for 
rapid  evaluation  of  trajectory  tables  is  included.  Methods  for  calculating 
trajectory  initial  conditions  are  provided  for  shells  fired  from  a  turreted, 
gatling  gun  in  a  maneuvering  aircraft,  and  the  problems  of  ballistic  and 
kinematic  prediction  are  discussed  briefly.  The  material  covered  he  -ein 
should  provide  personnel  in  the  Air  Force  and  in  industry  with  sufficient 
knowledge  of  exterior  ballistics  for  advanced  fire  control  system  design. 
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SECTION  I 


INTRODUCTION 

For  the  design  of  a  modern  airborne  gun  fire  control  system,  a 
knowledge  of  exterior  ballistics  is  a  necessity.  The  complicated  motion 
of  a  spinning  projectile,  fired  crosswind  from  a  gun  in  a  high-speed 
aircraft,  and  the  accurate  prediction  of  such  motion  requires  careful 
study. 

These  statements  are  borne  out  by  consideration  of  the  magnitudes 
of  aerodynamic  forces  and  moments  acting  upon  a  projectile  in  flight.  It 
is  almost  incredible  that  the  aerodynamic  drag  on  the  20-mm,  M56 
round  fired  at  sea  level  at  4000  ft/sec  is  about  20  lb,  and  that  the  left 
force  at  right  angles  to  the  velocity  vector  is  about  9  lb  if  the  angle  of 
attack  of  the  projectile  is  3  deg.  The  20-mm  round  itself  weighs  only 
0.22  lb.  The  effect  of  the  lift,  and  other  smaller  forces,  is  to  deflect 
the  motion  of  the  projectile  away  from  its  initial  direction  of  motion,  and, 
as  a  rule  of  thumb,  this  deflection,  called  the  windage  jump,  is  one  milli- 
radian  (mr)  per  initial  degree  of  angle  of  attack.  In  the  example  above, 
the  windage  jump  is  3  mr  and  one  might  ask  how  it  is  that  a  9-lb  force 
acting  at  right  angles  to  the  velocity  vector  of  a  0.22-lb  projectile  pro¬ 
duces  only  a  3-mil  deflection.  The  answer  is  that  the  lift  force  preces- 
ses  with  the  projectile  about  the  velocity  vector  and  nearly  cancels  out. 

That  a  careful  study  of  the  six-degree -of-freedom  motion  of  a  pro¬ 
jectile  is  necessary  should  now  be  apparent.  The  deviations  from  straight- 
line  motion  of  a  spinning  projectile  in  flight  are  due  to  complicated  inter¬ 
actions  of  aerodynamic  forces  and  moments,  gravity,  and  gyroscopic 
effects.  For  accurate  prediction  of  projectile  motion,  these  effects 
must  be  taken  into  account. 

This  report  is  limited  to  an  investigation  of  the  motion  of  a  spin¬ 
ning  projectile  fired  from  a  rifle.  Finned  missiles  and  rockets  are  not 
included. 
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Whereas  accurate  six -degree -of-freedom  projectile  trajectories 
can  be  calculated  directly  from  Newton's  laws  of  motion,  computations 
of  this  type  are  currently  not  believed  feasible  in  airborne  fire  control 
applications  because  excessive  computer  time  is  required.  Careful 
study,  however,  leads  to  the  selection  of  appropriate  algorithms  which 
yield  trajectory  data  within  the  error  budget  alloted  to  exterior  ballistics 
calculations,  and  which  can  be  utilized  in  airborne  applications. 

Because  of  the  9  or  1 0-year  lapse  of  interest  by  the  Air  Force  in 
exterior  ballistics,  it  is  felt  that  a  presentation  of  the  fundamentals  of 
the  subject  is  not  out  of  place  in  this  report.  New  personnel  in  the  Air 
Force  and  in  industry  are  not  likely  to  be  knowledgeable  in  this  subject, 
nor  is  the  necessary  literature  readily  available.  The  best  available 
book  on  exterior  ballistics  is  the  one  by  McShane,  Kelley,  and  Reno 
(Ref.  1),  but  this  book  is  not  well  suited  to  current  needs.  Other  useful 
information  is  to  be  found  in  reports  written  mostly  in  the  1950's  and 
earlier  which  may  be  hard  to  obtain,  such  as  in  Refs.  2  through  20. 

Most  of  this  work  is  based  upon  the  early  work  of  Fowler,  Gallop,  Lock, 
and  Richmond  (Ref.  21), and  upon  the  works  of  Nielsen  and  Synge  (Ref. 
22),  and  of  Maple  and  Synge  (Ref.  23).  A  useful  book  on  fire  control  is 
NAVORD  Report  1493  (Ref.  24),  and  the  book  by  Davis,  Follin,  and 
Blitzer  (Ref.  25)  provides  information  applicable  to  the  exterior  bal¬ 
listics  of  shells.  The  notation  of  Ref.  25  is  somewhat  different  than 
used  elsewhere,  but  the  book  provides  very  good  insight  into  projectile 
motion.  After  burning,  the  motion  of  a  spin-stabilized  rocket  is  the 
same  as  that  of  a  shell.  Two  other  books  on  exterior  balliscics  are  the 
ones  by  Moulton  (Ref.  26)  and  Bliss  (Ref.  27). 

The  aerodynamic  force  and  moment  system  in  current  use  is  semi 
empirical.  Whereas  the  functional  form  of  the  forces  and  moments  may 
be  partially  derived  from  considerations  of  dimensional  analysis  and 
projectile  symmetry,  aerodynamic  theory  has  failed  to  produce  accurate 
expressions  for  these  forces  and  moments.  Aerodynamic  data  used  in 
trajectory  computation  is  acquired  from  wind  tunnel  measurements  and 
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free -flight  tests.  The  arguments  of  dimensional  analysis  and  symmetry 
which  are  used  co  develop  the  aerodynamic  force  and  moment  system  in 
use  are  presented  in  Section  II.  Discussion  of  the  measurement  of  aero¬ 
dynamic  data  is  excluded.  Also  included  in  Section  II  is  a  comparison  of 
the  ballistic  K  notation  and  the  aerodynamic  C  notation,  physical  identifi¬ 
cation  of  the  various  aerodynamic  forces  and  moments,  a  discussion  of 
weaknesses  in  the  force  and  moment  system,  and  modeling  of  aerodynamic 
data  for  computer  utilization. 

Section  III  is  tutorial  and  should  provide  the  reader  with  the  neces¬ 
sary  background  to  understand  the  nature  of  the  angular  motion  of  a  spin¬ 
ning  projectile.  An  understanding  of  the  fundamental  causes  of  windage 
jump  follows  from  knowledge  of  projectile  angular  motion. 

Equations  of  motion  suitable  for  computation  of  trajectories  on  a 
computer  are  given  in  Sertion  IV.  Included  are  a  matrix  formulation 
and  an  Euler  angle  development  of  the  six-degree-of-freedom  equations. 
Also  included  is  a  set  ol  approximate  equations  which  are  useful  for 
rapid  generation  of  traj  sctory  data  when  computer  time  is  a  factor. 

The  basic  Siacci  method  is  presented  in  Section  V,  along  with 
methods  to  improve  its  accuracy.  The  Siacci  method  is  a  good  candi¬ 
date  for  onboard  utilization  because  it  provides  an  approximate  closed 
form  solution  to  the  trajectory  computation  problem. 

Equations  for  initial  conditions  for  trajectory  calculations  are  de¬ 
rived  in  Section  VI.  These  equations  apply  to  a  turreted,  gatling  gun 
fired  from  a  moving  aircraft.  They  are  possibly  more  complicated  than 
necessary,  but  they  may  be  simplified  as  any  application  permits.  A 
brief  discussion  of  the  problems  imposed  by  kinematic  prediction  is  also 
included  in  this  section. 
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SECTION  II 


AERODYNAMIC  FORCES  AND  MOMENTS 


1 .  General 

The  aerodynamic  forces  and  moments  acting  on  a  projectile  in 
flight  result  from  frictional  forces  and  pressure  distributions  over  the 
projectile  body  which  in  turn  result  from  the  motion  of  the  projectile 
through  the  air  mass.  A  complete  and  accurate  solution  for  projectile 
motion  would  therefore  involve  solution  of  the  equations  of  fluid  flow 
around  the  projectile.  In  practical  applications,  this  is  beyond  the  cap¬ 
abilities  of  any  available  digital  computer,  but  fortunately  this  is  not 
necessary.  The  alternative  is  a  semi-empirical  approach  whereby  aero¬ 
dynamic  forces  and  moments  are  measured  in  wind  tunnels  and  by  means 
of  free-flight  tests.  This  aerodynamic  data  is  modeled  in  a  form  suit¬ 
able  for  computer  use  and  employed  in  trajectory  calculations.  The 
means  by  which  this  data  is  modeled  are  treated  in  this  section. 

Dimensional  analysis  is  a  study  -t  certain  mathematical  relation¬ 
ships,  explainable  in  terms  of  the  d  .lensionality  of  measurements, 
which  exist  between  physically  measurable  quantities.  Physically  mea¬ 
surable  quantities  ha  e  dimensions  such  as  mass  (M),  length  (L),  Time(T), 
or  velocity  (LT  *),  and  physical  laws  are  independent  of  these  units  of 
measure.  Certain  restrictions  are  imposed  by  nature  upon  the  functional 
form  of  mathematical  relations  describing  such  laws,  and  dimensional 
analysis  identifies  these  restrictions.  In  the  semi-empirical  develop¬ 
ment  of  the  aerodynamic  force  and  moment  system  applicable  to  a  spin¬ 
ning  shell,  arguments  based  upon  dimensional  analysis  and  symmetry 
play  an  important  role. 

A  shell  possesses  two  types  of  symmetry,  namely  rotational  and 
mirror.  An  object  has  rotational  symmetry  if  "it  looks  the  same" 
when  it  has  been  rotated  through  a  given  angle  about  a  particular  axis. 
Similarly,  an  object  has  mirror  symmetry  if  it  looks  the  same  under 
certain  conditions  when  viewed  in  a  mirror.  Dimensional  analysis  and 
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symmetry  are  used  to  derive  the  functional  form  of  the  aerodynamic 
force  and  moment  system  in  terms  of  the  aerodynamic  or  ballistic  coef¬ 
ficients. 

Since  the  approach  used  to  mathematically  describe  aerodynamic 
forces  and  moments  is  semi-empirical,  it  has  certain  weaknesses. 

These  weaknesses  are  discussed  briefly. 

The  approach  herein  is  to  first  mathematically  model  the  force  and 
moment  system  and  then  to  identify  the  various  forces  and  moments.  One 
might  expect  the  analysis  to  proceed  in  the  other  direction,  but  use  of 
this  systematic  approach  uncovers  unsuspected  forces  and  moments. 

The  systematic  approach  should  not  trouble  the  reader  since  he  will  al¬ 
ready  be  aware  of  the  existence  of  certain  forces  such  as  the  drag,  the 
lift,  the  Magnus  force,  the  Magnus  moment,  etc.,  and  their  appearance 
in  the  equations  to  be  developed  will  be  no  surprise. 

This  section  also  contains  a  comparison  of  the  different  notations 
which  are  in  current  use  for  the  description  of  the  aerodynamic  forces 
and  moments  which  act  upon  a  projectile  in  tree  flight.  The  terminology 
used  in  the  literature  to  describe  the  aerodynamics  of  the  operational 
20-mm,  M56  round  is  the  K,  or  ballistic,  notation,  whereas  that  used  to 
describe  new  20-mm,  25-mm,  and  30-mm  rounds  under  development  is 
the  C,  or  aerodynamic,  notation.  Much  of  the  old  ballistics  literature 
which  employs  the  K  notation  is  still  of  interest,  but  the  new  C  notation 
is  used  by  wind  tunnel  and  free-flight  test  range  personnel  in  the  col¬ 
lection  of  new  ballistics  data.  Familiarity  with  both  notations  is  manda¬ 
tory. 

Choice  of  an  appropriate  notation  for  ballistics  is  subjective.  The 
C  notation  is  somewhat  awkward  in  use  in  mathematical  developments, 
although  it  does  exhibit  a  lucidness  of  meaning  which  is  absent  in  the  K 
notation.  The  two  systems  of  notation  are  mathematically  equivalent, 
and  the  K  notation  is  used  throughout  this  report,  since  most  of  the  lit¬ 
erature  upon  which  this  report  is  based  is  written  in  terms  of  the  K 
notation. 
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2.  Dimensional  Analysis 

The  subject  of  dimensional  analysis  is  an  important  part  of  the 
basic  theory  of  aerodynamics,  and  since  it  is  well  documented  (see,  for 
example,  Refs.  1,  28,  and  29),  it  will  not  be  the  purpose  of  this  section 
to  present  a  detailed  account  of  the  subject.  Rather,  results  derived 
from  the  theory  will  be  applied  to  the  problem  at  hand. 

If  a,  (3,  ...,  co  are  physically  measurable  quantities  (with  dimen¬ 
sions)  which  satisfy  an  equation  which  represents  a  physical  law,  such 
as 


cj>(a,  (3,  .  . .  ,  co)  =  0 

this  relation  should  also  be  satisfied  if  the  units  of  measure  of  a,  (3,  .  .  .  ,  to 
are  changed  (e.g.  ,  from  grams,  centimeters,  and  seconds  to  pounds, 
feet,  and  hours).  This  is  reasonable,  since  one  would  not  expect  the 
functional  form  of  <j)  to  change  under  such  a  transformation.  A  theorem 
of  dimensional  analysis,  known  as  the  Buckingham  li  -  theorem,  follows 
from  this  assertion.  Crudely  put,  the  theorem  states  that  the  relation 
<j>(a,  (3, . .  .  ,  co)  =  0  can  be  replaced  by  another  dimensionless  relation 


where  n  is  the  number  of  parameters  in  the  set  of  measurable  quantities 
a,  (3,  ...»  to,  and  m  is  the  number  of  fundamental  dimensions  (pounds, 
feet,  seconds,  etc.).  The  ir,  for  i  =  1,  2,...,  n-m,  are  dimensionless 
products  of  the  measurable  quantities  of  the  form 

a.  b.  z. 

1  o  1  1 

n\  =  a  (3  ...  co 

One  would  expect  intuitively  the  existence  of  a  function  such  as  4*. 
The  nature  of  a  physical  law  is  independent  of  the  units  of  measure  by 
which  it  is  described,  and  a  function  such  as  <{j  should  be  unaffected  by 
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changes  in  units.  Since  the  products  u\  (for  i  -  1,  2,  n-m,are 

dimensionless,  they  do  not  change  value  with  changes  in  units  and  it 
follows  that  .Jj  is  unaffected  by  such  changes.  The  function  <f>  is 
essentially  a  rearrangement  of  t|*. 

As  an  example  of  the  use  of  the  Buckingham  n  -theorem,  consider 
the  motion  of  a  body  with  acceleration,  a,  and  initial  velocity,  v.  The 
body  is  observed  to  move  a  distance,  s,  in  time,  t,  and  this  motion  is 
described,  of  course,  by  the  formula 

1  2 

<|>(s,  v,  a,  t)  =  s  -  vt  -  —  at  =0 

But  suppose,  for  the  sake  of  argument,  that  this  formula  is  not  known, 
and  that  the  object  of  an  experiment  is  to  discover  the  functional  form 
of  4>.  This  can  be  done  by  varying  v  and  a  in  some  systematic  manner 
and  by  measuring  s  as  a  function  of  t.  Now,  there  are  four  measurable 
quantities,  s,  v,  a,  and  t,  and  two  fundamental  dimensions,  length,  L., 
and  time,  T,  and,  according  to  the  Buckingham  fl  -  theorem,  there  are 
two  dimensionless  products,  and  ir^.  These  are 


sa 


and  there  must  be  a  dimensionless  function 


0 


which  describes  the  motion  of  the  body  in  question.  It  is  obvious  that  a 
function  of  two  parameters  is  simpler  than  a  function  of  four  parameters, 
and  so  it  is  easier  to  fit  if*  to  the  data  than  to  fit  <j>.  It  is  found  that 


7 


That  4)  is  equivalent  to  <j>  may  be  proven  by  multiplying  «j>  through  by  a/v  . 
In  general,  The  Buckingham  ri-  theorem  is  used  to  reduce  the  number  of 
independent  parameters  and  thereby  simplify  formulae  representing 
physical  laws. 

In  the  present  application,  the  measurable  quantities  which  describe 
the  motion  of  a  projectile  through  the  air  and  their  units  are  as  follows: 


Components  of  projectile  vilocity 

u 

L/T 

Components  of  projectile  angular 
velocity 

(0 

T"1 

Air  density 

P 

M/L3 

Viscosity  of  the  air 

^a 

M/LT 

Speed  of  sound  in  air 

vs 

L/T 

Projectile  diameter 

d 

L 

Force  on  the  projectile 

F 

ml/t 

There  is  a  total  of  n  =  11  measurable  quantities  (three  each  for  u  and  «), 
a  total  of  m  =  3  fundamental  dimensions  (M,  L,  and  T),  and  it  follows 
that  there  should  be  n  -  m  =  8  unique  dimensionless  products.  With  u 
instead  of  u^  they  are  as  follows: 


F 

Dimensionless  force 

,2  2 

pd  u 

s 

II 

Mach  number 

s 

r  =^i 

Reynolds  number 

u>.  d 

— —  Dimensionless  spin 

u 


u 


03 

u 


o),d 

and - 

u 


Components  of  cross  velocity 


Components  of  cross  angular 
velocity 


Thus,  there  is  a  dimensionless  relation  which  is  a  function  of  these 

2  2 

dimensionless  products.  If  this  relation  is  solved  for  F/pd^u  ,  an 
equation  of  the  form 


F 


,2  2„ 
pd  u  K 


is  obtained  in  which  the  function  K  is  dimensionless. 


Although  it  was  stated  that  the  components  of  u  are  among  the 
measurable  quantities,  u  is  used  herein  instead  of  u^.  The  principal 
references  for  this  development  are  Refs,  i,  3,  and  10.  In  Refs.  1  and 
10,  Uj  is  used,  whereas  u  is  used  in  Ref.  3.  The  arguments  of  dimen¬ 
sional  analysis  are  in  no  way  affected  b/  this  choice,  and  while  the  use 
of  Uj  is  advantageous  for  some  purposes,  the  use  of  u  is  more  common. 


This  development  ignores  the  fact  that  the  aerodynamic  force  can 
be  resolved  in  three  orthogonal  directions.  There  are  actually  three 
relations  such  as  the  one  above,  and  there  are  also  three  equations  for 
the  aerodynamic  moment  of  the  form 


G  =  pd3u2K> 


to^d  w^d 
u  *  u 


in  which  K  and  K'  are  known  as  aerodynamic  or  ballistic  coefficients. 

It  will  be  shown  in  the  following  sections  that  the  components  of  F  and  G 
can  be  represented  as  sums  of  expressions  similar  to  those  above  con¬ 
taining  dimensionless  functions  which  are  also  known  as  aerodynamic 
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coefficients.  This  will  follow  from  arguments  applicable  to  projectile 
symmetry. 

3.  Rotational  Symmetry 

In  the  preceding  subsection,  dimensional  analysis  was  employed 
to  ascertain  information  as  to  the  nature  of  the  functional  form  of  the 
aerodynamic  forces  and  moments  which  act  upon  a  projectile  in  free 
flight.  In  this  subsection  and  the  next,  more  information  of  the  same 
type  will  be  obtained  by  examination  of  the  properties  of  symmetry  of 
a  projectile.  Rotational  symmetry  will  be  treated  in  this  subsection 
and  mirror  symmetry  in  the  next  (Refs.  1,  3,  8,  10,  22,  and  23). 

In  vague  terms,  a  projectile  is  described  as  being  round.  What  is 
meant,  of  course,  is  tha.  the  appearance  of  the  projectile  is  unchanged 
by  a  rotation  about  its  longitudinal  axis.  It  follows  under  rather  general 
assumptions  as  to  surface  finish,  etc.  that  aerodynamic  forces  and 
moments  are  also  unchanged  by  such  a  rotation.  A  projectile  with  these 
properties  is  said  to  possess  rotational  symmetry  and  these  properties 
may  be  exploited  to  gain  additional  information  as  to  the  nature  of  the 
aerodynamic  forces  ana  moments. 

It  is  assumed  that  the  projectile  in  question  has  an  angle,  0,  of 
rotational  symmetry,  where  0  <  0  <  it,  and  that  aerodynamic  forces  and 
moments  are  unaffected  by  rotating  the  projectile  through  the  angle  0. 

It  is  noted  in  passing  that  the  analysis  to  follow  is  applicable  not  only  to 
shell,  but  to  missiles  with  three  or  more  symmetrically  located  fins. 

It  is  necessary  to  define  a  coordinate  system  and  parameters 
needed  in  the  development.  Aerodynamic  forces  and  moments  are  de¬ 
fined  with  respect  to  a  right-handed,  rectangular,  missile-fixed  x^, 
x^,  x^  coordinate  system  with  the  origin  at  the  center  of  mass  of  the 
projectile,  and  with  the  x^  axis  directed  toward  the  nose  along  the 
longitudinal  axis.  The  and  x^  axis  directions  are  fixed  in  the  projec¬ 
tile.  Projectile  velocity,  angular  velocity,  aerodynamic  force,  and 
aerodynamic  moment  are  represented  by  symbols  u,  w,  F,  and  G,  re¬ 
spectively,  and  components  of  these  vectors  along  a  designated  axis  are 
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denoted  by  the  subscript  denoting  that  axis  (e.  g. ,  u^  is  the  component 
of  u  along  x^). 

In  the  development  of  expressions  for  aerodynamic  forces  and 
moments,  it  is  convenient  to  use  a  complex  number  notation.  The  com¬ 
plex  number  notation  provides  a  convenient  way  of  designating  directions 
of  aerodynamic  forces  and  moments,  and  it  is  useful  in  exploiting  rota¬ 
tional  and  mirror  symmetry.  The  x^»  x^  plane  is  taken  to  be  the  com¬ 
plex  plane,  with  x^  the  imaginary  axis.  Let 


U ~  *  1 u. 

X  =  £  * 

u 

(Cross  velocity) 

(1) 

(u>,  +  i  go~)  d 

Ls  J 

(Cross  angular  velocity) 

(2) 

f=F2+iF3 

(Cross  force) 

(3) 

g  =  G2  +  iG3 

(Cross  moment) 

(4) 

2  2 


From  the  results  of  the  preceding  section,  it  follows  that  F  /pd'u 


,2  2 


Y 


and  f/pd'u  are  functions  of  the  Mach  number,  M,  and  Reynolds  number, 
Rg,  the  dimensionless  spin,  v,  and  also  of  X,  X.,  p,  and  p ,  where  a  bar 
over  a  symbol  denotes  the  complex  conjugate.  It  is  observed  in  passing 
that  the  notation  employed  here  differs  slightly  from  that  of  Refs.  1,  8, 


2  2 


and  10  in  which  u  replaces  u  in  Eqs.  (1)  and  (2)  and  in  the  ratios  F  /pd  u 

2  2  * 

and  f/pd  u  .  This  in  no  way  affects  the  analysis  to  follow  and  leads  to 

the  definition  of  aerodynamic  coefficients  consistent  with  current  usage. 


At  this  point,  an  approximation  is  made.  Since  usually  u  and  go 
are  almost  parallel  to  x^  (the  nose  of  the  projectile  nearly  points  along 
the  direction  of  motion),  it  follows  that 


|  X  |  =  |  X  |  «  1 ,  |p  j  =  |  p  |  «  1 
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V»>* in? *3 A.Jtjvpwf *»=n 


and  it  is  reasonable  to  expand  the  aerodynamic  forces  in  Taylor's  series 
expansions  in  X,  X,  p,  and  p  containing  only  constant  and  first  order 
terms,  i.  e.  , 

X  +  (jl  +  bjX  +  b2  jjl  (5) 

X  +  c2(i  +  di  X  +  d2il  (6) 

Since  F^  is  real,  the  right  side  of  the  first  equation  is  real.  The  coef¬ 
ficients  a  ,  b  ,  etc.,  are  functions  of  M,  R  ,  and  v,  and  for  the  present 
are  assumed  to  be  independent  of  X,  X,  p,  andp.  A  different  interpreta¬ 
tion  will  be  given  later  and  this  requirement  will  be  relaxed,  somewhat. 

Forces  F^  and  f  result  from  the  cross  velocity,  X,  and  the  cross 
angular  velocity,  p,  and  any  changes  in  X  and  p  will  cause  changes  in 
F^  and  f.  Rotational  symmetry  will  be  exploited  by  rotating  X  and  p, 
without  change  in  magnitude,  through  the  angle  of  rotational  symmetry, 

9.  This  will  not  change  F^,  but  it  will  induce  a  similar  rotation  of  f, 
and  it  is  equivalent  to  rotating  the  projectile  through  an  angle,  -0. 

If  £  is  a  vector  in  the  complex  x  ,  x_  plane,  it  can  bt  written  as 

C  =  r  e1Y  (7) 

If  £  is  subjected  to  a  positive  rotation  through  the  angle  0  about  the  x^ 
axis,  it  is  transformed  into 

i(Y+0) 

=  re 


,2  2  "  ai  *  bl 
pd  u 


,2  2  a2  +  C1 
pd  u 


or 


(8) 
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Hence,  rotation  of  X  and  p  through  9  results  in  the  following  transforma¬ 
tions: 


CD 

•«H 

o 

H 

(9) 

p*  =  pe10 

(10) 

f*  =  fel9 

(ID 

and 

F^  =  F1  (12) 

where  a  prime  denotes  parameters  associated  with  the  new  positions  of 
X  and  p.  The  truth  of  the  equation  for  f'  follows  from  rotational  sym¬ 
metry.  The  orientation  of  f'  with  respect  to  X1  and  p'  is  the  same  as 
that  of  f  with  respect  to  X  and  p.  It  follows  that  Eqs.  (5)  and  (6)  must  be 
satisfied  by  F|,  f',  X',  andp',  and  so 


—tt  =  ai +  V +  V +  V +  b2^ ' 

pd  u 


2T2  =  a2  +  ClV  +  C2^'+dlX'  +  V 


(13) 

(14) 


d  u 


Substitution  of  Eqs.  (9)  through  (12)  for  X',  p f',  and  F^  shows  that 


pd  u 


1  .  i 9  .  ,  i9  ,  —  r  -i  9  .  r-  —  -i 9 

2~z  ~  aj  rbjXe  +  b^p  e  +  b^Xe  +  b^p  e 


fe10  ,  .  i9  ,  i0  J  ,  r  -i9  .  ,  —  -i9 

- 5— y  =  a  +c  Xe  +  c_p  e  +  d  Xe  +  dpe 

pd  u 


and  comparison  of  these  relations  with  Eqs.  (5)  and  (6)  reveals  that 


bi(‘  -eie)  =  0 
b2(l  -a»)=  0 


(‘  -  = 


=  0 


d.  (i  -  e"2l0)=  0 


1 


-z(*  -*-2>  » 


Since,  by  definition,  0  <  8  <  tt,  it  follows  that 


'l  =  b2  '  a2  *  dl  =  d2  =  ° 


and 


.2  2 

pd  u 

f 

rr 


pd  u 


=  a. 


=  c]  X  +  c2  M- 


(15) 

(16) 


Completely  analogous  arguments  apply  to  the  aerodynamic  moments  and 
hence 


GI 

TT 


pd  u 


(»d 


=  e. 


-  c^K  +  p. 


(17) 

(18) 
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The  expressions  above  were  derived  for  a  system  of  coordinates 

fixed  in  the  projectile.  Actually,  they  are  valid  if  the  x^  and  x^  axes 

are  allowed  to  rotate  with  respect  to  the  projectile,  or  vice  versa.  To 

4  >!< 

prove  this  statement,  suppose  that  anx^,  x^,  x^  system  is  defined  such 

*  4c 

that  x.  coincides  with  x  but  the  angle  between  x  and  x  is  <{>,  where  <j> 

1  1  u  u 

may  be  a  function  of  time  (Fig.  1).  In  the  x^,  x^  system,  the  vector  £ 
may  be  written 


r,  =  r  e 


♦  * 


whereas  in  the  x^,  x^  system,  its  representation  is 

=  r  ei(v+4>) 


so  the  transformation  equation  between  the  two  systems  is 


i*  =  te1* 


Then 


f*  =  fe 


* 

g  =  ge 


\  =  \e 


* 

^  =  ^e 


i<t> 

icj, 

i4> 

i4> 


Obviously  F  and  G  are  independent  of  coordinate  systems,  and  multi- 

*  •  i  X 

plying  Eqs.  (16)  and  (18)  through  by  e  ™  yields 

f*  * 

—nr  =  cix  +  V 

pd  u 


-3  2 
pd  u 
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Figure  1 


The  Relation  Between  the  xltx2,x3  System  and  the 


x?,  xt,  x"  System;  x  =  x  Is  out  of  the  Page 


jjg  jjg  jjg 

Thus,  the  force  and  moment  equations  are  valid  in  the  x^,  x^,  system. 
This  is  an  important  simplification. 

4.  Mi;  ror  Symmetry 

One  might  suppose  that  if  the  image  of  a  projectile  as  seen  in  a 
mirror  were  real,  its  motion  would  be  governed  by  the  same  laws  of 
physics  as  those  which  determine  the  motion  of  the  projectile  itself.  Jf, 
iu  addition,  the  projectile  possesses  mirror  symmetry,  as  will  be  de¬ 
fined,  it  should  be  possible  to  infer  the  equations  of  motion  of  the  "image 
projectile"  in  terms  of  those  of  the  projectile  itself.  By  this  process, 
additional  information  can  be  gained  concerning  the  aerodynamic  forces 
and  moments  (Refs.  1  and  18). 

The  definition  of  mirror  symmetry  is  as  follows  (Ref.  1).  A  pro¬ 
jectile  possesses  mirror  symmetry  if  there  is  a  plane,  which  contains 
the  longitudinal  axis  of  the  projectile,  such  that,  if  each  point  of  the 
projectile  is  moved  to  the  point  on  the  opposite  side  of,  <_  ad  equidistant 
from,  the  plane,  then  the  projectile  exactly  covers  itself. 

It  is  assumed  that  the  projectile  in  question  possesses  a  plane  of 
mirror  symmetry,  in  which  case  the  axis  must  be  in  tha:  plane,  and 
it  is  also  convenient  to  define  the  x  axis  as  being  in  tha1;  plane. 

Consider  the  following  "thought  experiment"  (Fig.  2)  in  which  a 
projectile  is  fired  past  a  mirror.  It  is  arranged  so  that  the  projectile 
plane  of  mirror  symmetry  is  parallel  to  the  mirror  at  a  given  instant 
of  time  when  observations  of  the  projectile  and  its  image  are  made.  One 
can  choose  to  imagine  the  surface  of  the  mirror  and  the  plane  of  mirror 
symmetry  are  coincident  but  this  is  not  actually  necessary.  It  will  be 
supposed  that  u  and  u  represent  the  motion  of  the  projectile  whereas 
u'  and  u)'  represent  that  of  the  image.  Similarly,  F  and  G  represent  the 
aerodynamic  force  and  moment,  respectively,  acting  on  the  projectile, 
whereas  F1  and  G'are  associated  with  the  image. 

Since  x^  and  x^  lie  in  the  plane  of  mirror  symmetry  and  x^  is  per¬ 
pendicular  to  it,  it  follows  that  x!  and  x’  are  parallel  to  x  and  x  , 

k  L*  i  L> 
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respectively,  and  that  x^  is  measured  in  the  direction  opposite  to  x^i 
where  x^,  x^,  x^  is  the  coordinate  system  associated  with  the  image. 

Expressions  for  the  components  of  F',  G1,  u',  and  w'  in  the  x^, 
x^»  x^  coordinate  system  are  required  in  terms  of  the  components  of 

F,  G,  u,  and  w.  The  linear  motions  of  the  projectile  and  its  image 
are  the  same  in  the  1-  and  2 -directions,  and  so  are  the  forces  which 
cause  them;  in  the  3-direction,  the  motions  and  forces  are  opposite. 
Hence 


On  the  other  hand,  the  angular  motions  and  moments  are  opposite  in 
the  1-  and  2-directions  and  are  the  same  in  the  3-direction,  so 


"1  =  ’W1 


G1  =  -G1 


U2  ~w2 


G2  "  -G2 


w3  =  “3 


G3  =  G3 


It  follows  that 


=  F. 


£■  =  F'2tiF‘  =F2-iF3=  f 


-G, 


G2+*G2  =  ”G2  +  ’,G3  = 


-g 


(19) 

(20) 
(21) 
(22) 
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X' 


(22) 


u2+iu3  _U2~iu3 
u  u 


p’ 


Expressions  developed  in  the  preceding  section  for  the  aerodynamic 
forces  and  moments  contain  coefficients  c^  c^,  etc.,  which  are  functions 
of  M  R  ,  and  v.  The  values  of  M  and  R  are  unaffected  by  reflection 
through  the  plane  of  mirror  symmetry,  whereas  v  changes  sign.  Accord¬ 
ingly,  for  the  projectile,  we  write 


(co'2  +  iu!j)d  (-0J2+iu.j)d 


gj'j  d 

u 


Wjd 

u 


u 


=  -v 


=  "P 


(24) 

(25) 


2  2  "  al 
pd  u 


- Y  =  Cj(v)  X  +  c 2(v)p. 

pd  u  ' 


pd  u 


3  2"  el*v* 


(26) 

(27) 


(28) 


— f-F=  c3(v)X+  c4(v)p  (29) 

pd  u 


and  for  its  image, 


Fi 

r  =  ai(v,) 

pd  u 


f' 

TT 


pd  u 


(30) 


-  Cl(v«)X'  +  c2(vV 


(31) 


r 


pd  u 


=  ej(v') 


pd  u 


j—  =  c3{v*)X.'  +  c4(v')|i' 


Substitution  of  Eqs.  (19)  through  (25)  into  Eqs.  (30)  through  (33)  yields 


j  t 

pd  u 


=  a.(-v) 


A  2  =  cl("w)X‘  “  c2<“vhl 

pd  u 


j  J  * 

pd  u 


=  e1('v) 


pd  a 


2  =  C3(~v)^  “  c(jV-,/)p’ 


Comparison  of  this  set  with  Eqs.  (26)  through  (29)  leads  to  the  con¬ 
clusion  that 


ai(-v)  =  a^v) 

e^(~v)  =  _e  j (v) 

Cjt-v)  =  C^v) 


c?(-v)  =  -cJv) 


c,(-v)  =  -C,{v) 


c^'v) =  ca(v) 
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Evidently  a^,  the  real  parts  of  c^  and  c^,  and  the  imaginary  parts  of 

and  are  even  functions  of  v,  whereas  e^,  'he  imaginary  parts  of 

c,  and  c,t  and  the  real  parts  of  c,  and  c_  are  odd  functions  of  v. 
14  2  3 

5.  The  Ballistic  K  Notation 

In  the  last  three  subsections,  arguments  of  dimensional  analysis 
and  symmetry  were  used  to  gain  insight  into  the  functional  form  of  the 
aerodynamic  foices  and  moments  a.cting  on  a  projectile  in  flight.  It 
was  found  that  the  aerodynamic  forces  and  moments  can  be  represented 
by  equations  of  the  form 


where 


a2  2 

F,  =  pd  u  a. 
l  r  i 

F 2  +  iF ^  =  pd2u2  Jc^K  +  c^pj 

-  ,3  2 

Gj  =  pd  u  e^ 

+  iG^  =  pd^u^jc^X  + 


X 

v 


+  iu3 

U 

(u>2  +  iw^Jd 
u 


Wjd 

u 


and  where  the  dimensionless  coefficients  a  ,  c, ,  c_,  etc.,  are  functions 

l  i  £ 

of  M,  Re,  and  v ,  and  satisfy  Eqs.  (38)  through  (43)  of  the  last  subsection. 
The  only  restriction  imposed  upon  the  validity  of  these  equations  is  that 
the  first-order  Taylor's  series  expansions  in  X,  X,  p,  and  p  are  adequate, 
and  higher  order  terms  are  negligible.  For  the  moment,  the  truth  of 
this  hypothesis  is  assumed  so  that  the  K  notation  may  be  introduced. 


In  the  equations  above,  it  is  convenient  to  replace  the  coefficients 

a  ,  c  ,  c - ,  etc.,  by  a  new  set  of  coefficients  in  which  each  term  is 
112 


divided  into  its  real  and  imaginary  parts,  and  in  which  all  coefficients 
are  even  functions  of  v(e.  g.  ,  if  a(v)  is  odd  and  vb(v)  =  a(v),  then  b(v) 
is  even).  The  equations  above  are  written  as 


F  j  =  -pd2u2KDA  (44) 

F ^  +  iF 3  =  pd2u2[|  +  ivKjJx  +|  vKxp  +  iKgjp)  (45) 

Gj  =  -pd3u2vKA  (46) 


G2  +  iG3  =  Pd3u2[|-vKT-iKM|\  +(-Kh  +  ivKXT]p]  (47) 

where  the  K's  are  aerodynamic  coefficients  and  are  functions  of  M,  Re, 
and  v  and  are  even  in  v .  The  K  notation  is  known  as  the  ballistic 
notation.  Signs  were  chosen  so  all  the  K's  are  positive  under  normal 
conditions,  but  experience  has  shown  that  Kj  is  usually  negative. 

Separation  of  real  and  imaginary  parts  yields  the  following 
equations: 


2  2 

F2  =  -pd  -  pd  uvKFu3 

3  3 

+  pd  uvKxpti>2  -  pd  uKgUXj 

Fj  =  -pd2uKNu3+pd2uvKFu2 

3  3 

+pd  uvKX(;,oj3  +pd  uKg<*>2 

3  3 

G2  =  -pd  uvKtu2  +  pd  uKMu3 

4  4 

-pd  uKHw2  -  pd  uvKXTa>3 

3  3 

G3  =  -pd  uvKtu3  -  pd  uKmu2 

4  4 

-pd  uKpj(*>3  +pd  uK^^.w2 


(48) 


(49) 


(50) 


(51) 


23 


6.  Physical  Identification  of  Forces  and  Moments 

Each  term  in  Eqs.  (44)  through  (47)  bears  a  name  and  a  physical 

interpretation.  Nomenclature  is  given  in  Table  I.  The  axial  drag, 

2  2 

-pd  u  is  obviously  the  component  along  the  projectile  axis  of 

symmetry  of  the  aerodynamic  motion-retarding  force,  whereas  the 

3  2 

spin  decelerating  moment,  -pd  u  vK^,  is  the  frictional  moment  which 
arises  from  the  spin,  and  which  tends  to  damp 

The  normal  force,  -pd2u2KNX,  is  perpendicular  to  the  projectile 
longitudinal  axis  (Fig.  3),  and  is  in  the  direction  opposite  to  the  cross 
velocity,  X,  From  the  definition  of  X  it  follows  that 

X  X  =  sin2  6 

where  6  is  the  yaw  angle,  and  Xu  is  the  component  of  projectile  velocity 
normal  to  the  projectile  longitudinal  axis.  The  plane  of  yaw  is  the 
plane  containing  u,  Xu,  and  the  projectile  longitudinal  axis,  so  the 
normal  force  is  in  the  plane  of  yaw. 

R  is  the  vector  sum  of  the  normal  force  and  the  axial  drag  (Fig.  3). 
Resolution  of  R  into  components  perpendicular  and  parallel  to  u  yields 
the  lift,  L,  and  drag,  D,  respectively,  where 

L  =  pd2u2KL  sin  6  (52) 

D  =  -pd2u2KD  (53) 

These  two  forces  are  equivalent  to  the  normal  force  and  the  axial  drag, 
and  are  more  convenient  for  calculational  purposes.  It  follows  from 
the  definitions  that 


KL  =  Rn  c°s  6  -  Kda 

(54) 

KD  =  KN  sin2  6  +  Kda  cos  6 

(55) 
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Table  I 

Ballistic  Nomenclature 


,2  2IT 

-Pd  u  Kda 

Axial  drag 

-pd^2!^^ 

Normal  force 

ipd2u2vKp\ 

Magnus  force 

•  a1  2tr 

lpd  u  Kgp 

Cross  force  due  to  cross  spin 

,2  2  __ 
pd  u  vKXFp 

Magnus  cross  force  due  to 
cross  spin 

,3  2  __ 

-pd  u  vKa 

Spin-decelerating  moment 

.3  2  . 

-pd  u  vK,pX 

Magnus  moment 

-ipd3u2KMX 

Overturning  moment 

-pd3u2Kj^p 

Damping  moment 

ipd“’u2vKx^,pi 

Magnus  cross  torque  due  to 
cross  spin 
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Normal  Force 


u 


Figure  3 

Forces  in  the  Plane  of  Yaw 


The  lift  and  drag  are  usually  used  in  trajectory  computations  rather 
than  the  normal  force  and  the  axial  drag. 

3  2 

The  overturning  moment,  -ipd  u  K^j\,  is  the  torque  which  re¬ 
sults  from  the  action  of  the  normal  force, and  its  direction  is  that  of 
-i\.  It  is  perpendicular  to  the  plane  of  yaw  and  is  out  of  the  page  in 
Fig.  3  (not  shown).  For  a  spin- stabilized  shell,  the  center  of  pressure 
of  the  normal  force  is  usually  ahead  (toward  the  nose)  of  the  center  of 
mass,  and  the  overturning  moment  usually  acts  to  increase  the  yaw 
angle,  6. 

2  2 

The  Magnus  force,  ipd  u  v  KpX,  is  perpendicular  to  the  plane  of 
yaw  and  is  into  the  page  in  Fig.  3  (not  shown).  It  results  from  the  pro¬ 
jectile  spin  about  the  longitudinal  axis. 

3  2 

The  Magnus  moment,  -pd  u  v  K,pX,  is  the  torque  which  results 
from  the  action  of  the  Magnus  force.  Its  direction  is  the  same  as  that 
of  the  normal  force  provided  K^,  is  positive. 

2  2 

The  cross  force  due  to  cross  spin,  ipd  u  Kgp,  results  from  the 
component  of  cross  velocity,  p,  as  the  name  implies.  Its  direction  is 
that  of  ip,  and  is  perpendicular  to  p  and  the  projectile  longitudinal 
axis.  The  cross  force  due  to  cross  spin  is  physically  small  and  is 
assumed  to  have  negligible  influence  on  projectile  motion.  Its  presence 
will  be  ignored  in  trajectory  calculations. 

3  2 

The  damping  moment,  -pd  u  K^p,  *s  t*ie  torque  which  results 
from  the  action  of  the  cross  force  due  to  cross  spin.  Its  direction  is 
opposite  to  that  of  p,  and  its  action  is  to  damp  the  cross  angular  velocity. 
It  has  an  important  influence  upon  projectile  motion. 

2  2 

The  Magnus  cross  force  due  to  cross  spin,  pd  u  vK^pp,  has 
negligible  influence  upon  projectile  motion  and  is  ignored. 

3  2 

Likewise,  the  Magnus  cross  torque  due  to  cross  spin,  ipd  u^K^-j.u, 
which  is  the  moment  associated  with  the  Magnus  cross  force  due  to 
cross  spin,  is  negligible. 
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The  reason  for  retaining  the  cross  force  due  to  cross  spin,  lhe 
Magnus  cross  force  due  to  cross  spin,  and  the  Magnus  cross  torque 
due  to  cross  spin  in  the  analysis  is  for  logical  consistency.  If  these 
terms  wery  deleted,  a  shift  of  the  projectile  center  of  mass  would  lead 
to  a  logical  contradiction.  Otherwise,  these  terms  could  be  excluded 
from  further  consideration.  The  effect  of  a  shift  in  center-of-mass 
location  is  treated  in  the  following  subsection. 

7.  Center-of-Mass  Location 

Different  versions  of  a  particular  shell  (e.  g.  ,  HEI,  AP,  etc.) 
often  have  the  same  external  shape,  but  have  different  center-of-mass 
locations.  In  order  that  the  same  aerodynamic  data  can  be  used  for  all 
versions  of  a  projectile,  it  is  necessary  to  have  appropriate  formulae 
for  handling  a  shift  in  the  center  of  mass.  Such  relations  are  derived 
in  this  subsection. 

Aerodynamic  forces  and  moments  act,  respectively,  through  and 
about  the  center  of  mass  of  a  projectile,  and  it  follows  that  a  shift  in 
the  center-of-mass  location  will  produce  a  different  projectile  motion. 
On  the  other  hand,  it  is  clear  that  aerodynamic  r^essures  and  frictional 
forces  are  dependent  only  on  the  exterior  surface  of  the  projectile  and 
its  motion  relative  to  the  air,  and  do  not  depeid  upon  the  center-of- 
mass  location.  These  facts  are  used  to  derive  the  required  formulae. 

At  a  given  instant  of  time,  it  is  assumed  that  F,  G,  u,  and  u>  are 
known  for  a  given  projectile.  If  the  center  of  mass  is  moved  forward 
a  distance  r  along  the  longitudinal  axis  to  a  new  position,  and  if  primes 
denote  all  parameters  defined  in  a  new  xj, ,  x^,  coordinate  system, 
which  is  parallel  to  the  old  system  but  with  its  origin  at  the  new  center- 
of-mass  position,  then,  with  respect  to  the  new  system 

F»  =  F 

<5'  =  Cf  -  r  Xj  X  F 

=  x1G][  +  x2(G2  rF3)  +  x3(G3-rF2) 
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l*»"  >jM.-’y->.-> 


w'  =  0) 

u1  =u  +  (oX  ritj 

=  +  x2(u2  +  ru>3)  +  (u3  -  ra>2) 

Axial  components  of  F  and  G  are  unchanged,  so  and  are  in¬ 

dependent  of  center-of-mass  position. 

In  terms  of  complex  numbers,  the  transverse  components  be¬ 
come 


X.'  = 


u'  +  iu! 


u 


=  X.  -  ihfi 


(co'2  +  iw'3  )d 
u 


=  V- 


f.  =  f>2  +  if «  =  f 

g1  =  G'2  +  iG^  =  g  -ihdf 


where  h  =  r/d  (h  is  r  measured  in  calibers).  It  is  assumed  u'  =  u  since 
the  difference  is  negligible.  The  transverse  force  and  moment  equa¬ 
tions  for  the  new  system,  Eqs.  (16)  and  (18),  become 


f' 

T~Z 


pd  u 


=  c'jX.'  +  c'2fi' 


pd  u 


3  2 -  =  C3X'  +  c4^' 


Then 
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c’jfA.-ityi)  +  c^jx  - - 2~~2 -  =  c  \  +  c2(Jt 

pd  a 

c'3(\-ihp)  +  c^p  =  ~^— —  =  c  A  +c.p  -  ih(cA  +  c,p) 
pd  u  i  c 

and  it  follows  from  a  comparison  of  the  coefficients  of  X.  and  p  that 


Ci  =  C1 


C2  =  c2  +  ihcj 


C3  =  c3  '  ihcl 


c4  ~  c4  -ih(c 2-c ^ )  +  he 


Expressions  for  the  new  coefficients  in  terms  of  the  old  may  be  derived 
by  use  of  Eqs.  (45)  and  (47).  They  are: 


kda 

=  kda 

KX 

II 

“n 

=  kn 

km 

=  KM-hKN 

kf 

II 

t* 

*3 

K't 

=  KT-hKF 

Ks 

=  VhKN 

Kk 

=  Kn-MKg  +  km)  +  hZKN 

kxf 

=  KXF-hKF 

KXT 

=  KXT-h(KXF+KT)+h2KF 

8.  Weaknesses  of  the  Aerodynamic  Force  and  Moment  System 

The  aerodynamic  system  of  forces  and  moments  which  is  described 
herein  suffera  from  certain  weaknesses.  There  are  at  least  three  areas 
where  one  can  f.nd  fault,  namely  in  arguments  based  upon  projectile 
symmetry,  in  the  Taylor's  series  assumption,  and  the  projectile  past 
history.  These  will  be  discussed  in  turn. 


The  assumption  that  every  projectile  possesses  perfect  rotational 
and  mirror  symmetry  is  obviously  incorrect.  This  is  a  physical  im¬ 
possibility.  One  would  expect  slight  deviations  from  rotational  sym¬ 
metry,  and  rifling  marks  on  a  shell  weaken  the  assumption  of  mirror 
symmetry.  On  the  other  hand,  this  is  not  a  very  serious  problem  and, 
in  practice,  such  deviations  from  the  ideal  should  produce  negligible 
error  on  the  average.  These  effects  probably  give  rise  to  a  round-to- 
round  dispersion  which  may  even  be  desirable.  The  effect  of  slight  con¬ 
figurational  asymmetries  is  treated  in  Ref.  1  3. 

The  error  introduced  by  the  assumption  that  aerodynamic  forces 
and  moments  can  be  represented  by  a  first-order  Taylor's  series  is 
more  important.  The  Taylor's  series  assumption  implies  that  the 
aerodynamic  forces  and  monfents  are  continuous  functions  with  all  de¬ 
rivatives,  but  the  characteristics  of  fluid  flow  tend  to  discount  this 
assumption.  Instantaneous  transitions  in  flow  are  observed,  such  as 
a  breaking  away  of  the  flow  from  a  particular  location  near  the  nose  of 
the  projectile  as  the  yaw  angle  increases.  Tnis  suggests  a  discontinuous 
nature  for  aerodynamic  forces  and  moments.  Also,  the  assumption 
that  forces  and  moments  are  linear  in  X  is  violated  in  practice. 

There  is,  fortunately,  a  way  around  these  difficulties.  A  close 
examination  of  the  arguments  of  symmetry  reveals  that  there  is  no 
contradiction  if  the  aerodynamic  coefficients  (the  K's)  are  allowed  to 
vary  with  products  of  the  form  XX,  Xp,  Xp,  and  pp  .  These  products 
are  invariant  under  rotation  of  coordinates  about  the  projectile  longi¬ 
tudinal  axis  and  so  are  the  aerodynamic  coefficients.  Measurements, 
in  fact,  show  a  variation  of  the  aerodynamic  coefficients  with 
XX  =  sin  6,  but  no  significant  variation  with  the  three  other  products. 

It  is,  accordingly,  convenient  to  drop  the  Taylor's  series  assumption, 
but  to  retain  the  equations  derived  from  it  as  a  convenient  notation  in 
which  the  aerodynamic  coefficients  are  allowed  to  vary  with  sin  6. 
Should  the  need  arise,  additional  correction  terms  which  possess  the 
appropriate  symmetry  may  be  added. 
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The  assumption  that  aerodynamic  forces  and  moments  depend  only 
on  velocity,  u,  and  angular  velocity,  u,  is  subject  to  question.  A  pro¬ 
jectile  may  conceivably  reach  particular  values  of  u  and  «  by  means 
of  a  different  history  of  motion.  In  this  instance,  the  flow  pattern  of 
the  air  around  the  projectile  and  the  resulting  aerodynamic  forces  and 
moments  may  be  different. 

Aerodynamic  forces  and  moments  are  dependent  upon  the  pattern 
of  flow  of  the  air  mass  past  a  projectile  whereas  the  pattern  of  flow  is 
in  turn  a  function  of  the  history  of  motion  of  the  projectile.  Presumably, 
the  past  history  of  projectile  motion  can  be  reconstructed  from  knowl¬ 
edge  of  u  and  u)  and  all  time  derivatives  of  u  and  w  at  a  given  instant. 

If  this  is  so,  it  follows  that  the  pattern  of  flow,  and  hence  the  aerody¬ 
namic  forces  and  moments,  are  functions  of  u  and  to  and  all  time  deriva¬ 
tives  of  u  and  oj  .  Projectile  acceleration  and  angular  acceleration  are 
included  in  the  generalized  force  system  developed  in  Ref.  10.  This 
system  contains  the  following  torms: 
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J2  is  the  component  of  the  angular  velocity  of  the  x  ,x_,x_  coordinate 

X  i  b  J 

system  in  the  x^  direction,  and  new  ballistic  coefficients  bear  the  addi¬ 
tional  subscript  A  for  acceleration.  These  equations  differ  from  those 
of  Ref.  10  in  that  u  replaces  u^. 


Reasonable  assumptions  as  to  the  size  of  the  acceleration  coef¬ 
ficients  (Ref.  10)  lead  to  the  conclusion  that  they  were  all  negligible 
except  for  K 


MA* 


In  fact,  an  approximate  solution  of  the  projectile  equa¬ 
tions  of  motion  reveals  that  free -flight  test  range  measurements  for  Kpj 
actually  yield  rather  than  K^.  There  is  no  way  to  separate 

K TJ  and  K.  ,A  in  such  measurements. 

H  MA 


In  conclusion,  it  should  be  remarked  that  these  weaknesses  in  the 
development  do  not  necessarily  imply  that  ballistics  calculations  are  of 
insufficient  accuracy  to  be  useful.  The  point  is,  rather,  that  ballistics 
calculations  may  not  be  as  accurate  as  one  might  expect  if  one  is  un¬ 
aware  of  the  weaknesses. 


9.  The  Aerodynamic  Notation 

As  an  aid  to  the  free  exchange  of  information  between  the  aero- 
dynamicist  and  the  ballistician,  aerodynamic  notation  has  been  adapted 
to  projectiles  with  rotational  and  mirror  symmetry  (Refs.  18,  19,  and 
20).  The  orthogonal  right-handed  coordinate  system  usually  employed 
by  the  aerodynamicist  is  a  system  fixed  in  the  aircraft  or  missile  with 
the  X  axis  along  a  principal  axis  of  inertia  (the  longitudinal  axis)  and 
the  Y  axis  parallel  to  the  span  of  the  principal  lifting  surface  (out  the 
right  wing).  Components  of  the  aerodynamic  force  and  moment  in  this 
system  are 

L  =  Ip  V2SdC„ 

2 r  £ 

M  =  Ip  VZSdC 
2 r  m 

N  =  |PV2SdCn 


X  .  |Pv2scx 
V  =  |pv2scy 
Z  =  Ip  V2scz 
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where  X,  Y,  and  Z  are  components  of  force,  and  L,  M,  and  N  are 
components  of  moment.  Y  is  the  projectile  speed,  S  is  a  reference 
area,  and  d  is  a  reference  length.  Only  one  reference  length  is  used 
for  missiles,  contrary  to  standard  usage  for  aircraft.  Iu  the  present 
treatment,  the  reference  length,  d,  is  the  projectile  diameter  and  the 
reference  area  is 

S  .  Jd2 

Standard  notation  for  components  of  velocity  and  angular  velocity  in  the 
X,  Y,  and  Z  directions  are,  respectively,  u,  v,  and  w,  and  p,  q,  and  r. 
The  orientation  of  the  velocity  V  in  the  X,  Y,  Z  coordinate  system  is  de¬ 
fined  by  the  angle  of  attack,  a,  and  the  angle  of  sideslip,  (3,  where  a  is 
measured  about  the  Y  axis  and  (3  is  a  rotation  about  the  Z  axis.  For 
small  angles, 

w  .  v 

a  s  and  (3  s  y 

Nomenclature  for  the  drag  and  lift  forces  is 

Drag  =  i  pVZSCD 

Lift  =IpV2SCT  a 
Z  r  La 

and  that  for  the  spin  deceleration,  or  damping-in-roll,  is 

X  =!pV2Sd|^C<p 

The  normal  coefficients  Cv,  C„,  C  ,  and  C  are  assumed  to  be 

Yam  n 

functions  of  a,  (3,  p,  q,  r,  of  time  derivatives  a,  (3,  q,  and  r,  and  are 
expanded  in  first-order  Taylor's  series  in  terms  of  these  variables. 

The  linear  force  and  moment  coefficients  which  have  counterparts  in 
the  K  uotation  are  defined  by  the  following  relations: 


twh 


-*-**•***•■ 
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Terms  which  have  been  left  out  of  these  expansions  are  eliminated  by 
the  same  arguments  of  symmetry  used  in  the  development  of  the  K 
notation,  so  i*  is  convenient  to  ignore  them  here. 

The  notation  herein  differs  somewhat  from  that  of  the  references, 
and  indeed  the  notation  is  not  in  agj  .ement  between  references.  The 
nomenclature  herein  has  been  chosen  so  that  definition  of  coefficients, 
particularly  signs,  are  in  agreement  with  those  of  Arnold  Engineering 
Development  Center.  This  is  discussed  below. 

The  aerodynamic  coefficients  above  are  sometimes  referred  to 
as  aerodynamic  derivatives;  for  example,  a  notation  such  as 


3C  y 

CYpa  ~  9(pd/2V)  9a 


is  often  seen.  This  notation  is  strictly  correct  only  if  Cy  is  indeed  a 
linear  function  of  (pd/2V)  and  of  a. 

The  Y  axis  was  defined  as  being  parallel  to  the  span  of  the  princi¬ 
pal  lifting  surface  for  a  missile  which  has  such  a  surface.  For  a  mis¬ 
sile  with  an  angle,  0,  of  rotational  symmetry,  there  may  be  no  such 
unique  surface,  and  for  a  shell,  G  can  be  any  angle.  The  Y  axis  direc¬ 
tion  is  not  uniquely  defined. 

If  it  is  assumed  that  the  projectile  has  90-deg  rotational  symmetry 
(0  =  90°),  tiii symmetry  can  be  exploited  by  means  of  two  coordinate 
systems,  e.  g.  ,  If,  Y,  Z  and  X',  Y',  Z'  differing  by  a  rotation  of  90  deg 
about  the  X  axis.  It  follows  that  the  equations  above  are  valid  for  both 
systems  and  that 


F  ' 
Z 


q' 

r 1 


Fz 

V 

ii 

a 

N 

-F 

Y 

G  ' 

z 

II 

r 

a1 

=  -p 

q 

P* 

=  a 
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-..  re  relations  between  forces,  moments,  angles,  and  angular  rates  de¬ 
fined  in  the  two  coordinate  systems.  Utilization  of  these  relations  re¬ 
veals,  lor  example,  that 


C 


Za 


c 


YfJ 


Table  II  contains,  a  list  of  all  such  relations.  Because  of  the  lack  of 
unique  Y  and  Z  axis  directions,  subscripts  Y  and  Z  are  often  replaced 
by  the  symbol  N  (for  normal)  and  m  and  n  are  replaced  by  the  symbol 
M.  This  notation  is  also  included  in  Table  II,  and  is  called  the  aero- 
ballistic  notation  in  Ref.  19. 

1 0.  Comparison  of  Ballistic  and  Aerodynamic  Notations 

To  compare  ballistic  and  aerodynamic  nomenclature,  it  is  con¬ 
venient  to  use  complex  number  notation.  It  follows  from  Table  II  that 


f 


CY  +  lCZ 


r  j  Pd  c 

^Na  ‘  2V  Npa 


tj  (P+ia)  - 


pd  r  _  iC 

2V  ^Npq  Nq 


(q  +  ir) d 
2V 


L  .  pd  (P  +  ia)d  pd  iC  .  Atiill 

_  |CNa  “  X'cV  CNpa  2V  |^2V  CHpq  NqJ  4V 2 


(58) 


C  +  iC 
m  n 


f  Pd 

2V 


-j  »“ 

J<p+i»)  +  [' 


CMpa"iCMa  ^  +  ia)+  CMq_12V  CMpq  \~~ZV 


■  Pd  r  (q  +  ir)d 

2\ 


+  '—  CMpa"iCM4 


a  Z7  Mq  2V  Mpq  4y2 

•J  L  -J 


(59) 
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Table  II 


Equivalent  Aerodynamic  Coefficients 


CZa 

= 

CYP 

- 

CNa 

Cv 

Zq 

= 

"CYr 

= 

CNq 

GZa 

= 

CYP 

= 

Cxr* 

Na 

czq 

= 

'CYr 

= 

CNq 

= 

~CYpa 

- 

CNpa 

CZpr 

= 

CYpq 

- 

CNpq 

C  • 

zpp 

= 

"cyp& 

= 

G  Npa 

°ZPr 

-= 

CYpq 

= 

CNpq 

“ma 

= 

"CnP 

= 

CMa 

Cmq 

= 

Cnr 

= 

CMq 

C 

ma 

= 

"Cnp 

= 

C.,  • 
Ma 

C  • 

r  q 

= 

C  • 
nr 

= 

C  ♦ 
Mq 

Cmpp 

- 

C 

npa 

= 

CMpa 

C 

mpr 

= 

-C 

npq 

= 

C.. 

Mpq 

Cmp(3 

= 

C  • 
npa 

= 

C..  • 
Mpa 

C  • 
mpr 

= 

-C  • 
npq 

= 

c.,  • 

Mpq 
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Also,  it  is  evident  that,  for  small  angles, 

(3  +  ia  =  X. 


and 


(P  +  ia)d  ^  \d  X.' 

2V  =  2u  =  2 

where  the  prime  denotes  differentiation  with  respect  to  arc  length,  s, 
where 


Also, 


and 


s  .  f  > 

0 

(q  t  ir)d  <"2  +  iu3,d  .  U 
2V  2H  "  2 


(q  -r  ir)dZ  _  p 1  q.  p  V1  ^  p 

„,r2  4  2  V  =  4 


By  comparison  of  the  equations  above  with  the  corresponding  Eqs.  (56) 
and  (57),  involving  the  ballistic  notation,  a  one-to-one  correspondence 
can  be  made  between  the  aerodynamic  coefficients  and  the  ballistic 
coeffieients.  This  one-to-one  correspondence  is  shown  in  Table  III. 
Corresponding  coefficients  are  not  necessarily  equal;  if  the  C  notation 
is  used,  it  should  be  adopted  in  toto. 

1 1.  Arnold  Engineering  Development  Center  (AEDC)  Notation 

Assume  measurements  of  aerodyramic  forces  and  moments  are 
made  in  a  wind  tunnel  in  which  the  mot;on  of  the  projectile  model  is 
constrained  such  tha*-  there  is  no  sideslip  and  there  is  no  rotation  about 
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Table  HI 

Comparison  of  Ballistic  and  Aerodynamic  Notations 
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the  Z  axis,  i.e.  ,  such  that  (3  =  0,  r  =  0,  and  a  =  q.  (The  X  axis  is  taken 
to  be  at  an  angle  a  above  the  relative  wind  vector,  the  Y  axis  is  taken 
to  be  horizontal  and  out  the  side  of  the  wind  tunnel,  and  the  Z  axis  is 
directed  downward  and  at  an  angle  a  from  the  vertical.)  Under  these 
constraints  and  with  negligible  terms  excluded,  the  aerodynamic  equa¬ 
tions  become 


Pd 
Y  pa  2  V 


C„  ~sr—  a 


C  a 
Za 


C  a  +  (C  +  C  .) 
ma  mq  ma 


c  -gU 

npa  2V 


2V 


With  one  exception,  these  relations  define  the  notation  used  at  AEDC, 

where  CAT  is  used  in  place  of  C_  .  It  is  also  evident  that  C  and 
Ma  Za  mq 

C _ .  cannot  be  measured  separately  in  a  wind  tunnel.  A  similar  con- 

ma  1 

elusion  has  already  been  mentioned  with  regard  to  free-flight  test 
measurements  since 

K  -  K. , ,  =  — j- r —  ( C  +  .  \ 

H  Mi).  16  \  mq  ma  f 

(See  Subsection  8. ) 

It  was  mentioned  that  the  notation  herein  was  chosen  to  be  in 
agreement  with  that  of  AEDC.  The  following  comments  seem  appro¬ 
priate:  AEDC  reports  (Refs.  30  through  38)  available  to  the  author 
contain  data  for  CM  ,  C..  ,  C  ,  C  +  C  .,  and  C  .  Data  for 

the  other  coefficients  (besides  drag  and  damping -in-roll)  is  lacking 
and  for  them  no  definite  sign  convention  is  available.  These  missing 
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coefficients  are  considered  to  be  negligible,  however,  and  no  sign  con¬ 
vention  is  needed.  Numerical  values  of  C and  C  are  positive, 

iNa  ma 

whereas  C  +  C  .  is  negative.  Signs  of  Magnus  coefficients  Cv 

mq  ma  x  pa 

and  C  may  be  in  error  since  it  is  difficult  to  determine  their  sign 
npa  3 

conventions  from  AEDC  reports.  It  is  noted  that  AEDC  does  not 
strictly  adhere  to  these  definitions.  One  should  be  very  careful  to  ob¬ 
tain  correct  information  as  to  definitions  of  force  and  moment  coef¬ 


ficients  from  wind  tunnel  and  free -flight  test  range  personnel.  Their 
definitions  may  differ  somewhat  from  those  herein,  and  it  is  very  easy 
to  get  a  sign  wrong  or,  say,  be  off  by  a  factor  of  two.  Some  people  use 
pd/V  in  their  definitions  rather  than  pd/2V. 


Reference  39  (Advisory  Group  For  Aerospace  Research  and 
Development,  AGARDograph  121)  contains  methods  of  obtaining  aero¬ 
dynamic  data  by  use  of  wind  tunnels.  The  notation  herein  agrees  with 
that  defined  in  Ref.  39  except  for  which  is  negligible. 

The  notation  herein  differs  somewhat  from  that  of  Ref.  40 
(Ballistics  Research  Laboratory  Memorandum  Report  No.  2192).  The 
major  difference  is  in  the  use  of  pd/V  instead  of  pd/2V  for  the  dimen¬ 
sionless  spin. 


12.  Modeling  of  Aerodynamic  Data  for  Trajectory  Computations 

In  previous  subsections,  expressions  have  been  developed  for 
aerodynamic  forces  and  moments  in  terms  of  the  aerodynamic  coef¬ 
ficients,  but  nothing  has  been  said  about  the  form  of  these  c  .  ficients 


other  than  that  they  are  functions  of  M,  Rg,  v,  and  6.  Attempts  have 
been  made  to  derive  expressions  for  these  coefficients  from  fundamental 


theory,  but  in  practice,  it  is  necessary  to  measure  them  for  various 
combinations  of  M,  Rg,  v,  and  8.  Such  measurements  are  made  in 
wind  tunnels,  or  by  means  of  free-flight  tests  (Refs.  30  through  47). 
Data  thus  obtained  is  tabular  and  a  problem  arises  in  the  utilization  of 
this  data  in  trajectory  computation  on  a  digital  computer.  Means  of 
multi-dimensional  interpolation  between  data  points  is  required.  To 


ascertain  the  extent  of  the  problem,  it  is  advantageous  to  examine  the 
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tabular  data  and  determine  as  much  as  possible  about  the  functional 
form  of  the  aerodynamic  coefficients. 


12.  1  The  Functional  Form  of  the  Aerodynamic  Coefficients  -  As 
has  been  stated,  the  aerodynamic  coefficients  are  functions  of  Mach 
number,  M,  Reynolds  number,  Rg,  dimensionless  spin,  v,  and  the  yaw 
angle,  6.  Fortunately,  variations  with  R^  are  slight  and  usually  can  be 
ignored.  The  expression  for  R^  is 


R 

e 


pud 

^a 


MV  d 
s 


where  r|  is  the  kinetic  viscosity  and  V  is  the  speed  of  sound.  Since 
3>  S 

both  V  and  u  are  functions  of  altitude,  R  varies  with  altitude  and 
s  ra  e 

attempts  have  been  made  to  model  certain  aerodynamic  data  as  such 
(Ref.  48).  Unless  extreme  variations  in  altitude  along  a  trajectory  are 
anticipated,  however,  Reynolds  number  variations  may  be  ignored. 

This  may  be  regarded  as  the  best  policy  unless  experience  proves 
otherwise. 

The  effect  of  dimensionless  spin,  v,  upon  the  drag,  lift,  and  over¬ 
turning  moment  coefficients,  Kj-j,  K^,  and  K^,  respectively,  is  be¬ 
lieved  to  be  negligible.  That  there  is  a  slight  effect  is  demonstrated 
in  Ref.  7.  Until  proven  otherwise,  Kp,  K^,  and  are  assumed  to  be 
independent  of  v.  The  damping  moment  coefficient  is  also 
usually  modeled  as  such.  On  the  other  hand,  the  Magnus  coefficients, 
Kp,  and  Kp,  are  both  functions  of  v  and  should  be  modeled  as  such.  The 
spin  deceleration  coefficient,  K^,  is  presumably  also  a  function  of  v, 
but  is  hard  to  measure  and  little  data  may  be  available.  This  poses 
somewhat  of  a  dilemma  since  lack  of  knowledge  of  in  trajectory 
computation  implies  a  lack  of  knowledge  of  v.  Extensive  modeling  of 
Kp,  and  Kp  in  terms  of  v  is  then  of  questionable  value. 

All  coefficients  are  functions  of  M  and  6,  although  adequate  data 
for  modeling  may  be  hard  to  obtain.  This  is  particularly  true  of  Kp 
and  Kp  since  they  are  functions  of  three  variables,  M,  v,  and  5. 
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In  summary,  the  aerodynamic  coefficients  may  be  written  in 
terms  of  their  independent  variables  as  follows: 


kd 

=  kd 

(Mj 

5) 

kl 

'  kl 

(M, 

6) 

km 

=  Km(M, 

6) 

s  kh 

(M, 

6) 

kf 

=  kf 

(M, 

v, 

6) 

Kt 

=  Kt 

(M, 

v, 

6) 

ka 

=  ka 

(M, 

V, 

6) 

All  other  aerodynamic  forces  and  moments  are  considered  to  be  negli¬ 
gible,  and  are  excluded  from  consideration. 

12.2  Polynomial  Curve  Fits  -  Modeling  of  aerodynamic  data  for 
computer  use  can  be  accomplished  by  fitting  low -order  polynomials  in 
two  or  three  variables  to  the  data.  The  drag  coefficient  can  usually  be 
expressed  by  the  relation 

KD  =  Kd  (M)  1  +  Kd  (M)  sin2  5 

o  L  6 


in  which  Kp  (M)  is  the  zero-yaw  drag  and  the  term  in  the  square  brack - 


o 

ets  accounts  for  non-zero  yaw.  For  small-yaw  applications,  sin  6  is 

sometimes  replaced  with  6,  whereas  for  high-yaw  situations,  more 
4  6 

terms  in  sin  6,  sin  6,  etc.  can  be  included  if  they  are  needed.  (M) 

o 

and  Kn  ,(M)  may  be  fitted  in  sections  in  powers  of  M  or  1/M.  If  M  , 

6  1 


M  are  numbers  such  that 
n 


for  example, 


M^  <  M^  <  .  . .  <  M 


n 
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fffgwCEgffy? 


a.  + 

i 


b. 

1 

M 


provides  a  good  fit  to  K^  over  the  interval 

o 

M.  <  M  =s  M. , , 
i  l+l 


if  the  interval  is  not  too  large.  On  the  other  hand,  an  expression  valid 
for  high  Mach  numbers  has  been  found  useful.  It  is 


(a  +  bM)^ 
M2 


1 


This  can  be  obtained  by  curve  fitting  the  "Qfunctiou"  (Ref.  49). 


Q  =  /  1  +  M2Kd  =  a  +  bM 

N  O 

The  lift  coefficient,  K^,  may  be  handled  in  a  similar  fashion.  In 

the  interval  M.  <  M  <  M.  ,  ,  and  sin  6.  <  sin  6  <  sin  6.  ,  ,  one  can  write, 
l  l+l  l  l+l 

for  example 


KL  =  ^a0  +  aiM  +  *  *  *  +  ai  ) 

+  (a«  +  b,M  +  •••  +  b  Mm)  sin2  6 
'  0  1  m 

+  .  . . 

+  (cn  +  c,M+  •••  +  c  Nn)  sin2^  6 
1  0  i  n 

Similar  expressions  can  be  written  for  K.,  and  KIT. 

M  H 

Kp,  Krp,  and  are  functions  of  three  variables:  M,  v,  and 
sin26.  Polynomials  may  be  constructed  for  them  in  much  the  same 
manner. 


45 


2 

^ical  plots  of  Kq  vs  M  and  Kq  vs  sin  6  with  M  constant  are 
o 

shown  in  Fig.  4.  Plot*5  of  other  aerodynamic  coefficients  may  be 
found  in  the  literature.  See,  for  example,  Refs.  30  through  40.  Three- 
dimensional  plots  of  some  of  the  coefficients  are  shown  in  Ref.  50. 
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SECTION  III 


A  STUDY  OF  PROJECTILE  ANGULAR  MOTION 
1 .  General 

The  system  of  aerodynamic  forces  and  moments  has  been  sub¬ 
jected  to  a  comprehensive  examination  in  the  previous  section  and  it 
is  convenient  at  this  point  to  investigate  the  effect  of  these  forces  and 
moments  upon  the  motion  of  a  spinning  shell.  There  is  no  hope  of  ob¬ 
taining  a  closed-form  analytical  solution  for  this  motion,  so  it  is  appro¬ 
priate  to  solve  simplified  sets  of  equations  to  gain  insight  into  the  be¬ 
havior  of  projectile  motion.  The  starting  place  is  to  solve  the  equations 
of  motion  for  the  case  where  there  are  no  aerodynamic  forces  and 
moments  acting  upon  the  shell.  The  main  interest,  here,  is  in  the 
angular  motion,  since  the  translational  motion  is  simple.  The  next 
step  is  to  solve  the  equations  of  angular  motion  for  the  case  where  the 
only  torque  acting  is  an  overturning  moment.  These  two  simplified 
examples  aid  in  understanding  the  processional  and  nutational  motions 
of  a  projectile.  The  second  case  involving  the  overturning  moment  is 
mathematically  equivalent  to  the  motion  of  a  spinning  top  acted  upon  by 
gravity,  as  is  treated  in  many  text  books  on  Mechanics  (Ref.  51,  for 
example).  The  next  step  is  to  compare  the  motion  of  a  projectile  acted 
upon  by  an  overturning  moment  with  epicyclic  motion,  which  is  some¬ 
times  used  to  approximate  projectile  motion.  The  addition  of  damping 
to  epicyclic  motion  completes  the  picture,  more  or  less,  and  gives 
rise  to  an  understanding  of  dynamic  stability. 

With  this  background  in  mind,  it  is  possible  to  obtain  an  approxi¬ 
mate  solution  for  the  angular  motion  of  a  spinning  shell  by  use  of  the 
complex  notation  of  Section  II.  The  six-degree-of-freedom  equations 
of  motion  are  set  up  in  the  complex  notation  and  solved  under  simplify¬ 
ing  assumptions  as  to  the  importance  of  certain  terms  with  regard  to 
size.  Slowly  varying  terms  are  assumed  to  be  constant,  so  the  solu¬ 
tion  should  be  valid  only  along  a  short  section  of  the  trajectory.  The 
solution  thus  obtained  is  epicyclic. 


-K  -s-r- 
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This  section  is  closed  with  a  discussion  of  windage  jump  and  drift, 
both  of  which  result  from  the  angular  motion  of  the  projectile  and  forces 
normal  to  the  trajectory.  Windage  jump  must  be  considered  in  situa¬ 
tions  where  the  angle  betweenthe  gun  and  the  aircraft  velocity  vector  is  large. 
Drift  is  not  usually  considered  in  air  firings,  but  it  may  be  important 
for  new  rounds  under  development  at  long  ranges. 

The  material  in  this  section  is  entirely  tutorial.  The  equations 
derived  using  the  complex  notation  are  not  suitable  for  numerical  inte¬ 
gration,  and  equations  appropriate  for  that  purpose  are  derived  in 
Section  HI.  A  knowledge  of  this  material  is  essential  to  the  under¬ 
standing  of  projectile  motion,  windage  jump,  and  drift,  however. 

2.  The  Equations  of  Motion 

The  ?ix-degree-of-freedom  equations  of  motion  are  derived  from 
.  .->n's  laws,  which  are 


m 


du 
eft*  = 


F  +  mg 


where  m  is  the  projectile  mass,  g  is  the  acceleration  due  to  gravity, 
and  His  the  angular  momentum.  The  x^,x^,x^  coordinate  system  of 
the  projectile  is  chosen  in  an  incompletely  specified  manner  so  that 
the  longitudinal  axis  of  the  projectile  coincides  with  the  x^  axis,  but 
the  projectile  is  free  to  rotate  with  respect  to  the  x^,  x^  axes.  The 
angular  velocity  of  the  x^,x^,x^  coordinate  system  is  ft  and  it  follows 
that 


0) 


2 
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The  axial  and  transverse  moments  of  inertia  of  the  projectile  are 
A  and  B,  respectively,  and,  with  respect  to  the  x^x^.x^  coordinate 

system 

H  =  Xj  AtOj  +  x2B  +  x3Boj3  ^0) 

where  x  ,  x^,  and  x3  are  unit  vectors  in  the  indicated  directions.  The 
time  rate  of  change  of  a  vector  C  in  a  rotating  coordinate  system  with 
angular  velocity  ft  is  given  by 

t  =  C  +  ft  X  S  (61) 

where  a  dot  above  the  arrow  indicates  a  time  derivative  measured  m 
inertial  space  and  a  dot  below  the  arrow  represents  a  time  derivative 
observed  in  the  rotating  coordinate  system  (Ref.  51).  It  follows  that 


m(u1 

+  Uzu3  . 

'  u3Uz’  = 

+  mgL 

(62) 

m(u2 

+  aju,  ■ 

•niU3>  = 

F2  +  mg2 

(63) 

m(u3 

+  Vz 

•  Vi1  = 

F3  +  mg3 

(64) 

Awj 

=  G1 

(65) 

B-2 

+  WjAu^ 

-  f2^Bo>3 

=  °2 

(66) 

Bi3 

+  niBu2 

-  w2A“l 

=  G3 

(67) 

Use  has  been  made  of  =  u>z  and  ft3  =  u >y  The  set  of  angular  motion 
equations  will  be  solved  for  the  case  where  G  =  0,  and  for  the  case 
where  G  =  G£  =  0  and  G$  =  M  sin  0.  This  will  be  followed  by  a  dis¬ 
cussion  of  epicyclic  motion,  damped  epicyclic  motion,  and  dynamic 
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stability.  These  solutions  should  give  the  reader  sufficient  insight 
into  the  angular  motion  of  a  projectile  tc  understand  later  develop¬ 
ments. 

3.  Torque -Free  Motion 

It  is  of  interest  in  passing  to  ascertain  the  motion  of  a  projectile 
when  the  aerodynamic  torque  is  zero.  In  this  instance,  the  angular 
momentum  vector  H  is  constant  in  magnitude  and  direction,  and  it  will 
be  convenient  to  choose  a  coordinate  system  such  that  the  x  axis  lies 
along  H  (Fig.  5).  In  such  a  system,  the  motion  of  the  projectile  is 
particularly  simple. 

In  Fig.  5,  0  is  the  angle  between  the  projectile  x^  axis  and  x  axis, 

and  <j>  is  the  angle  between  the  y  axis  and  the  plane  containing  the  x^  and 

x  axes.  The  angular  velocity  of  the  projectile  is 
c* 

oj  =  x,  vb  +  l  *  +  x_  e 
1  T  x'  3 

where  iji  is  the  rate  of  change  of  the  orientation  angle  of  the  projectile 

measured  about  the  x,  axis,  and  1  is  a  unit  vector  in  the  x  direction. 

1  x 

But 


1  =  x,  cos  0  -  xn  sin  0 

x  1  2 

and  the  components  of  cj  in  the  x^.x^.x^  system  are 


Wj  =  41  +  cos  6 

(68) 

• 

=  -  4>  sin  0 

(69) 

• 

GO,  =  0 

(70) 

The  components  of  angular  momentum  are 
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H1  * 

Aw 

(71) 

H2  = 

Bw2 

(72) 

H3  = 

Bco3 

(73) 

Since  is  perpendicular  to  H,  it  foUows  that 

H3  =  Bco3 

=  B0  =  0 

(74) 

so 


=9=0 

and  0  is  constant.  Components  of 


H  =  XjAo,  + 

resolved  along  the  x  and  y'  directions  are 

AcOj  cos  0  -  Bu^  sin  0  =  H  (75) 

Auj  sin  0  +  Bu?  cos  0  =  0  (76) 

Since  0  is  constant,  and  co^  are  also,  because  this  set  of  equations 
can  be  solved  simultaneously  for  co^  and  co^  in  terms  of  0  and  H.  Since 
is  constant,  Eq.  (69)  shows  that 

4>  =  - 

is  constant,  and  it  follows  from  Eq. 
stant. 


u. 


sin  9 


(77) 


(68)  that  4j  is  also,  since  co^  is  con- 
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The  analysis  is  almost  complete  for  a  projectile  which  is  acted 
upon  by  no  aerodynamic  torques.  It  precesses  about  its  angular  mo¬ 
mentum  vector  with  constant  rate,  <{>.  An  expression  for  in  terms 
of  H  is  required,  however.  Since  =  0,  it  follows  from  Eq.  (67) 


niBw2  ~  CJ2‘^‘t*>l  =  a,2^*l  '  =  0 


and  since,  in  general,  f  0 


B  "l 


But  is  the  angular  velocity  of  the  x^.x^.x^  coordinate  system  about 


Xj  and 


=  4<  +  flj 


so,  from  Eq.  (68) 


=  <j>  cos  9 


FromEqs.  (75)  and  (76) 


Aoj^  =  H  cos  0 


FromEqs.  (78),  (79),  and  (80),  it  follows  that 


1  _  H 

4»  -  -g- 


It  also  follows  that 
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(82) 


W 


1 


_H 

A 


cos  9 


CO 


2 


H  . 
--rsin 


9 


(83) 


so  tOj  and  are  given  as  functions  of  9. 

4.  Motion  Under  the  Action  of  an  Overturning  Moment 

A  second  example  of  interest  is  the  motion  of  a  projectile  under 
the  influence  of  an  overturning  moment.  Figure  5  and  the  notation  of 
the  last  subsection  can  be  used,  but  H  is  no  longer  along  the  x  axis  as 
shown.  The  overturning  moment  is  always  perpendicular  to  the  plane 
containing  x^  and  x  and  is  taken  to  be 

=  M  sin  9  (84) 

where  M  is  a  positive  constant. 

This  problem  is  mathematically  equivalent  to  the  motion  of  a  spin¬ 
ning  top  under  the  influence  of  gravity,  in  which  case  M  =  mgf  ,  where  i 
is  the  distance  from  the  top  center  of  mass  to  the  pivot  (Ref.  51). 

By  Eq.  (65),  with  =  0 

=  Au>l  (85) 

is  constant.  Also,  since  G^  is  perpendicular  to  x,  the  component  of  H 
in  the  x  direction, 


H  =  Aw,  cos  9  -  Bw„  sin  9 
x  1  2 


(86) 


is  constant.  These  last  two  equations  yield 


W2 


H  cos  9  -  H 
1 _ _x_ 

B  sin  8 


(87) 
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and  from  Eq.  (69) 


H,  cos  0  -  H 

i  =■■  -  -t - , - 2-  (88) 

B  sin40 


FromEqs.  (68),  (85),  and  (88) 


H. 


+  =  -r 


cos  0 


H 


B  sin^0 


cos  0 


(89) 


4  « 

These  expressions  for  4>  and  4)  are  dependent  upon  0  alone.  An  expres¬ 
sion  for  0  in  terms  of  0  can  be  derived  from  consideration  of  the  energy. 
The  rate  of  doing  work  is  G  •  u,  and  from  Eq.  (61), 


fi  =  H  +  OXfi  =  G 


It  follows  that 


since 


H 


=  (H  +  SX  H) 


Ci) 


(a) 


u  •  (0  X  H) 


W1 

U2 

W3 

a 

^2 

w3 

Bw2 

Bw^ 

=  0 


But,  according  to  Eq.  (60) 
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H  •  co  = 


-  I  (A“l  +  B“2  +  Ba4) 


and  from  Eqs.  (70)  and  (84),  and  since  -  0 


G  •  w  =  G .a)-  =  0  M  sin  0  =  -  -t—  M  cos  0 


Then 


±  i  (A.2  +  Bu2  +  B„2  )  = 


-  dT  M  cos  0 


and  it  follows  that 


i(A“i 


1  +  B“2  +  Bw3 


)-■ 


M  cos  0  +  E 1 


where  E'  is  a  constant.  Since  is  constant,  this  expression  may  be 
written  as 


iB  (“ 


^  +  U3)  +  M  cos  0  =  E'  -  jAwj  =  E 


where  E  is  a  constant.  E',  of  course,  is  the  total  rotational  energy  of 
the  projectile,  anu  E  is  the  energy  associated  with  the  transverse 
motion.  Solution  of  this  expression  for  u!?  yields 


2  2E  2M  2 

“3  ~B - B~  C0S  9  -  w2 


and  from  Eqs.  (70)  and  (87), 
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cos  0 


(92) 


•2  2E  2M 
6  =  B~  '  B~ 


fH  1  C°S  9  -  HxV 

l  B  sin  0  I 


This  is  the  desired  expression  for  0  in  terms  of  0.  With  the  substitution 


w  =  cos  0 


(93) 


and  with  rearrangement,  it  can  be  put  in  the  form 


•  2 

w  = 

(1  -  w2)(a  -  pw)  -  (aw  -  b)2  =  F  (w) 

(94) 

where  a,  p,  a, 

and  b  are  constants  given  by 

a 

2E 
=  B 

(95) 

P 

2M 
"  B 

(96) 

a 

Hi 
"  B 

(97) 

b 

H 

X 

=  B 

(98) 

It  follows  that 

t  = 

w 

dw 

w  4  F(w) 

o 

(99) 

A  solution  for  t  may  be  found  in  terms  of  elliptic  functions  accord¬ 
ing  to  Ref.  51,  but  such  a  solution  would  serve  no  particular  purpose 
here.  It  is  sufficient  to  make  certain  observations  concerning  the  varia¬ 
tion  of  9  with  time.  F(w}  is  a  cubic  and  may  be  plotted  as  shown  in 
Fig.  6.  For  large  |  w  [,  the  dominant  term  in  F(w)  is  (3w  ,  and  since 
(3  is  positive,  F(w)  is  positive  for  large  positive  w  and  negative  for 
large  negative  w.  For  w  =  ±  1,  F(w)  equals  -(±a  -  b)  and  is  negative 
(unless  b  =  a  or  b  =  -a).  It  follows  that  the  plot  of  F(w)  vs  w  must  have 
the  characteristics  shown  in  Fig.  6,  with  two  roots  in  the  region 
-1  <  w  <  1,  and  a  third  root  in  the  region  w  >  1.  The  physical  motion 

of  the  projectile  can  occur  only  if  F(w)  is  positive,  and  this  occurs 

•  2  • 

between  w^  and  w^.  Otherwise,  0  would  be  negative  and  0  imaginary, 
or  else  cos  0  would  be  greater  than  zero,  which  is  impossible.  Thus 

Wj  =  cos  0^  <  cos  0  <  cos  0,,=  w^ 
or 


e2-6Sei 

The  angles  0^  and  0^  are  "turning  angles"  at  which 

0  =  ±sl F(cos  0) 

fc 

changes  sign.  When  0  reaches  0^,  0  changes  sign  and  approaches  0^; 

0  changes  sign  again  at  0^  and  0  approaches  0^.  This  motion  is  cyclic, 
and  typical  examples  are  sketched  in  Fig.  7  for  different  values  of 
parameters  a,  (3,  a,  and  b.  These  sketches  represent  0(t)  vs  4>(t)  in 
polar  coordinates  as  t  varies.  The  projectile  is  observed  to  move  with 
a  relatively  slow  precessional  motion  about  the  x  axis.  On  top  of  this 
precessional  motion,  a  faster  "nutational"  motion  is  observed  as  0 
moves  back  and  forth  between  0^  and  0^. 
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(c)  b  -  a  cos  0  changes  sign 
between  G  ^  and  0^ 

Figure  7 

Typical  Examples  of  Precessional  and  Nutational  Motion 
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The  motion  shown  in  the  three  examples  of  Fig.  7  is  explained 
as  follows.  Equation  (88)  may  be  written 


where  a  and  b  are  determined  by  initial  conditions.  In  Fig.  7(a),  is 

positive,  4>  continuously  increases,  and  the  curve  is  tangent  to  both 

the  inner  and  outer  circles.  In  this  case,  b  -  aw  >  0  for  all  w  between 

w  and  w  (or  for  9  between  0  and  0  ).  In  Fig.  7(c),  4>  changes  sign 

between  0  and  0  .  Clockwise  motion  is  exhibited  at  0  but  counter- 
1  ^  1 

clockwise  motion  occurs  at  0^.  In  this  case,  b  -  aw  =  0  has  a  root 
between  w,  and  w  ;  i.  e.  ,  w,<  b/a  <  w  .  In  Fig.  7(b),  b  -  aw  =  0  at  0  ,  in 

i  6.1  L.  1. 

which  case  both  4>  and  0  are  zero  at  the  same  time,  the  motion  stops 
instantaneously,  and  a  cusp  touches  the  inner  circle.  This  is  not  an 
exceptional  case,  as  one  might  think, since  the  values  of  a  and  b  are 

determined  by  initial  conditions.  At  the  instant  a  shell  is  fired  from  a 

•  • 

stationary  gun,  c)>  =  0  =  0  and  the  expected  motion  is  that  of  Fig.  7(b). 


A  type  of  motion  which  approximates  that  described  in  the  last 
subsection,  and  which  is  sometimes  used  to  approximate  the  angular 
motion  of  a  projectile,  is  epicyclic  motion.  If  a  wheel  is  attached  to 
the  rim  of  another  wheel,  as  shown  in  Fig.  8,  a  point  P  on  the  rim  of 
the  first  wheel  executes  epicylic  motion.  From  Fig.  8  it  is  seen  that 
the  coordinates  of  P  are 

0  =  A  cos  w  t  +  A  cos  w  t  =  I  0  I  cos  d> 

y  p  p  n  n  '  1 


9  =  A  sin  u  t  +  A  sin  t  =  8  sin  d> 

z  p  p  n  n  1  1 


In  terms  of  complex  numbers 


loj  t  iu)  t  . 

0=0  +  iG  =  A  e  P+Ae  "  =  I  0  I  e1<p  (100) 

y  Z  p  n  11 

Subscripts  p  and  n  stand  for  precession  and  nutation,  respectively. 

A  better  approximation  to  projectile  motion  is  provided  if  A^  and 
A^  are  damped,  that  is,  if 


A 


no 


e 


e 


-P 


t 

P 


t 

n 


so 


0.  = 


A  e 
po 


-(it+iojt 

P  P 


+  A 


no 


(101) 


In  this  case,  the  arms  A^  and  A n  continually  get  shorter  and  the  motion 
is  similar  to  that  sketched  in  Fig.  9.  This  motion  is  typical  of  a  dynam¬ 
ically  stable  projectile.  The  criteria  for  dynamic  stability  are 


Up  >  0  and  >  0  (102) 

Damping  results  from  the  inclusion  of  the  damping  moment 

Gh  =  -  Dwt  (103) 

where  D  is  assumed  to  be  constant  and 

+  X3U3  (104) 

The  differential  equations  for  an<^  “3  with  the  overturning  and  damping 
moments  included  are 
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l 


Bu>2  +  (Aojj  -  BQ^)  = 

-  Du>2 

(106) 

Bd>2  -  (Ay^  -  BOj)  to2  = 

-  Du ^  +  M  sin  0 

(106) 

Obviously,  the  damping  moment  opposes  and  tends  to  reduce  the  angular 
motion. 

6.  Complex  Notation 

The  equations  of  motion  can  be  written  in  terms  of  the  complex 
notation  of  Section  II,  and  an  approximate  solution  can  be  obtained 
(Refs.  1,  2,  3,  and  16,  for  example).  As  in  Section  II 


X  = 


U2  +  iu3 


P 


(w2  +  ico^)d 


Time  derivatives  are  transformed  into  derivatives  with  respect  to  arc 
length,  s,  measured  in  calibers,  where 


s 


(107) 


Differentiation  with  respect  to  s  will  be  denoted  by  a  prime,  e.g.  , 


u 


! 

1 


dul 

ds 


(108) 


Upon  utilization  of  Eqs.  (44)  through  (47),  and  Eqs  (62)  through  (67), 
the  equations  of  motion  in  complex  notation  become 


U1  +  i  P  X  -  p,  X 
u  2 


m 


KDA  + 


gld 


(109) 
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Aerodynamic  forces  and  moments  which  are  considered  negligible  have 
been  deleted. 

It  is  convenient  to  utilize  the  equation  of  translational  motion  which 
corresponds  to  the  direction  of  u.  It  is 


,2  2  „  . 

mu  =  -  pd  u  KD  +  mg^ 

where  g^  is  the  component  of  g  in  the  direction  of  u. 
stitution 


(113) 

With  the  sub- 


u 


•  u' 


Equation  (113)  becomes 

u_' 

u 


kd  + 


(114) 


(115) 
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•*  1-^jwnr 


This  equation  can  be  used  to  eliminate  u'/u  from  the  equations  for  X' 
andji1.  Upon  making  the  substitutions 


£ 


u 


cos  6 


(116) 


and 


(117) 


in  which  X  is  any  of  the  subscripts  D,  N,  T,  etc.  ,  Eqs.  (110)  and  (112) 
become 


(118) 


(119) 


Also,  the  equation  for  v'  becomes 


v 


(120) 


7.  Magnitude  of  Aerodynamic  Forces  and  Moments 

It  is  of  interest  to  estimate  the  magnitude  of  the  aerodynamic  forces 
and  moments  to  acquire  a  feel  for  their  effect  upon  projectile  motion. 
Accordingly,  the  aerodynamics  of  the  20-mm,  M56  round  will  be  used 
(Refs.  40jand  52  through  54).  For  the  20-mm  round, 


M  =  0.  22  lb 
A  -  0.  00013  lb  ft2 
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B  =  0.  00096  lb  ft2 
d  =  0.  06562  ft 

If  firings  are  made  at  sea  level, 

p  =  0.  076475  lb/ ft3 

and  parameters  needed  in  the  equations  of  Subsection  6  are 


Pd 

m 

=  0.  982  X  10 

md2 

A 

=  7.28 

md2 

B 

=  0.986 

A 

B 

=  0.135 

Aerodynamic  data  at  M  =  3  for  small  yaw  is 

Kd  «  0.  14 
—  1.2 
kf  -  0.1 

km  -  °-7 

Kh  &  2.0 
Kt  ^  -0.05 

ka  *  0.01 
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Except  for  Kp  and  K^,  these  values  were  taken  from  Ref.  53.  The 
value  for  Kp  was  calculated  from  the  curve  fits  in  Ref.  54.  No  data  is 
readily  available  for  K^,  so  a  more  or  less  typical  value  is  used.  The 
value  for  Kjj  is  the  value  obtained  when  the  ball  rotor  fuze  (Refs.  52 
and  53)  is  removed.  The  measured  value  with  the  ball  rotor  in  place, 
at  M  =  3,  is  Kpj  “  0.  05. 

The  value  of  projectile  velocity  corresponding  to  M  =  3  at  sea  level 
is 

u  =  3350  ft/sec 

A  value  for  v  at  the  gun  muzzle  can  be  obtained  from  a  knowledge 
of  the  projectile  muzzle  velocity,  V^,  and  the  twist  of  the  rifling,  n, 
where 

VM  =  3300  ft/sec 

n  =  25.586  calibers  per  turn 

(Ref.  54).  The  projectile  turns  through  one  revolution  when  it  travels 
a  distance  nd  in  the  barrel,  so 


1  _  2ir 

VM  '  ‘S'r 

and  the  initial  value  of  is 


Uj  =  12, 400  rad/sec 

The  initial  value  of  v  is 


v 


u 


0.24 


(121) 


The  terms  due  to  gravity  are  less  in  absolute  magnitude  than 
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2 

= 

1. 88  X  10" 

u 

while  values  of  or  vJ^  are 

JD 

1.4  X  10‘5 

JN 

1.2  X  10‘4 

vJF 

2.4  X  10'6 

JM 

7  X  10"5 

JH 

2  X  10'4 

vJT 

-1. 2  X  10‘6 

O 

1 

O' 

Using  these  typical  values,  it  is  seen  from  Eq,  (120)  that  v'  is 
small: 


rad2 

A 


JD  + 


'u 


0.  7  X  10'5 


The  solution  to 


v  *  —5 

—  =  -0.7X  10  3 
v 

should  be  valid  at  least  for  small  values  of  s  and  is 

-0.7X  I0"5s 

v  =  v  e 
o 

It  is  clear  that  v  does  not  change  much  over  a  trajectory. 

In  the  equation  for  X1,  Eq.  (118),  the  leading  terms  in  the  coef¬ 
ficient  of  X  are 
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.  nld 

JN  “  lv  JF  +  1  ~ 


JD  is  about  10%  of  and  the  gravity  term  is  completely  negligible; 
ivJp  is  retained  because  real  and  imaginary  terms  are  treated  sepa¬ 
rately.  The  leading  terms  in  the  coefficient  of  jj.  in  the  jx 1  equation, 
Eq.  (119),  are 


md* 

B 


.  .A  nid 

JH  -  ‘BV  +  ‘  — 


The  gravity  term  is  negligible,  and  is  small  compared  to  the  damp¬ 
ing  term  containing  J^.  The  set  of  equations 


(122) 


(123) 


with  v  constant  as  a  candidate  for  approximate  solution. 

8.  An  Approximate  Solution 

To  solve  Eqs.  (122)  and  (123)  with  v  constant,  a  complete  speci¬ 
fication  of  the  rotating  x^.x^.x^  coordinate  system  must  be  made  so 
thatflj  may  be  defined  (see  Subsection  2).  This  can  be  done  in  a  number 

of  ways,  but  for  present  purposes,  it  is  convenient  to  define  the  x.,x_,x_ 

i  “  ' 

system  so  =  0.  The  details  are  not  required,  however,  since  only  a 
partial  solution  of  projectile  motion  is  desired.  Approximate  equations 
for  X.  and  p,  will  be  derived,  whereas  the  development  of  equations  for 
the  direction  cosines  relating  the  x^x^.x^  coordinate  to  inertial  space 
will  be  omitted. 
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Since  the  influence  of  gravity  over  a  short  distance  s  along  a  tra¬ 
jectory  is  slight,  the  term 


( g2  +  ig3)d 

2 

u 

will  be  omitted. 

Also, 

i  =  cos  6 

is  approximately  equal  to  1  for  reasonably  small  yaw  angles,  e.g. , 
cos  6  =  0.978  at  8  =  12°,  and  so  we  set 

i  =  l 

The  rest  of  the  coefficients  are  slowly  varying  and  it  will  be  assumed 
that  they  are  constant.  The  solution  thus  obtained  will  not  be  exact, 
but  it  will  show  the  character  of  the  true  solution.  The  equations  to  be 
solved  may  be  written 


where 


+  axX  +  a2*x  =  0 

p '  +  b^X  +  b^p  =  0 


=  JN-ivJj 


a,  = 


-  i 


md^ 

"B 


(124) 

(125) 

(126) 

(127) 

(128) 

(129) 
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The  solution  is  straightforward  (see  Ref.  55).  Set 


\  =  X.  e 
o 


ks 


t*  =  >\,e 


ks 


(130) 

(131) 


where  \q,  (jlq,  and  k  are  constants.  Substitution  of  these  expressions 
into  the  differential  equations  yieldc 


(k  +  a  )X  +  a  a  =  0 
1  o  2ro 


bl\>  +  (k  +  b2*o  =  ° 


(132) 

(133) 


and  nonzero  solutions  exist  for  X.  and  n  only  if  the  determinant 

o  o 


k  +  aj  a2 


bl  k  +  b2 


(134) 


equals  zero  or  if 


k  +  ^al  +  b2*k  +  aib2  "  a2bl  =  ° 


(135) 


The  solution  of  this  quadratic  equation  is 


k  = 


-<al+b2>±M 


(al  +  V  -  4(aibZ  -  a2bl' 


(136) 
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Now 


and 


, ,  ,  T  md  T  .  A 

al  +  b2  =  JN  "  1  vJF  +  ~~B  JH  '  1  B“  V 


T  ,  md  T  .  A 
~  JN  +“B~  JH  "  1B’  V 


(a1  +  b/-4(a1b2-a2b1)  =  (a,  -  b/ 4  ia^ 


fT  md^  T  ,  .  A  "]  2  .  md^  I  T  .  _  ! 

“  [JN  "  B  JH+1BVj  "  4l-B  [vJT  +  lJMj 

/  md2  T  .A  .21 

“(4—  JM'{BV)  \ 


+  Ziv< 


A  A  md2  md2  j  1 

B"  JN  ”  g-  “B-  JH  ‘  £"~W~  Tf 


(137) 


(138) 


where  negligible  terms  have  been  deleted.  With 


T  ,  md  T 

3-  —  J  -NT  ^  “  J  r 


N  B  H 


b  =  tt-  v 


c  = 


4  md  J 
B 


(A  V 
\-BV) 


(• 


.  A  T  A  md2  T  .md2  T 

e  2v  <  B  JN  "  B  B  JH  2"B~  JT 


} 


(139) 

(140) 

(141) 

'142) 


we  have 


2k  =  -a  +  ib  ±  7c  +  ie 


(143) 
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According  to  DeMoivre’s  theorem  (Ref.  56) 


-  _  «/ 2 

'  c  +  ie  =  c  +  e 


2  I  a  ...  a 

[cos  2  +  1  slnT 


(144) 


where 


cos  a  = 


(2  2 
c  +  e 


sin  a  = 


2  2 
c  +  e 


(145) 


(146) 


.a  /  1  +  cos  a 
’  2  =  «/  2 


(147) 


a 

T  = 


1-ccs  a 
2 


(148) 


and  so 


'c  +  ie  = 


(n I  c2  +  e2  +  c)  +  i  Jj  (\fc^  +  e  -  c)  =  h  +  iq  (149) 


(It  has  been  assumed  that  a  is  in  the  first  quadrant;  if  not,  appropriate 
changes  must  be  made  in  the  signs  of  the  radicals  in  Eqs. (147)  and  (148).) 
It  follows  that 


2k  =  -a  +  ib  ±  (h  +  iq) 


(150) 


and  that 


-a  s+i8  s  -a  s  +  i8  s 

\  \  P  P  _i  i  n  rn 

X  =  X  e  +  X  e 

p  n 


(151) 
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where  X.  and  X  are  constants,  and 

P  n 

“p  =  i<a-h> 

(152) 

an  =  |  (a  +  h) 

(153) 

Pp  =  \  (»  +  q) 

(154) 

p„  -  i  (b  -  q) 

(155) 

A  similar  expression  exists  for  p,  and  it  is  seen  that  the  motions  of  X  and 
p  are  epicyclic  in  the  variable  s.  For  dynamic  stability,  it  is  required 

that 

a  >0  and  a  >  0 

P  n 

(156) 

It  is  observed  that  |  X  |  =  sin  6,  where  6  is  the  angle  between  the 
projectile  spin  axis  and  the  velocity  vector,  and  the  precessional  and 
nutational  motion  is  about  the  velocity  vector.  The  plane  of  yaw  is  the 
plane  containing  the  spin  axis  and  the  velocity  vector  ar„d  so  it  precesses 
with  the  projectile.  Since  the  lift  force  is  in  the  plane  of  yaw,  and  the 
Magnus  force  is  perpendicular  to  it,  these  forces  also  precess  about  the 
velocity  vector. 

The  analysis  of  this  section  is  approximate,  and  the  results  are 
tutorial.  A  better  approximation  is  given  in  Ref.  16  along  with  a  develop¬ 
ment  of  criteria  for  static  and  dynamic  stability.  A  treatment  similar  to 
that  of  Ref.  16,  with  slightly  less  general  results,  may  be  found  in  Ref.  1. 
In  passing,  it  is  observed  without  proof  that  the  condition  for  gyroscopic 
stability  is 

s*  >  1  (157) 
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where 


s*  =  (Av/B)2 

4(md2/B)JM 


A2  2 
A  u>^ 


4Bpd3u2K 


M 


(158) 


(Refs.  1  and  16). 

9.  Windage  Jump  and  Drift 

Windage  jump  and  drift  may  both  be  interpreted  as  deflections  of 
the  projectile  trajectory  away  from  its  initial  direction  of  motion  due  to 
aerodynamic  forces.  Windage  jump  is  caused  by  the  precessional  and 
nutational  motion  of  the  projectile  near  the  gun  muzzle,  while  drift  re¬ 
sults  from  the  effects  of  gravity  at  long  ranges.  Whereas  drift  has  not 
been  found  important  in  the  past,  it  may  be  significant  for  the  new,  heavy, 
high-muzzle -velocity  rounds  under  development. 

The  windage  jump  arises  from  the  side  forces  on  the  projectile 
and  the  precessional  motion  of  the  projectile  about  the  velocity  vector. 

The  side  forces  are  the  lift  force  in  the  plane  of  yaw,  and  the  Magnus 
force  perpendicular  to  the  plane  of  yaw.  These  forces  would  be  zero  if 
the  angle  of  attack,  6,  were  zero.  As  the  spinning  projectile  moves  down 
its  trajectory,  the  aerodynamic  moments  cause  it  to  precess  about  its 
direction  of  motion  like  a  top  under  gravity  and  so  the  plane  of  yaw  pre- 
cesses  with  the  projectile.  The  side  forces  are  carried  with  Ihe  plane  of 
yaw  and  the  changing  direction  of  the  forces  moves  the  projectile  first  in 
one  direction  and  then  in  another.  For  a  dynamically  stable  projectile, 
the  angle  of  attack  decreases  (on  the  average)  with  time  and  the  side 
forces  decrease.  The  center  of  mass  of  the  projectile  moves  downrange 
along  a  "spiral"  of  continuously  decreasing  radius.  The  angle  of  attack 
decreases  to  effectively  zero  usually  within  the  first  1000  ft  and  the  spiral¬ 
ing  motion  stops.  The  net  result  is  that  the  direction  of  motion  of  the 
projectile  is  changed  by  the  amount  of  the  windage  jump.  Equations  for 
the  windage  jump  are  given  in  Sections  IV  and  V.  An  explanation  of  the 
drift  is  as  follows. 
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5=rs 


At  long  ranges,  and  in  the  absence  of  gravity,  the  yawing  motion 
of  the  projectile  would  completely  die  out.  Gravity,  however,  causes 
curvature  in  the  trajectory  and  the  continually  changing  direction  of  the 
velocity  vector  induces  a  yaw  angle,  called  the  yav  ;f  repose,  between 
the  projectile  spin  axis  and  the  velocity  vector.  A  balance  of  aerody¬ 
namic  and  gyroscopic  moments  causes  the  projectile  to  move  with  its 
nose  pointed  slightly  up  and  to  the  right  of  the  trajectory,  and  the  effect 
is  to  deflect  the  trajectory  up  and  to  the  right.  A  treatment  of  drift  is 
given  in  Ref.  57. 


SECTION  IV 


EQUATIONS  OF  MOTION  FOR  COMPUTER  UTILIZATION 

1 .  General 

The  complex  number  representation  of  the  six-degree -cf-freedom 
equations  of  motion  derived  in  Section  III  is  useful  for  analytical  studies 
of  the  behavior  of  a  projectile  in  flight,  but  it  is  not  suitable  for  rapid 
computation  on  a  digital  computer.  Two  different  formulations  of  the 
equations  of  motion  are  presented  in  this  section.  The  first,  a  matrix 
formulation,  can  be  advantageously  programmed  on  a  digital  computer 
for  the  rapid  generation  of  trajectory  tables,  whereas  the  second  method, 
an  Euler  angle  approach,  is  more  adaptable  to  numerical  studies  and  to 
approximation.  Derivations  for  the  matrix  formulation  and  the  Euler 
angle  representation  follow  *u  Subsections  2  and  3.  In  Subsection  4, 
approximate  equations  are  dev  eloped  from  the  Euler  angle  set  for  the 
approximate  computation  of  trajectory  tables.  These  equations  are  use¬ 
ful  when  it  is  permissible  to  sacrifice  accuracy  in  favor  of  speed  in  com¬ 
putation.  Initial  conditions  are  not  give'  in  this  section.  For  computer 
studies,  these  parameters  may  te  chosen  arbitrarily,  whereas  they 
must  be  calculated  in  airborne  a ppli cations.  Calculation  of  initial  con¬ 
ditions  in  airborne  applications  is  treated  in  Section  VI. 

2.  A  Matrix  Formulation 

A  right-handed  X,  Y,Z  inertial  coordinate  system  is  defined  as 
follows:  X  and  Y  are  horizontal  and  Z  is  vertical.  The  X,Y,  and  Z 
coordinates  refer  to  the  center  of  mass  of  the  projectile  and  unit  vectors 
X,  Y,  and  Z  are  defined  in  the  indicated  directions. 

A  moving  x^.x^.x^  system  is  defined  in  Section  II  with  its  origin 
at  the  center  of  mass  of  the  projectile.  The  x^  axis  is  along  the  longi¬ 
tudinal  axis  of  the  projectile  and  is  positive  toward  the  nose.  The  x^ 
and  x^  axes  are  normal  to  x^  but  are  not  fixed  in  the  projectile.  The 
Xj.x^.x^  system  is  right-handed  and  x^,  x^,  and  x^  are  unit  vectors 
in  the  indicated  directions.  The  orientation  of  the  x^,x?,x^  system 
with  respect  to  the  X,  Y,Z  system  is  given  by  the  relation 
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The  angular  velocity  of  the  x^,x?,x^  system  is 

a  =  x1n1  +  x,j22  +  x3n3 


(159) 


while  the  linear  and  angular  velocities  of  the  projectile  are 

u  =  xu  +  X?U2  +  X3U3 

and 

"  =  *l“l  +  *ZW2  +  *3U3 
The  angular  momentum  of  the  shell  is 

H  =  x,Au  +  x.Bu,  +  x-Bw, 

1  1  C.  C.  J  J 

where  A  and  B  are,  respectively,  the  axial  and  transverse  moments  of 
inertia  of  the  projectile.  The  equations  of  motion  are 


F  =  m  u 


m(u  +  ft  X  u ) 


and 

G  =  £f  =  H  +  ft  X  H 
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A  dot  over  the  arrow  refers  to  a  time  derivative  measured  in  the  X,  Y,Z 
system,  whereas  a  dot  under  the  arrow  denotes  a  time  derivative  mea¬ 
sured  in  the  moving  system,  F  is  force,  G  is  torque,  and  m  is  the  pro¬ 
jectile  mass.  The  equations  of  motion  in  matrix  form  are 


The  restriction  that  the  axis  coincides  with  the  projectile  longitudinal 
axis  implies  that 

n2  =  «2  (162) 


(163) 

may  be  chosen  arbitrarily.  The  choice  is  inconvenient  for 

numerical  integration  because  a  very  small  step  size  (of  the  order  of 
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0. 1  X  2ir/w^)  will  be  required.  Also,  we  have  no  interest  in  the  angular 
orientation  of  the  projectile  about  the  x^  axis.  A  second  choice  is  to 
constrain  the  x^.x^.x^  coordinate  system  to  move  in  a  manner  such 
that  the  x^  axis  remains  in  the  horizontal  plane.  This  choice  is  con¬ 
venient  since  x  is  always  in  a  vertical  plane  and  £1  is  small.  A  third 
^  1 
convenient  choice  is  £1^  =  0. 

Forces  and  moments  are 


F.  ■ 

-pd2u2KDA  -  mgf  3 

(164) 

F2  = 

2  2 

-pd  uKNu2  -  pd  uvKpu3  -  mgm3 

(165) 

*1 

u> 

II 

-pd2  uKnu3  +  pd2  uvKpU2  -  mgn3 

(166) 

C1  * 

-pd3u2vKA 

(167) 

G2  ' 

3  3  4 

-pd  uvKTu2  +  pd  uKj^u3~pd 

(168) 

g3  = 

-pd^uvK^Uj  -  p  d^uKj^u^  -  p  d^uKjjt^ 

(169) 

where  g  is  ihe  acceleration  due  to  gravity.  The  aerodynamic  forces 
and  moments  are  the  same  as  those  of  Section  II,  Subsection  5  with 
negligible  terms  deleted.  If  aerodynamic  data  for  and  is  avail¬ 
able  rather  than  that  for  Kj-jA  and  K^,  the  following  expressions  may  be 
used 


DA  = 

kd 

cos  6 

-  Kt  sin26 

J-i 

(170) 

ii 

2 

W 

kl 

cos  6 

+  kd 

(171) 

Also 


cos  6 


v 


a 


(172) 

(173) 


u  = 


(174) 


The  aerodynamic  coefficients  (the  K's)  are  functions  of  Mach  number, 
sin  6,  v,  and  possibly  Reynolds  number  (see  Section  II). 


The  matrix  equation  for  the  velocity  is 


and  the  matrix  equation  for  the  time  rates  of  cnange  of  the  direction 
cosines  is 


(176) 


This  set  of  matrix  equations,  along  with  the  equations  for  the  aero¬ 
dynamic  forces  and  moments,  completely  describes  the  motion  of  a  pro¬ 
jectile.  For  numerical  solution,  it  lias  the  advantage  that  no  trigonometric 
functions  need  be  evaluated  from  angles  (cos  6  =  u^/u)  and  as  a  consequence 
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computations  are  fast.  It  also  has  the  advantage  that  no  division  by 
small  numbers  occurs,  such  as  division  by  sin  6  for  small  6.  The 
problem  of  division  by  small  numbers  often  occurs  in  formulations  in¬ 
volving  Euler's  angles.  But,  instead  of  the  expected  twelve  equations, 
there  are  eighteen!  The  redundancy  is  contained  in  the  equation  for  the 
direction  cosines.  One  might  use  the  properties  of  the  direction  cosine 
matrix  (orthogonality,  etc. )  to  reduce  the  number  of  equations,  but 
this  has  not  btcn  found  to  be  advantageous.  Instead,  all  eighteen  equa¬ 
tions  are  integrated  simultaneously.  (Actually  only  seventeen  need  be 
integrated  since  the  angle  of  rotation  of  the  projectile  about  its  spin 
axis  is  of  no  interest.)  Numerical  problems  occur  which  are  associated 
with  the  direction  cosine  matrix  not  remaining  orthogonal,  however. 

One  method  of  correcting  this  deficiency  is  to  replace  the  matrix  by 

j  (A  +  Aj  ) 

where  A  is  the  direction  cosine  matrix  and  A is  the  inverse  of  the 
transpose  (Ref.  58). 

3.  Euler  Angle  Development 

Large  portions  of  the  development  given  here  and  in  the  next  sub¬ 
section  is  taken  almost  verbatim  from  Ref.  59.  Changes  have  been 
made  as  necessary,  however,  to  clarify  and  adapt  Ref.  59  to  present 
purposes.  Portions  not  pertinent  to  present  needs  have  been  deleted. 

A  right-handed,  orthogonal,  rectilinear  £,t ],  £,  coordinate  system 
which  is  stationary  with  respect  to  the  ground  is  defined  with  £  mea¬ 
sured  down  range,  t]  vertically  up,  and  t,  to  the  right  as  seen  by  a 
person  facing  down  range.  A  second  right-handed,  orthogonal,  recti¬ 
linear  coordinate  system  1',  J’,  K'  is  superimposed  or.  £,q,£,.  The 
two  coordinate  systems  have  a  common  origin  O,  and  the  angular 
orientation  of  1',  J',  K'  with  respect  to  £,q,£  is  specified  by  the  two 
angles  a  and  0  (see  Fig.  10).  The  angles  a  and  0  are  by  definition  the 
azimuth  and  elevation  angles,  respectively,  of  the  bullet  velocity 
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*) 


Figure  lO 

Velocity  Vector  Coordinates 
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vector  u  as  measured  in  the  £,11,  £  coordinate  system,  and  1'  is 
always  along  u  .  The  azimuth  angle  a  is  measured  from  the  |  axis 
to  £',  the  projection  of  u  on  the  horizontal  plane,  and  the  elevation 
angle  9  is  measured  from  £'  to  u  or  to  1'.  The  angle  a  is  a  counter¬ 
clockwise  rotation  as  seen  from  above  (from  the  positive  i\  axis)  and 
0  is  a  counterclockwise  rotation  about  K*  as  seen  from  the  positive  K' 
axis.  Note  that  K'  and  £'  are  coplanar  with  E,  and  5  ar  d  are  horizontal, 
and  that  l'  and  J'  are  coplanar  with  £'  and  t|. 

Now,  define  a  1,  J,K  coordinate  system  parallel  to  1*,  J',  K'  but 
moving  witn  the  projectile  and  with  its  origin  at  the  center  of  mass  of 
the  projectile.  In  this  system,  define  a  right-handed,  orthogonal 
coordinate  system  1,  2,  3,  (see  Fig.  11).  The  1-2  plane  contains  the 
axis  of  symmetry.  A,  of  the  projectile  and  is  the  plane  of  yaw;  A  is 
directed  from  the  projectile  center  of  mass  toward  the  nose.  The  angle 
is  a  counterclockwise  rotation  about  1  as  seen  from  1  and  is  measured 
in  the  J -K  plane  from  J  to  2.  The  angle  6,  the  yaw  angle,  is  the  angle 
in  the  plane  of  yaw  between  the  projectile  velocity  vector  and  the  axis 
of  symmetry,  and  it  is  a  counterclockwise  rotation  about  3  as  seen 
from  3.  The  angle  4;  (not  shown  in  the  figure)  is  the  roll  angle  of  the 
projectile  measured  about  A,  and  it  is  assumed  to  be  a  counterclock¬ 
wise  rotation  as  seen  from  A;  is  the  roll  rate  of  the  projectile. 

It  is  observed  that  the  angles  6,  <}>,  and  ijj,  which  define  the  pro¬ 
jectile  orientation  with  respect  to  the  moving  1,  J,  K  system,  are 
Euler  angles.  Contrary  to  the  notation  of  Section  II,  the  1,  2,  3,  system 
is  attached  to  u  rather  than  to  the  projectile.  The  A,B,  3  coordinates 
used  here  correspond  to  the  x^,x^,x^  system  defined  in  Section  II. 

Figures  10  and  11  show  that  the  orientation  angles  of  the  1,2,3 
coordinate  system  are  a,  9,  and  <{),  and  that  1,2,  3,  constitutes  a  co¬ 
ordinate  system  with  one  axis  parallel  to  the  velocity  vector  u  and 
with  the  other  two  axes  rotating  at  essentially  the  precession  rate  of 
the  projectile.  The  angular  velocity  of  the  1,2,3  coordinate  system 
with  respect  tc  the  fixed  system  is 


87 


J 
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Figure  1 1 

Coordinates  of  Angular  Motion 
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where  T  ,  T.,,  and  T,  are  unit  vectors  along  tj,  K(or  K')  and  1,  respec- 

tj  K.  1 

tively,  and  a  dot  represents  differentiation  with  respect  to  time.  To  ob 

tain  components  of  ft  along  the  1,  2,  and  3  directions,  unit  vectors  1 

and  It,  are  resolved  along  the  1  and  2  directions  as  follows: 

K 


1R  =  1 £  sin<J>  +  1  ^  cos  6 


1  =  1 ,  sin  0  +  1  T  cos0 

T)  1  J 


=  l^sinP  +  (1^  cos<j>  -  l^sim},)  cos  G 

The  equations  for  Tj  and  are  obtained  directly  from  Fig.  li,  but 

1  is  obtained  from  Fie.  10.  Substitution  of  and  1  into  the  ex- 
ti  _  *1 

pression  for  ft  yields  the  1,  2,  3  components  of  ft  and  they  are 


ftj  =  <j>  +  a  sin  0  (177) 

ft^  =  a  costp  cos0  +  0  sin<j)  (178) 

ft^  =  -a  sincf>  cos0  +  0  cosij)  (179) 

Note  that  a  and  0  are  small  since  the  direction  of  the  velocity  vector 
does  not  cnange  much,  and  consequently  ft^  and  ft^  are  small,  ft^, 
therefore,  is  essentially  the  precession  rate  of  the  projectile. 

Figure  11  shows  that  the  angular  velocity  of  the  projectile  is 
w  =  ft  +  1a4j  +  1^6 


and  therefore,  the  components  of  angular  velocity  of  the  projectile 
resolved  along  the  A,  B,  3  axes  are 
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(180) 


to.  =  £2,  cosS  +  £1.,  sin6  +  tp  =  £2.  +  lit 
A  1  2  A 


(181) 

(182) 


u>  -  -Q .  sin6  +  £2  cos6  =  £2. 
13  1  L 

w3  =  +  6 


Components  of  the  angular  momentum  of  the  projectile  are 


ha 

=  a^a 

(183) 

=  Bcb 

(184) 

H3 

=  Bu)^ 

(185) 

where  A  and  B  are  th€  axial  and  transverse  moments  of  inertia  of  the 
projectile.  It  is  assumed  that  the  projectile  is  symmetrical  about  its 
longitudinal  axis,  and  that  the  moments  of  inertia  about  all  transverse 
axes  through  the  center  of  mass  are  equal. 

Newton's  second  law  for  rotational  motion  can  be  written  as 


5  =  H+fi'xfl  (186) 

where  G  is  the  aerodynamic  torque,  and 

£2'  =  u  - 

is  the  angular  velocity  of  the  A,B,3  coordinate  system;  that  is,  £2'  is 
equal  to  w  minus  the  axial  spin.  Components  of  £2'  along  A,  B,  and  3 
are 


flA  =  WA  '  *  =  flA 


£2 


B 


B 


fl3  =  u3 
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and  corresponding  components  of  Eq.  (186)  are 


GA  =  AiA 

(187) 

GB  =  B“b  +  <*>3A<l)A  '  nABa)3 

(188) 

g3  =  Bi3  +  oabWb  -  nBAWA 

(189) 

Equations  (187),  (188),  and  (189)  constitute  the  angular  motion  equa¬ 
tions  of  the  projectile,  and  GA,  Gg,  and  G3  are  the  aerodynamic 
moments. 

Newton's  second  law  for  linear  motion  with  velocity  and  force 
components  resolved  along  the  1,2,  3  directions  is 

F  =  m£u  +  12  X  u  j 

(190) 

and  also 

Fj  =  mu 

(191) 

F^  =  m£23u 

(192) 

F  3  =  -  m  12  ^  u 

(193) 

The  components  of  F  and  G  may  be  identified  with  the  proper 
aerodynamic  forces  and  moments,  from  Fig.  11  and  the  results  of 
Section  II,  with  careful  attention  to  geometry.  The  x  , x  ,x  directions 

1  u  J 

of  Section  II  may  be  identified  with  the  A,  B,  3  directions,  respectively, 
in  this  section,  and  in  Eqs.  (44),  (48),  and  (49),  F  becomes  F.,  F 

1  A  L 

becomes  Fg,  u^  becomes  Ug  =  -u  sin  6,  and  u^  is  zero.  If  gr  vity  is 
ignored  for  the  moment,  and  negligible  aerodynamic  terms  are  dis¬ 
carded,  it  follows  that 
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fa 

2  2 

=  - pd  u  KDa 

(194) 

fb 

7  2 

=  pd"u  sin  6 

(195) 

F3 

=  ~pd2u2  vKp  sin  6 

(196) 

But,  for  the  1  and  2  directions  of  this  subsection 


F,  -  cos  6  -  Fg  sin6  =  -pd2u2  (KgA  cos  6  +  KNsin2  6)  (197) 

F^  =  FA  sin6  +  Fg  cos  6  =  pd^u2(-KgA  +  Kj^  cos  6)  sin6  (198) 

From  Eqs.  (54)  and  (55),  and  when  the  gravity  terms  are  added  in,  it 
follows  that 


-pd2u2KD  -  mg  sine 

(199) 

pd2u2  Kt  sin6  -  mg  cos  0  cos  <j> 

JLi 

(200) 

2  2 

-p  d  u  vKp  sin  6  +  mg  cos  0  sin  <j> 

(201) 

In  Eqs.  (46),  (50),  and  (51)  for  the  aerodynamic  torques,  identification 
of  the  x.,x  , x_  system  with  the  A,B,  5  system  implies  that  G  becomes 
GA»  G^  becomes  Gg,  u^  becomes  Ug  =  -u  sin  6  (as  before),  u^  is  zero 
(as  before),  and«7  becomes  cjg.  When  negligible  terms  are  discarded 
it  follows  that 


ga 

=  -pd3u2vKA 

(202) 

gb 

3  2  4 

=  pd  u  vKj  sin6  -  pd  uKpjWg 

(203) 

G3 

=  pd3u2K^j  sin6  -  pd^uK^w^ 

(204) 
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.»* 


where 


v 


(205) 


Equations  for  a  and  0  can  be  obtained  in  terms  of  £2^  and  S2^  from 
Eqs.  (178)  and  (179). 


a  cos  0  =  cos  4>  -  f2,  smiji 


0  =  +  ^3  cos  4* 


From  Eqs.  (192)  and  (193)  it  follows  that 


a  cos  0  =  — - —  -  F,  cos  d>  -  F_  sin  6 
mu  1  3  Y  2  Y 


•  1  ; 

0  =  -  !  -  F,  sind>  +  F.  cos  d> 

mu  ^3  T  2  Y 


and  finally,  from  Eqs.  (200)  and  (201) 


.  pd 

a  =  — 


m  cos 


- -  r 

os  0  I 


vK„  cos  4>  -  Kt  sintj? 
x  L 


sin  6 


q  =  2.  —  _  r vK„  sin  A  +  Kt  cos  <j>  j  sin  6  -  —  cos  0 
m^FYL  Y  J  u 


(206) 

(207) 


(In  numerical  integration,  division  by  cos  0  causes  numerical  trouble 
when  0  is  near  ±  90°. )  From  Eqs.  (191)  and  (199) 


u  =  - 


,2  2 

p  d  u 

m 


Kn  -  g  sin  0 


(208) 


Equations  (206),  (207),  and  (208)  are  the  equations  for  u  in  polar  co- 
ordirates.  The  equations  for  projectile  position  can  be  obtained  from 
examination  of  Fig.  10.  They  are 
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u  cos  0  cos  a 

(209) 

u  sin0 

(210) 

-u  cos  0  cos  a 

(211) 

By  use  of  Eqs.  (187),  (188),  (189),  (202),  (203),  and  (204),  the 
angular  motion  equations  may  now  be  written  as 


a<!)a  =  -pd  uo>aKa 
b«Lb  +  Acja  -  flABa>2 


(212) 


4  4 

=  pd  uo^Krj-,  sin6  -  pd  uKHu>. 


HB 


(213) 


B<L3  +  nABuB  - 


3  2  4 

=  pd  u  K^sinS  -  pd  uK^w^ 


(214) 


where  v  has  been  replaced  by  the  right  side  of  Eq.  (205).  £2A  is  obtained 
as  follows:  use  Eqs.  (178)  and  (179)  to  calculate  £2^  and  then  £2^ 
may  be  obtained  from  Eq.  (181),  i.  e. 


£2 


1 


£2^  cos  6  -  w, 
sin  6  " 


(215) 


(Division  by  sin  6  causes  numerical  integration  problems  when  6  is 
small.)  Given  £2 ^  and  £2^,  Eq.  (180)  yields  £2A>  Also,  4>  and  6  may  be 
obtaint  d  from  Eqs.  (177)  and  (182). 

3.  1  Summary  of  Equations  -  It  is  convenient  to  collect  together 
the  complete  set  of  six-degree-of-freedom  equations.  They  are  as 
follows: 
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a  =  — P.4...1*  Ty  K_  coo  4»  -  K  sin  <J»1  sin  6 

m  cos  0  r  T  Li  J 

*  ■  #  [»  “f  sin  <j>  +  cos  4>J  sin  6  -  £ 


n. 


si. 


SI. 


,2  2 

-  £-1-2-  Kn  -  g  sin  0 
m  T)  6 

,  * 
a  cos  0  cos  $  +  0  sin  $ 

-a  cos  0  sin  <{>4-0  cos  <{> 

-<*>_,  +  SI ,  cos  6 
o _ i 

sin  5 


nA  =  cos  6  +  fl.,  sin  8 
A  1  2 


w. 


pd  U  XT 

Tr“AKA 


;  95 

j 

X 


cos  0 


-H  -  -  d  sin  8 

§  -  a  cos  0  cos  a 

rj  =  u  sin  0 
£  =  -u  cos  0  sin  a 

Initial  conditions  are  treated  in  Section  VI. 

4.  Approximate  Equations  for  Large-Yaw  Computations 

This  subsection  contains  a  derivation  of  the  approximate  equa- 
n»  '•f  .  ’otion  used  at  Eglin  Air  ‘,''>rce  Base  for  trajecto>*y  table 
ale il* tioi-  (c  the  M5i  round  An  equivalent  set  has  also  been 

•  •see  in  .he  ...imputa  ion  oi  ; u... le  j  b,;  LA-' •  (Re!.  3i).  Presumably,  the 
n».  tn  >d  c  n  l  adant:,-1  r  -e«v  rounds  unde  .,  uevdopment.  Another  der¬ 
ivation  of  r:  °  two  approximate  equations  of  angular  motion  can  be  found 
in  Ref.  6C.  This  derivation  is  quite  "odious,  however,  and  derivations 
of  the  other  equab-ons  are  lacking.  The  development  herein  is  that  of 
del.  39  \,Rh  minor  changes  lor  cinity  and  is  based  upon  Ine  equations 
r f  tt-.e  iast  sub'-ection.  f=*e  equations  of  li-.ear  motion  must  os  formu- 
atc-d  it,  tii--  irq,f,  system,  however. 

The  ••sason  for  using  these  approx.mate  equations  is  thev  _au  be 
evaluated  much  ster  on  a  digital  computer  than  can  the  six-degree-of- 
freedom  equations,  Tne  angular  motion  is  approximated  and  the  nata¬ 
tion  motion  (Section  TTt)  is  eliminated.  Since  the  fine  detail  of  the  nuta- 
tional  met'o?.  is  absent,  and  only  the  ^recessional  mot  o*"  is  left,  the 
numerical  step  s  can  be  increased  and  thereby  computer  time  is 
reduced.  The  derivation  proceeds  as  follows. 

Equation  (192)  is  combined  hero  withEq.  (200)  for  future  ref¬ 
erence: 
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(216) 


2  2 

mSl^u  =  pd  u  sin6  -  mg  cos0  cos4> 

Equations  (199),  (200),  and  (201)  transformed  to  the  1,  J,  K  system 
are 


F1 

2  2 

=  -pd  u  Kp  -  mg  rrin  0 

(217) 

FJ 

=  .  pd2u2  |k^  cos<j>  +  vKp  sin<j>Jsin6  -  mg  cos0 

.  (218) 

fk 

=  pd^u^^K^  sin<j>  -  vKp  coS(J>Jsin6 

(219) 

And  in  terms  of  the  £,q,r,  system 

F? 

=  F^  cos0  cosa  -  Fj  sin0  cosa  +  F^.  sina 

(220) 

F 

n 

=  Fj  sin0  +  Fj  cos0 

(221) 

=  -  F  cosO  sina  +  F_  sin  0  sina  +  F„  cos  a 

1  J  is. 

(222) 

IfEqs.  (217),  (218),  and  (219)  are  combined  with  Eqs.  (220),  (221), 
and  (222),  and  use  is  made  of  Eqs.  (209),  (210),  and  (211),  the 
following  expressions  are  obtained 


f£  = 

•s 

2  •  2  2  r  1 

-pd  uKp£  -  pd  u  cos<})  +  vKp  sin<})  sin6sin0cosa 

>  ->  r 

+pd“u~  sin<|)  -  vKp  cost}) 

j  sin6  sina 

(223) 

F  = 

n 

-pd^uKj^r)  +  pd^u^  cos<|>  +  vKp  sin<fij  sin6  cos0  -  mg 

(224) 

F£  = 

2  •  2  2  r  1 

-pd  uK^  +  pd  u  cos4>  +  vKp.  sin»j)  sin6  sin0  sina 

+pd^u^  sin<£  -  vKp  cos4> 

sin6  cosa 

(225) 
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F  inally , 


me  -  f4 

(226) 

m‘rj  =  F 

n 

(227) 

me  =  f5 

(228) 

4. 1  The  Siacci-Type  Approximation  to  the  Force  Equations  -  As 
in  the  derivation  of  the  Siacci  equations  (Section  V),  l3  is  defined  as  a 
vector  along  the  initial  velocity  vector  of  the  projectile  to  a  point 
(almost)  directly  above  the  shell,  and  Q  is  defined  as  a  vector  pointing 
vertically  down  from  the  tip  of  P  (almost)  to  the  projectile.  The  dis¬ 
tance  from  the  tip  of  Q  to  the  projectile  is  the  swerve  Sf .  The  swerve 
is  defined  to  be  the  displacement  of  the  shell  from  1?  +  Q;  i.e. ,  it  is 
the  displacement  of  the  shell  from  its  straight-line  path  due  to  forces 
other  than  gravity  (or  due  to  aerodynamic  forces  but  not  to  gravity). 

The  Siacci  equations  ignore  the  swerve,  but  it  is  accounted  for  later 
by  the  windage  jump.  The  procedure  here  is  similar;  differential  equa¬ 
tions  for  the  swerve  will  be  obtained  from  the  differential  equations  of 
motion  by  subtracting  out  the  Siacci  equations  for  P  and  Q. 


By  definition,  the  initial  velocity  vector  uq  is  taken  in  the  £,tj 

plane  so  that  a  =  aQ  =  0.  Since  P  lies  along  uq,  the  elevation  angle  of 

P  above  the  horizontal  (and  in  the  plane)  is  6^  (see  Fig.  12). 

Siacci-type  equations  for  P  and  Q,  and  for  swerve  components  S.,  S  , 

5  *n 

and  S  are  derived  as  follows.  In  terms  of  P,  Q,  St,  S  ,  and  S  ,  the 


coordinates  £,  tj,  and  £,  are 


£’ 


V 


4  are 

£  = 

P  cos  8  +  S 
°  £ 

(229) 

T)  = 

P  sin  8  -  Q  +  S 

(230) 

O  T) 

£  = 

st 

(231) 
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Thus,  the  force  Eqs.  (223)  through  (228)  become 


m(P  cos0  +  S’  =  -pd^uKr.(P  cos9  +  S..) 
'  O  5  U  O  § 


-pd2u2  cos<j>  +  vKp  sin<f>j  sin6  sin0  cos  a 

•f  pd^u2  sin<f>  -  vKp  cos<|>j  sin6  sin  a 

(232) 

m(P  sin0  -  Q  +  S  1  =  -pd2uKr,(P  sin0  -QrS  ) 
o  rj  r  D  o  q 

+  pd2u2  |k^  cos(j>+  vKpSin4>j  sin6  cos0  -  mg 

(233) 

mS^  =  -pd2uK^S^ 

+  pd2u2  |kl  cos<(>  + vKp  sin<j>j  sin6  sin0  sina 

+  pd2u2  sin<|>  -  vKp  coS(|>j  sin6  cos  a 

(234) 

By  definition,  the  Siacci  equations  for  P  and  Q  are 

mP  =  -pd2uKDP 

(235) 

mQ  =  -pd2uK^Q  +  mg 

(236) 

with 

u  =  (P2  -  2PQ  sin  eQ  +  qM  (237) 


They  are  obtained  from  the  equations  above  with  the  swerve  set  to  zero 
and  with  all  aerodynamic  coefficients  except  set  to  zero.  Upon 

substitution  of  Eqs.  (2  )5)  and  (236)  into  Eqs.  (232)  and  (233),  terms 

•  • 

containing  P  and  Q  cancel  and  the  equations  for  the  swerve  components 
S.  and  S  remain.  The  following  are  obtained: 
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mS . 


COSQ 


mS  = 


-pd2uKDSe 

-pd2u2  ^Kp  cosiji  +  vKp  sincfjj  sin6  sin0 

+pd2u2  j^Kp  sin<}>  -  vKp  cos^j  sin5  sina 

2  *  2  2  f  1 

-pd  uKpS^  +  pd  u  I  Kp  cosd>  4  vKp  sincjjJ  sin6  cosG 


(238) 

(239) 


The  equation  for  is  Eq.  (234). 

•  • 

The  equations  derived  up  to  this  point  would  be  exact  if  S^,  S^, 
and  were  included  in  Eq.  (237).  Little  error  results,  however,  if 
these  terms  are  ignored.  For  Eqs.  (238)  and  (239)  to  be  used  as 
written,  it  would  be  necessary  to  integrate  the  equations  for  a  and  0. 
But  since  0  and  a  do  not  change  much  from  their  original  values  of 
0  =  0Q  and  a  =  0°  for  those  trajectories  which  are  of  interest,  0  will  be 
replaced  by  9o,  cos  a  will  be  replaced  by  1  and  sin  a  by  0.  The  dif¬ 
ferential  equations  for  the  swerve  become 

mS^  =  -pd  uKj-jS^  -  pd  u  I  Kp  coS(j>  +  vKp  sincjjJ  sinb  sin0Q  (240) 

mS^  =  -pd2uKj-jS^  +  pd2u2  costf)  +  vKp  sin(j)j  sin6  cos9q  (241) 

••  2  *  2  2  r  1 

mS^  =  -pd  uKpS^  4  pd  u  ^KpSin<J>  -  vKp  cos^J  sin6 
Equations  (240)  and  (241)  can  be  combined.  Take 


(242) 


S.  =  -S,.  sin  0  +S  cos  0 
I  |  o  f|  c 


S  =  S.  cos  0  +S  sinO 
II  §  o  q  o 


'243) 

(244) 


where  andS^  are  in  the  £,q  plane  and  are  respectively  perpendi¬ 
cular  and  parallel  to  S  is  positive  above  P,  and  it  will  be  shown 

••  •• 

that  is  identically  zero.  The  equations  for  and  as  derived 
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'•npl  TTW-Ji  Tgf**t 


from  Eqs.  (240)  and  (241)  are 

••  2  •  2  2  r  1 

mS^  =  -pd  iiK^S^  +  pd  u  cos$  + vKp  sin<j>  sin6  (245) 


mS..  =  -pd^uKnS, 


(246) 


m  w 

Since  at  t  =  0,  S^  =  S^  =  0,  it  is  evident  that  is  identically  zero  for 

all  t.  Values  of  S,.  and  S  ,  then,  are 
S  9 


St  -  -S  sin  0 
i  L  o 


(247) 


S  =  S  cos  0 
r\  i  o 


(248) 


and  it  is  observed  that  the  three  differential  equations  for  s-jverve  have 
been  reduced  to  two,  Eqs.  (242)  and  (245). 

If  Eqs.  (242)  and  (245)  are  used  in  their  present  form  to  com¬ 
pute  firing  tables,  the  initial  precession  angle  <j>o  will  have  to  be  one 
of  the  table  entry  parameters;  if  one  more  coordinate  transformation 
is  made,  however,  this  can  be  avoided  and  the  volume  of  the  tables 
will  he  greatly  reduced.  Components  of  swerve  in  and  perpendicular 
to  the  plane  of  initial  yaw  are  used.  See  Fig.  13  in  which  Pis  direct¬ 
ed  into  the  page.  The  new  swerve  components  are 


S2  =  S1  COS*o  +  S£  sin*o 


(249) 


S^  =  -S^  sin4>o  +  S^  cos<|>c 


(250) 


with  S^  in  the  plane  of  yaw  and  normal  to  it.  The  differential  equa¬ 


tions  for  S^  and  S^  are 


mS^  =  -pd^uK^S^  +  pd^u^  ^K^cos^1  +  vKp  sin<)>'j  sin6  (251) 

mS^  =  -pd^uKjjS^  +  pd^u^  ^K^sin<j>'  -  vKpCOS<j>'j  sin6  (252) 
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where 


4>‘  -  4>  -  4>0 

and  <j>'  =  0  at  t  =  0.  Given  4>q,  0q,  and  values  of  and  S^, 

S.  ,  S  ,  and  S  can  be  found  from 
g  n  b 

values  of 

=•  -(S^  cos<j>o  -  S3  sin4>o)  sin0Q 

(253) 

=  (S^  r.os4>o  -  S3  sin4>Q)  cos0q 

(254) 

S£,  =  S2  sin4,o  +S3cos<i>o 

(255) 

If  angular  windage  jumps  (in  milliradians)  are  desired,  they  are 


J2  =  lo3:VP  (256) 

J3  =  103S3/P  (257) 

In  Ref.  59,  S  and  S  were  denoted  by  the  symbols  S  and  S  , 

c  j  x  y 

respectively.  The  notation  used  here  is  preferable,  since  and  S3 
are  in  the  2  and  3  directions,  respectively,  of  Fig.  11  when  viewed 
at  the  time  of  firing.  Calculations  show  that  and  J3  approach  more 
or  less  constant  values  as  the  projectile  moves  down  its  trajectory. 
Accordingly,  windage  jump  can  be  accounted  for  by  calculating  the 
point-mass  trajectory  defined  by  Eqs.  (235),  (236),  and  (237),  and  by 
moving  the  resultc.nt  hit-point  vector  I?  +  Q  through  angles  and  J3> 
in  and  perpendicular,  respectively,  to  the  initial  plane  of  yaw,  as  seen 
from  the  point  of  fire.  The  windage  jump  correction  amounts  to  a  slight 
change  in  the  direction  of  P  or  u. 

The  set  of  equations  derived  thus  far  is  complete  and  can  be  used 
to  compute  firing  tables.  If  the  spin  is  sufficiently  large,  however, 
the  equations  of  angular  motion  can  be  simplified  and  the  computation 
time  can  be  reduced  considerably.  The  following  subsection  contains 
the  simplified  equations  and  the  justification  for  simplification. 
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4.2  The  Approximate  Equations  of  Angular  Motion  -  The  crite¬ 
rion  for  simplifying  the  equations  of  angular  motion  is  that  the  spin 
must  be  large  enough  for  the  motion  of  the  shell  to  be  essentially  gyro¬ 
scopic.  In  this  case,  the  angular  momentum  vector  H  will  be  nearly 

...  .  M  *— » 

parallel  to  1^,  and  the  component  of  H  perpendicular  to  1A  will  be 
small. 

If  s  is  a  unit  vector  in  the  direction  of  H, 

H  =  Hs  (258) 

and  Newton's  second  law  for  angular  motion  becomes 

G  =  il  =  Hs  +  Hs  (259) 


and  so 


oXG  =  Hs  Xs  =  Hujj 


(260) 


where 


w 


H 


=  s  X  s 


(261) 


is  the  angular  velocity  of  H,  as  can  be  seed  from  Fig.  14.  Under  the 
assumption  that  Aw A  »  j«b  (see  Fig.  15),  it  follows  that  H  is  nearly 
parallel  to  1^,  or  that 


S  1 


(262) 


and  that  the  angular  velocity  of  H,  that  is  u^,  is  approximately  equal 
to  the  transverse  angular  velocity  of  the  projectile, 


u> 


H 


—  u) 


T 


=  ‘b"b 


+  l3u3 


(263) 
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t  =  lim  s  (t  +  At)  -  s  (t)  _  lim  -AS 
At  —  0  At  “  At  —  0  T  At 

=  s'  X  s'  -  s'  X  T  p  -  p 

s  is  a  unit  vector  along  H 

t  is  a  unit  vector  along  As 

v  is  a  unit  vector  pointing  out  of  the  page 


Figure  14 
„  «. 

Proof  that  s  X  s 


Then  H  s*AtoA  and  Eq.  (260)  becomes 

AwA  (TBuB  +T3w3>  "TAX  ^ 


(264) 


It  follows  from  Eqs.  (?03)  and  (204)  that 


4  .  4 

u^^Abij^  a*  pd  uw^K,p  sin 6  -  pd  uKyiOg 


(265) 


3  2  4 

-ujgAto^  =  pd  u  Kj^sin6  -  pd  uK^w^  (266) 

This  approximation  is  equivalent  to  deleting  terms  containing  B  in  Eqs. 
(213)  and  (214).  The  advantage  is  the  elimination  of  terms  containing 
w  B  and  w^.  This  negates  the  need  for  integrating  to  obtain  and 
and  also  smoothes  out  the  nutational  motion  as  is  explained  by  Reed 
(Ref.  60).  The  smoothing  of  the  nutational  motion  alleviates  the  require¬ 
ment  for  a  small  step  size  (At)  in  numerical  calculations. 

•  * 

The  approximate  equations  for  <{>  and  6  are  obtained  by  solving 
Eqs.  (265)  and  (266)  simultaneously  for  u>g  and  oo^. 


OJ 


B 


U 


3 


-AuApd3u2KMsin5  +  p2d8u2uAKHKT 
(Awa)2  +  (pd4uKH)2 

4  ?  2  7  3 

Apd  uwAK,pSin6  +  p  d  u  sin6 

(A(oa)2+  (Pd4uKH)2 


sin6 


(267) 


(268) 


4  2 

For  high  spin,  (pd  aKJ  is  usually  much  less  than  one  percent  of 

2  4  2 

(A^a)  ;  hence,  (pd  uK^j)  is  dropped  from  the  denominator  and  these 
equations  become 
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Tfm**&V*  r* 


—~*&fT_  <2»^‘rT>7'7J  ST  ^-sp?.?- ’ 


“B  *  •^[KM‘1TKHKI 

"3  -  ^  [KT  *  f£  Vm] 


sin  6 


sin  6 


(269) 

(270) 


» 

The  equation  for  6  can  now  be  derived  from  Eqs.  *182),  (2i6),  and  (270). 
It  becomes 


6  =  - 


Ik  .s5!/kt  + 

m  L,  A  \  T 


Pd~ 

a  2 

Av 


khkm 


}]si 


sin6  +•£  cos0  cos<|>  (271) 


The  gravity  term  must  be  dropped  from  this  equation,  since,  for  very 
small  yaw,  the  aerodynamic  term  will  be  small  compared  to  the  gravity 
term  and  erroneous  results  will  be  obtained  in  numerical  integration 
attempts.  The  term  containing  is  probably  negligible.  Calcula¬ 

tions  show  this  to  be  so  for  the  20-mm,  M56  round,  but  a  judgment  will 
be  needed  for  new  rounds  under  development.  To  obtain  the  time  de¬ 
rivative  of  cos  6,  Eq.  (271)  is  multiplied  through  by  -  sin  5.  Then 


d  COS6  =  [kl  - ^  (kt  +  *£  KhKm|J  si„26  (272) 

The  equation  for  4>  is  derived  from  Eqs.  (177),  (178),  (181),  and 

(269)  under  the  assumption  that  a  and  0  are  negligible  in  comparison 

•  •  *  • 

with  <{>.  If  a  and  0  are  neglected,  the  equation  for  <j>  is 


.3  2 

•  =  pdj^  K 

T  Acj  .  M 


(273) 


The  term  containing  was  dropped  because  it  is  negligible. 
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4.3  Summary  of  Equations  -  The  approximate  equations  are 
summarized  as  follows  for  the  convenience  of  the  reader: 

mP  =  -pd^uKj-jP 


mQ  =  -pd  uKnQ  +  mg 

•• 

mS., 

U 

mS„ 


D 

-pd^uKpS^  +  pd2u2  |"kl  cos<j>'+  vKpSin^'j  sin  5 


d  cos  6 

— at- 


-  -pd  uKdS3  + 

2 

pd  u 

K  - 

m 

L 

,3  2 

pd  u 

-  K. 

Aw  A 

M 

u2  |* 
!»*[* 


sii 


^  Ikt + ^4  vm}1  si"2  6 

l  Av  j 


“  ■  ft 


P2  -  ZPQ  sin  9  +  Q2 
o 


V  - 


WAd 


4>' ~  <!>  -  4>, 


£  =  P  cos  9  +  Se 
*  o  § 


n  =  P  sin  0  +  S  -  Q 
1  O  q 


6  =  Sr 


S.  =  -  (S,  ;.os  d>  -  S,  sin  A  )  sin  6 
£  2  To  3  To  o 


Cj  =  (S„  cos  A  -  S.  sin  A  )  cos  0 

q  '2  ~o  3  To  o 


S,  =  S_,  sin  A  +  S_  cos  <j> 
t,  2  o  3  o 
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SECTION  V 


THE  SIACCI  METHOD 


1 .  General 

The  Siacci  method  (Refs.  I,  24,  27,  and  61)  is  a  means  of  obtaining 
approximate  solutions  to  projectile  trajectories  by  use  of  tables.  It  is 
applicable  in  many  situations  where  high  computational  accuracy  is  not 
essential,  such  as  in  preliminary  design  and  in  new  concept  studies,  and 
in  particular,  i*  is  valid  for  relatively  short  trajectories  where  gravity 
drop  is  not  appreciable  and  yaw  is  small.  The  calculational  ease  with 
which  the  method  can  be  utilized,  e.  g.  ,  in  hand  computations,  or  in  an 
airborne  fire-control  computer,  makes  the  Siacci  method  invaluable  in 
many  modern  applications  despite  its  early  origin.  The  original  method 
was  devised  sometime  around  1880  by  F.  Siacci  of  Italy.  The  treatment 
given  here  is  that  of  Ref.  61  with  minor  changes. 

With  the  development  of  the  20-mm,  M56  round  in  the  late  1950's, 
the  trend  in  ballistics  calculations  was  away  from  use  of  the  Siacci 
method  and  toward  more  sophisticated  calculations  since  modern  com¬ 
puting  equipment  was  becoming  available.  Sophisticated  methods  are 
now  well  suited  to  ground-based  ballistics  investigations,  but  for  airborne 
fire-control  calculations,  onboard  computers  are  still  somewhat  limited 
in  capacity.  The  Siacci  method  is  still  a  candidate  for  airborne  fire- 
control  calculations. 

A  description  of  the  basic  Siacci  method  follows  and  a  derivation 
is  given  in  Subsection  2.  In  its  basic  form,  the  trajectory  as  given  by 
the  Siacci  solution  is  as  follows: 


t 


C 

aao 


T(u/a  )  -  T(u  /a 
o  o 


) 

o 


(274) 
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[s(u/ao)  -S(uQ/ao)J 


(27  6} 


Q  =  (^7) 2  [A(u/ao} "  A(uo/ao)  -  W  i- p]  (276) 


in  which  S,  T,  I,  and  A  are  trbulated  functions  of  U  =  u/a  :  t  is  the 

'  o 

time  of  flight,  P  is  the  "pseudorange1'  along  the  initial  velocity  vector 

imA 

u  »  and  Q  is  the  gravity  drop  (Fig.  16).  The  parameter  <r  is  the  relative 
air  density  at  the  firing  altitude,  a  is  the  ratio  of  the  speed  of  sound  at 
the  firing  altitude  to  that  at  sea  level,  and  C  is  the  ballistic  coefficient 
given  in  terms  of  the  projectile  mass  m,  in  pounds,  and  diameter  d,  in 
feet,  as  shown  below; 


C  =  — — -y  (lb/in3  ) 


144d 


The  parameter  u  is  defined  as  follows 


and  is  the  independent  variable. 

The  tabulated  functions  S,  T,  I,  and  A  are  calculated  by  numerical 
integration  of  the  following  equations; 


1 

G(U) 


(277) 


dT 

~3U~: 


UG'(U) 


(278) 


dl 

~3U~ 


U2G(U) 


(279) 
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dA. 

W  ' 


I(U) 

"^itn 


(280) 


in  which 

g  =  32.  174  ft/sec2 
is  the  acceleration  due  to  gravity,  and 


G(U) 


UKD  (U/Vso) 
1883 


(281) 


K  is  the  zero-yaw  drag  coefficient,  and  V  =  1,116.45  ft/sec  is  the 
lD  so  ' 

speed  of  sound  at  sea  level.  The  numerical  integration  is  carried  out 
backwards  (i.e. ,  the  integration  proceeds  from  larger  to  smaller  U) 
from  some  arbitrarily  chosen  value  of  U,  e.  g.  ,  V,  such  that  V>  u/aQ 
for  all  u  of  interest.  Initial  conditions  are  S(V)  =  0,  T(V)  =  0,  A(V)  =  0, 
and  I(V)  —  0.  The  solution  given  by  Eqs.  (274),  275),  and  (276)  is  inde¬ 
pendent  of  the  choice  of  V.  The  Siacci  functions  bear  the  following  names: 


S  Space  function 

T  Time  function 

I  Inclination  function 

A  Altitude  function 


At  the  date  of  writing,  tables  have  been  prepared  for  the  7.  62-mm 
NATO  round,  the  20-mm,  M56  round,  and  the  40-mm,  Mk  2  round,  (Refs.  62, 
63,  and  64). 

For  computer  applications  (airborne  or  ground  based)  the  Siacci 
functions  can  be  curve  fitted. 

The  Siacci  method  in  its  basic  form  is  useful  in  most  situations 
where  point-mass  (particle)  trajectories  are  applicable,  that  is,  situa¬ 
tions  where  the  projectile  angle  of  attack  (yaw)  is  small  and  all  aero¬ 
dynamic  forces  except  the  drag  are  negligible.  This  situation  occurs 
for  forward  fire  from  fixed-gun  fighters  and  sometimes  in  tail  defense 
of  bombers.  For  example,  the  Siacci  method  is  applicable  to  fire- 
control  problems  for  the  F-lll  and  for  use  with  tracer-line-type  gun 


114 


sights,  although  it  may  be  just  as  convenient  to  perform  onboard  particle 
trajectory  computations  in  these  instances. 


For  fire  from  flexible  guns  and  from  fixed  gun  fighters  in  high  angle- 
of-attack  situations,  the  angle  of  attack  of  the  projectile  can  be  large  and 
other  aerodynamic  forces  besides  the  drag  become  important.  In  this 
instance  the  Siacci  method  in  its  basic  form  becomes  inadequate.  On¬ 
board  trajectory  computations  are  no  longer  reasonable,  however,  when 
six-degree-of -freedom  computations  are  necessary.  Techniques  for  im¬ 
proving  the  accuracy  of  the  Siacci  method  for  the  large-yaw  situation  are 
given  in  Subsection  3. 

2.  The  Basic  Siacci  Method 

The  basic  Siacci  trajectory  is  an  approximation  of  a  point-mass 
trajectory.  Corrections  can  be  added  to  improve  its  accuracy  as  will 
be  discussed  later  m  Subsection  3.  The  equations  of  motion  of  a  mass 
point,  as  can  be  derived  from  Fig.  16,  are 


mx  =  -Ex 


(282) 


my  =  -E  y  -  mg 


(283) 


where 


E  ^  pd^  VKd(V/Vs) 


(284) 


V  is  the  oroiectile  velocity,  d  is  the  projectile  diameter,  m  is  the  pro¬ 
jectile  mass,  g  is  the  acceleration  due  to  gravity,  is  the  drag  coef¬ 

ficient,  and  p  and  V  are  respectively,  the  air  density  and  the  speed  of 

s 

sound  at  the  altitude  of  the  firing  point.  Both  p  and  Vg  are  assumed  con¬ 
stant  over  the  trajectory.  This  is  not  a  bad  assumption  for  the  air-to- 
air  case  and  for  close  air  support. 

Coordinates  x  and  y  are  to  be  replaced  by  Siacci  coordinates  P 
and  Q  shown  in  Fig.  16.  Substitution  of 


x  =  P  cos  0 


(285) 
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(286) 

(287) 

(288) 


(289) 

(290) 

(291) 


V  =  / P2  -  2PQ  sin  e  +  Q2  (292) 

N  /  O 

*  • 

For  short,  flat  trajectories,  Q«  P,  and, to  a  good  approximation, 
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V  *  P  U 


(293) 


This  is  the  Siacci  approximation  which  enables  the  trajectory  computation 
to  be  reduced  to  tabular  form  since,  from  Eq.  (289), 


dt 

mdu 

P  d2  u2Kd(u/  Vs) 

(294) 

can  be  integrated  directly. 

Before  solution,  it  is 

tions: 

desirable  to  make  the  following 

substitu- 

t  = 

C  [t  -  T  1 

*ao  L  °J 

(295) 

P  = 

(296) 

Q  = 

(297) 

u  = 

a  U 
o 

(298) 

where  C,  cr,  a  ,  T  ,  and  S  are  constants  and  T,  S,  q,  and  U  are  new 
o  o  o 

variables  which  replace  t,  P,  Q,  and  u,  respectively.  C  is  the  ballistic 
coefficient  given  by 

C  =  (299) 

(12d; 

2 

Units  of  C  are  lb/in.,  (for  consistency  with  past  usage)  when  m  is  in 

pounds  and  d  is  in  feet;  <r  is  the  ratio  of  the  air  density  at  altitude,  p, 

to  that  at  sea  level,  p  , 

o 

o-  =  -£—  (300) 

^o 
3 

where  pQ  =  0.076474  lb/ft  .  The  parameter  aQ  is  the  ratio  of  the  speed 
of  sound  V  ,at  altitude,  to  that  at  sea  level,  V  =  1 ,  1 16.  45  ft/ sec , 


a  =  TT 

O  V 


(301) 


With  these  substitutions,  the  equations  of  motion  become 


=  -  UG<U) 


(302) 


(303) 


(~tT  -  t! 


(304) 


where 


G(U)  = 


UK_(U/V  ) 
D'  '  so 


(305) 


The  reason  for  the  substitution  is  now  clear;  these  equations  are  indepen¬ 
dent  of  Vg  and  p  and  hence,  the  Siacci  tables  are  independent  of  altitude. 

From  the  first  equation 


dT  =  - 


auim 


and  hence 

U 

o 

m  ^  r  du 

T  ■  To  "  J  UGTdT 

u 

where  U  =  u  /aQ,  and  u  is  the  initial  projectile  velocity.  If  T  is  defined 


V 

C  du 


T<u>  •  j  mu] 

U 


(306) 
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where  V,  for  all  applications  of  interest,  is  any  convenient  number  such 
that  V  >  U  ,  it  follows  that 


T  -  T  =  T(U)  -  T(U  ) 
o  o 

V  V 

f  dU  C  dU 

“  3  UG(U)  “  j  UGlUJ 

U  U 


That  T  -  Tq  is  independent  of  V  ~s  obvious. 

The  equation  for  S  may  be  solved  in  a  similar  manner  by  eliminat¬ 
ing  dT. 


S  is  defined  as 


dS  =  U  dT  = 


S  -  S 

o 


U 


dU 

G(U) 


dU 

mi 


V 

<■* 


S(U)=\ 

J 

u 


dU 

mi 


therefore 


(307) 


S  -  S  =  S(U)  -  S(U  ) 
o  o 

For  the  thud  expression,  we  have 


gdU 

U*G(U) 


119 


1  dq 

tt  -at 


U2G(U) 


since  dq/dT  =  0  at  T  =  T^.  If  a  function  I(U)  is  defined  as 


v 

'Is 


'G(U) 


(308) 


we  have 


Thus 


tt4?  =  ««)-««„) 


dq  =  UI(U)dT  -  UI(UQ)dT 
I(U)  dU  A  I(Uo)  dU 

—  -  1  /rfr'  '  - a.  if  f  i 


I(^dU  -  I(U J  ds 


and  it  follows  that 


I(U)  dU 
'  G(U) 


-  i(uo)  (s  -  so) 


If  a  function 


U0 


.vr; 


V 

A(U)  .  ]" 

u 


I(U)  dU 

G(U) 


(309) 


is  defined,  then 


q  =  A(U)  -  A(U  )  -  I(U  )  |s(U)  -  S(U  )\ 
o  o  (  o  | 


It  follows  that  the  relations  for  t,  P,  and  Q  are 


t  = 


—  [■' 

ora 

o  L 


<ra.  ,  -  T< 


u  /a  )1 
o'  o 


(310) 


P  =  ~  [s(u/a  )  -  S(u  / a  )1 

cr  o  o  o  J 


(311) 


Q  = 


(~\  |"a(u/ a  )  -  A(u  /a  )  -I(u  /a  )  f-  p] 

1  cra^y  I  o  o  o  o  o  C  I 


(312) 


where  T,  S,  I,  and  A  are  given  by  Eqs,  (306)  through  (309),  respectively. 
In  application,  the  functions  S,  T,  A,  and  I  are  tabulated  once  and  for 
all  by  numerical  integration.  Given  initial  and  final  bullet  velocities  uq 
and  u,  respectively,  t,  P,  and  Q  can  be  calculated  by  use  of  the  tables. 

3.  Corrections  to  the  Siacci  Method 

This  section  contains  techniques  for  improving  the  accuracy  of  the 
Siacci  method  when  some  of  the  basic  assumptions  upon  which  the  method 
is  founded  are  violated.  Corrections  are  developed  for  a  variable  atmos¬ 
phere,  yaw  drag,  and  windage  jump.  This  section  is  not  complete  in  the 
sense  that  the  methods  are  fully  developed.  Rather  it  is  meant  to  serve 
as  an  indication  of  how  correction  techniques  are  derived.  Indeed,  tech¬ 
niques  depend  somewhat  upon  individual  rounds  and  are  best  investigated 
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in  conjunction  with  other  calculational  methods  such  as  six-degree-of- 
freedom  calculations.  Correction  terms  needed  for  one  type  of  pro¬ 
jectile  may  not  be  needed  for  another.  Also,  correction  terms  may  be 
derived  by  a  number  of  different  methods.  Thus,  different  formulae  for 
the  correction  terms  derived  herein  can  be  found  in  Refs.  65  and  66,  A 
comparison  of  correction  techniques  is  given  in  Ref.  67  for  some  ob¬ 
solete  rounds. 

3.  1  A  Variable  Air-Density  Correction  -  The  Siacci  theory  in  its 
basic  form  is  founded  upon  the  approximation  that  the  air  density  is  coii- 
stant  along  the  projectile  trajectory.  To  compensate  for  this  approxima¬ 
tion,  correction  formulae  may  be  derived.  The  following  comments  are 
pertinent,  however. 

In  problems  involving  air-to-air  combat,  the  target  and  attack 
aircraft  are  usually  at  about  the  same  altitude  and  the  air-density  cor¬ 
rection  is  unnecessary.  In  low-altitude  air-to-ground  problems,  ex¬ 
perience  has  shown  that  the  calcul  ted  hit  point  is  not  much  affected  by 
the  constant  air-density  approximation,  whereas  the  calculated  time  of 
flight  is  affected.  But,  in  air-to-ground  situatiors,  the  time  of  flight  is 
usually  not  needed  to  any  great  accuracy  provided  the  target  is  not  mov¬ 
ing  rapidly  and  the  wind  is  not  strong.  It  is  concluded  that  in  most  cases 
an  air-density  correction  is  not  needed,  but  nevertheless  correction 
formulae  will  be  developed  for  situations  where  they  may  be  needed.  The 
derivation  follows. 

It  will  be  assumed  that  the  expression 


P 


p  or  e 
r  o 


-hAy 


(313) 


adequately  represents  the  air-density  variation  at  points  between  the  gun 

and  the  target,  where  p  is  sea  level  air  density,  <r  is  the  relative  air 

■**3  “1 

density  at  the  gun,  h  is  a  constant  (h  =  3.  158  X  10  ft  ),  and  Ay  is 
the  altitude  variation  along  the  trajectory.  Evidently,  from  Fig.  16, 


Ay  =  P  sin  0  -  Q  ^  P  sin  0 

o  o 

Now  consider  the  unvawed  equations  of  motion  (289),  (290),  and  (291). 
Eq.  (289)  becomes 


-hP  sin  £ 


mu  =  -  p  ere 
r  o 


°  d2u2KD(u/Vs) 


under  the  Siacci  approximation  u  =  V.  It  follows  that 


C  dU 

o-a  UG(U) 
o 


-  hP  sin  0 


and  that 


c  du  -hr-  sin  0O 

—  cw  =  e  d 


Integration  of  the  second  expression  leads  to 


-9-  js(U)  -  S(UQ)J  =  -L: 


-hP  sin  0 


h  sin  0 


(314) 


-hP  sin  0 


Solving  for  e 


yields 


-hP  sin  0 


=  l _ h  sin  0  |S(U)  -  S(JQ) 


This  expression  is  of  the  form 


Hence 


-x  , 
e  =  1  -  y 


eX  =  jly  =  1  +  y  +  y2  +  yJ  + 


provided  [  j  <1.  For  an  air-to-ground  attack  from  an  altitude  of,  say, 


P  sin  0  =  5,  000  ft 

o 


x  =  hP  sin  0 

=  (3.  158  x  10'5)  (5,  000)  =  0.  16 


It  follows  that 


e  x  =  0.852 


y  =  0.  148 


y  =  0.0219 


Then,  to  within  about  2  percent 


e  -  1+y 


hP  sin  0 


=  x  +  -£_  hsiu  0  [  S(U)  -  S(U  )1 

cr  o  L  o  J 


The  expression  for  dt  becomes 


hP  sin  0 


dt  =  -  e 


o  C  dU 

<ra  UG(U) 
o 


C  h  sin  8o  [S(U)  -  S(Uo)]  j 


and  hence 


t 


c 

<rao 


[T 


(U) 


h  sin  0 


|[H(U) 


h(Uo) 


]  -S(Uo)  [T(U)-T(Uo)]| 


(315) 


where 

H<U»  -  j  ^  (316, 

u 

H(U)  may  be  tabulated  once  and  for  all  along  with  S,  T,  I,  and  A. 

The  numerical  example  above  applies  to  maximum  P  and  hence^to 
minimum  U.  As  U  runs  from  Uq  to  U,  y  varies  from  zero  to  0.  148. 

This  implies  that  the  integration  is  probably  much  more  accurate  than 
2  percent.  If  accuracy  is  found  insufficient,  e.  g.  ,  for  higher  altitude 
air-to-ground  fire,  more  terms  can  be  taken  in  the  expansion  for  (1  -  y) 

The  analysis  proceeds  in  a  straightforward  manner. 

The  correction  term  for  Q  may  be  derived  as  follows: 

d(  A  )  =  £.  dt 

'  u  r 


a  U 
o 


1  + 


Ah.i„eo[s,u,-S(uo>]|_^  ^ 


Then 


Q  _  C 


u  2 

<ra 

o 


[l(U,  -  l(Uo)] 


(317) 


h  sineo  [jS(U!  -  S(Uo>(  I(U)  -  jA(U)  -  A(Uo>j] 
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If,  for  short,  Q/u  =  F,  then 


dQ  =  Fa  dt 


F  <  1  +  —  h  sin  0 
I  O'  o 


[S(U)  -  3(Uo)]| 


C  dU 

T  CRXTT 


If  the  second-order  term  in  h  sin  is  discarded,  since  it  is  obviously 


small,  then 


dQ  = 


-(ttJ 

3 

-/■£)  aohsin94<S|U)  -S(Uo>}I<u>  -)A(U)  -  A>uo,(]<$rr 

-  (£~f  aohSin.u[l(U)  -  I(Uo>]  [»•!)  -  S(0o)] 


and  finally, 


Q  = 


(^t)2[|A(U)  '  A(U°’j ■  S(uo’!] 

+  (^)  ‘ohsineo[jW1lU)-Wl<Uo>| 


-  2S(Uo)  J  A(U)  -  A<Uq)  |  -  l(Uo)  jW2(U)  -  W2(Uo) 


+  JA(Uq)  +  I(Uq)  S(Uq)[  S(U)  -  S(Uc 


(318) 


where 


Vu»=i' 


2S(U)  I(U)  -  A(U) 


(319) 
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W2<U>  =  I  CM  dU  <320) 

u 

The  correction  term  is  quite  complicated,  and  hopefully,  is  negligible 
in  most  cases. 

3.2  A  Yaw-Drag  Correction  -  In  more  exact  trajectory  computations, 
other  aerodynamic  forces  in  addition  to  the  zero-yaw  drag  force  are  in¬ 
cluded,  and  the  projectile  is  treated  as  a  rigid  body  which  executes  angular 
motion  about  its  center  of  mass.  An  important  additional  force  is  the  yaw 
drag  which  arises  from  the  angle  of  attack  (yaw)  of  the  projectile.  The 
yaw  angle  6  is  measured  between  the  projectile  body  axis  and  the  projectile 
velocity  vector.  The  drag  force  for  the  yawed  projectile  can  be  written 
as 


Kd(M,6) 


kd  (M) 

o 


1  +  K 


V 


(321) 


where,  for  present  purposes,  Kn  ?  is  considered  to  be  a  constant  and 


M  =  V/Vg  is  the  Mach  number.  This  expression  for  takes  the  place 
of  Kj-j  in  the  equations  of  motion,  i.  e.,  inEqs.  (2 8Z);  (283),  and  (284), 
and  also  inEqs.  (287)  and  (288).  If  E  retains  its  meaning,  Eqs.  (287)  and 
(288)  become 


mP  =  -E  I  1  +  K 


[1  +  Kvf2^ 


mQ  =  -E  1  +  K 


1  +  KD  ?  62j  Q  +  mg 


If  the  derivation  proceeds  as  before  with  the  Siacci  approximation  V  -  P, 
Eqs.  (289),  (290),  and  (291)  show  that  the  following  relations  will  be 


obtained: 
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•  j  i. 

mu  =  -  p  d 


“2  Kelts' [>  t^S2! 


P  =  u 


d 

St 


(4)-f 


Substitution  of  aQU  for  u,  pQ<r  for  pt  144d  C  for  m,  and  use  of  Eq.  (305) 
for  Kj-j  yields  the  following  result  for  the  first  relation  above 

C  dU 


1  ♦  K  62' 
.  o 


dt  =  - 


0- aQ  UG(tF) 


It  follows  immediately  that 


r 

o-  a 

o  L 


T(u/ao) -  T( 


t 

V»0>1  -  KDr2  I  62  dt 
J  ®  0 


(322) 


With 


dt  « 

u 


dP 

a~TT 

o 


it  follows  that 


[‘ +  Kv  6Z] 


1  +  kDs2  62J  dP  =  -•£-  ^ 


and 


P  = 


[s(u/ao)  -  S,uo/ao)]  -  KDf  l  j"  62  dP 


(323) 


In  order  to  derive  a  correction  term  for  Q  it  is  convenient  to  write 
the  differential  equation  for  Q  in  a  different  form  as  follows 
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or 


Let 


so  that 


Then 


or 


d 

cTT 


fo\  =  ±(2\ 

d 

dQ 

dP 

d 

MQ\ 

w  dtUJ 

1  dt 

dP 

dt 

dP 

Up) 

=  u 


d2Q 

dP2 


D  = 


d2Q 

dP2 


_£._ 

2 

u 


dQ 

dP 


1 

u 


dD  _  g 

dP  “  “T 
u 


[i+v62]dD  =MltKv62]dP 


g  £  dU 

=  ’  2  cr  GltTT 


C  gdU 


dD  *  -V62  ^dP’  „*Glm 

o  u  era  U  G(U) 


and  it  follows  that 


D  -  -KD r2  I  6 


2  -h  dP  +  -%  [««>  -  IIU„I| 

o  ~  u  era  L  J 

0  o 


(324) 


129 


Also 


[‘  +  *v] 


dQ  =  D|wKDg2|dP 


-  -v[l  +  Kv1PJ 


dP'  dP 


~r  [1(U)  -  I(Uo»]  f  wit 

o 


It  follows  that 


Q  =  -  KD  y  a2  dQ  -  K  J  [l  ♦  K  a2]  j  62  \  dP" 
13  0  0  0  L  0 

[|A<U)  -  a<uo»|  -«u.»|s‘u»  -S(uo'|] 


D  .2  2 

o  cr  a 

o 


I(u  /a  )  f  52  dP 
'00  J 


Q 

-K  j  52dQ 


(^r)tA(u/a° 


)  -  A(uo/ao)  -  I(u0/a0)  £  pJ 


dP"  dP’ 


(325) 
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*■' •re~trr^A*y  w  „  ,. 


It  is  seen  that  equations  for  t,  P,  and  Q  can  be  written  as  before 
but  with  correction  terms 


Atl  = 


AP, 


t 

-KDfi2  I  62  dt 
6  0 

p 

-kd »  £  **** 
6  0 


(326) 


(327) 


AQ1  = 


-KDs2  -m  «  W  J  6'  dP 

6  ,ao  0 


Q 

[D  2  J  6 


,  J  6  dQ-K 
0 


°62 


j  [1+Klv62]jV4 


dP"dP' 


(328) 


These  equations  can  be  evaluated  only  if  6  is  a  known  function. 
According  to  an  approximate  analysis  by  Sterne  (Ref.  9)*  6  is  given 
approximately  by  the  relation 


62  =  62  S°  "  VZ  e-2(rcP  (329) 

o  s  -  1 
o 

where  6  is  an  average, squared  yaw,  represents  the  static  stability 
factor  which  is  given  by 


s 

o 


A2N2 


4Bpd3u2KM 


(330) 


and 


c  =  c 1  + 


(331) 
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where  N  is  the  axial  spin, 


and 


c 


ii 


(332) 


(333) 


When  Sterne's  analysis  was  derived,  there  was  no  convenient 
method  of  solving  the  general  six-degree-of-freedom  equations  since 
there  were  no  large-scale  digital  computers.  Now,  it  may  be  convenient 
to  solve  them  for  yaw  dependence  and  curve  fit  the  solution.  In  any  case, 
it  is  probably  true  that  7)  is  given  approximately  by 

Z  =  K6Qe'|3P  (334) 

where  K  and  (3  are  constants. 

From  the  approximate  theory,  the  equation  for  yaw,  Eq.  (272), 
may  be  written  as 


where  the  gravity  term  and  the  damping  term  have  been  deleted  and 
the  approximation  sin  1)  =  1)  has  been  used.  a.  Eq,  (3j4)  is  assumed 
for  6,  then 


%  =  -  (3  PK6q  e'pP  =  -pu6 

and  it  follows  that 
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it  rs  interesting  to  compare  this  result  with  that  obtained  from  the  Sterne 

theory,  Eqs.  (329)  through  (333),  in  which  p  =  crc.  Curiously  enough, it 

can  be  shown  that  the  two  are  in  reasonable  numerical  agreement,  at 

least  for  the  20-mm,  M56  round,  despite  the  different  emphasis  on  KTT  in 

ri 

the  two  theories. 

From  examination  of  Eqs.  (326),  (327),  and  (328),  it  is  evident  i 
only  AP^  can  be  integrated  directlv: 


?  =  -1/2  r 

=  -kds26„  -r- T-  J 

6  o  Q 


-20-CP.T3  o 

e  dP  =  -  — 

2crc 


[l  -  e2'cP] 


(335; 


where 


s  -  1/2  , 

k  =  -4 - T-  Kn  7  6 

o  s  -  1  u  2  o 

o  6 


The  equation  for  S(u/aQ)  is  often  written  as 


(336) 


S(u/a  )  =  S(u  /a  )  +  -fP  +  (337) 

o  o  o  C  2c  C 

2(r  c  P 

in  which  the  e  term  is  neglected  for  a  large  enough  P. 

An  expression  for  At^  can  be  derived  as  follows:  since  P  =  u, 
and  with  k  =  Zac, 

At  =  -  K 


1 


7  '  */z 

6  - 5 -  T  =  -k  T 

D  2  o  s  -  1  o 

5  o 


where 


f  • 

J  ^ 
0 


-kP 


a  t  \  — r* —  d  P 


The  expression  for  S(U)  with  the  AP^  correction  is 


S(U)=S(Uo)+^P+^.J,  - 


■  2<rcP 


(338) 
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By  use  of  this  expression,  a  series  in  P  can  be  derived  for  1/U.  It  is 


TT"  tr-  +  -^2—  j(1  +  ko>  rr  pj 


G'u0) 

2U2 


2cCk 


yf-  G(U  )  -  G'(u  ) - -2L- 

Uo  °  0  (1  +  k  )2 


|u+ko)rp| 


where 


G'(U  )  = 


dG(U  ) 
o 

TO - 


If  only  the  first  two  terms  of  this  series  are  used, 

P 
a 


V“  I 


1 


+  — 
o  U 


G(U  ) 

°-  (1  +  k  )  -£-P 


2  ’  "o'  TT 


-kP 

e  dP 


The  function  represented  by  the  first  two  terms  of  the  series  for  1/U  is 
a  straight  line  with  the  correct  magnitude  and  slope  at  P  =  0.  As  P 
increases,  it  deviates  from  1/U.  But  as  P  increases,  the  exponential 
term  decreases,  so  the  major  contribution  to  the  integral  is  obtained  for 
P  near  zero  where  the  approximation  to  1/U  is  more  accurate.  It 
follows  that 


at 


1 


-kP 

I  -  e 
ku 

o 


+ 


<1  +  ko} 

kV 

o 


-  {1  +  kP)  e 


If  this  expression  is  not  accurate  enough,  more  terms  can  be  retained. 
For  large  kP,  it  reduces  to 
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1  ,  (1  +  ko|G(uo/aol  ,rac 

atl  -  -ko  T3~  - 7T~2 - C~ 

o  k  a 

o 


(339) 


An  expression  can  be  derived  in  a  similar  manner  for  AQ^.  By 
use  of  Eqs.  (329)  and  (336),  Eq.  (328)  becomes 


C  i  lrP 

AQ,  =  -  -V  I(-„/a  )  k  \  e  **  dP 

1  ra2  °  °  0  J 

o  0 


-  ko  J  e  "kP  -3F  dP  -  ko  ff1  +  ko  e‘kP'l  I e_kP"  ^  dP"  dP' 

0  0  L  J0  u 


LetZ^,  Z^,  and  Z^  represent,  respectively,  the  three  integrals.  The 
first  is 


Z1  =  j  e~^  dP  = 


1  -  e 


The  second  may  be  integrated  by  parts  and  by  use  of  the  expression 


d  Q  =  _g_ 


the  following  expression  is  obtained 


a 

„  _  dQ  e~'  f  g  e” n‘ 

Z2  "  "  "cTF  +  J  2  ~ 

0  U 


where,  as  follows  from  Eq.  (324) 


P 

ar- =  7^2  ["I<u/ao'  -  I<uo/ao']'  ko  I' 
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The  third  integral  Z  may  be  integrated  by  parts;  it  is 


'  [p  "ir  °  P]  I e  kP  ~r  dp 


-I 

II  -' 


-2kP  g 
e  -VdP 
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With 


Z4  =  j  \  dP 
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Z5=  IPe_kP3 


■  I  *"2kP  5 dP 


it  follows  that 


‘1  “  k 


*-%  j 

o  v 


I{u/ao)e_kP  -  I(uo/ao)j 


-d+kp)  Z4  +  kZ5-koZ6 


Z^,  Z^,  and  Z^  can  be  evaluated  approximately  by  expanding  1/U  in 
series.  The  first  two  terms  of  this  series  are 


1=1+  2  G(U  )(1  +  k  )  Z-  P  + 

U2  U2  U3  °  °  ^ 

o  o 
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-■■* ■■'. - TSKZl — prr 


It  follows  that 


(340) 


2 

More  terms  in  the  series  for  1/U  may  be  retained  if  necessary.  On  the 
other  hand,  one  would  expect  to  drop  some  of  the  exponential  terms  when 
P  is  large. 

3.  3  A  Windage-Jump  Correction  -  A  spinning  projectile  fired  into 
a  crosswind  experiences  an  angular  deflection  of  its  direction  of  motion 
out  of  its  initial  plane  of  yaw.  This  deflection  is  numerically  equal  to 
about  5  percent  of  the  initial  angle  of  attack  (yaw).  For  an  initial  yaw  of 
;°  =  17.  5  milliradians,  the  windage  jump  is  about  0.  87  milliradians ,  or 
for  10°  yaw  it  is  about  8.7  milliradians.  For  the  20-mm,  M56  round,  the 
windage  jump  in  milliradians  is,  as  a  rule  of  thumb,  numerically  equal 
to  the  yaw  angle  in  degrees  (e.  g.  ,  5°  yaw,  5  milliradians  windage  jump). 
For  forward  fire  from  a  fixed  gun  in  an  aircraft,  windage  jump  is  usually 
ignored  since  the  dispersion  of  such  systems  is  of  the  order  of  5  milli¬ 
radians.  For  a  flexible,  side -firing,  gun  system,  the  windage  jump  must 
be  included,  however,  since  the  initial  yaw  may  be  as  high  as  20°. 

For  aircraft  fire,  the  initial  plane  of  yaw  is  the  plane  containing 
the  velocity  vector  of  the  aircraft  and  the  muzzle  velocity  of  the  gun 
VM.  The  initial  yaw  angle  6q  is  the  angle  between  VM  and  the  projectile 
initial  velocity  vector  uq,  where 
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If  the  angle  between  and  Vj^  (the  gun  angle)  is  A,  then,  as  can  be 
seen  from  Fig.  17,  6q  is  given  by  the  relation 


6  ~  sin  6 

o  o 


sin  A 
u 

o 


(341) 


where 


V  2  +  V  2 
VA  M 


+  2Vavm  cosA 


(342) 


Figure  17 

Projectile  Geometry  Showing  Initial  Angle  of  Attack 

6  in  the  Initial  Plane  of  Yaw 
o 


"  ~  -7  -  >  -v-  75* 


The  magnitude  of  the  windage  jump  in  radians  is  given  approximately 
in  Ref.  9  as 


e 


b6 
_ o 

u 

o 


(343) 


where 


.  AN 
b  "  ^ 


(344) 


A  is  the  axial  moment  of  inertia,  N  is  the  spin,  m  is  the  mass,  d  is  the 
diameter,  is  the  lift  coefficient,  and  is  the  overturning  moment 
coefficient.  (Note  that  this  expression  does  not  contain  the  Magnus 
moment  coefficients  Kp  and  K^,.) 

In  calculations,  the  windage  jump  is  usually  treated  as  a  small 

correction  to  u  ,  in  which  case  u  is  defined  as 
o  o 


u  =  VA  +  V,  +  J  (345) 

o  A  M 

where  J  =  uq€  and  J  is  in  tne  direction  of  VA  X  V^j.  F rom  Eq.  (341 ), 

Eq.  (343),  and  since  I  ^A  X  VM  |  =  VAVM  sin  A>  follows  that 


J 


(?A  X  hj 

IVhJ 


(346) 


Since  J  is  small,  the  magnitude  of  uq  is  not  changed  significantly.  The 
direction  of  P  is  taken  to  be  along  uq,  so 


p  =  — 2-  P  (347) 

u 

o 


An  alternate  treatment  of  the  windage  jump  can  be  based  upon  the  set 
of  approximate  equations  of  motion  of  Section  IV,  or  on  any  appropriate 
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six-degree-of-f reedom  code.  The  Eglin  code  R370,  based  upon  the 
approximate  equations,  is  used  to  calculate  windage -jump  components 
(among  other  things)  for  the  20-mm,  M56  round  in  and  perpendicular  to  the 
initial  plane  of  yaw.  These  components  may  be  used  instead  of  J  above. 
Experience  indicates  that  these  windage -jump  components  are,  for  large 
enough  P,  very  nearly  independent  of  all  initial  condition  parameters 
except  6  .  Appropriate  formulae  and  necessary  explanations  follow. 

Windage-jump  parameters  are  J 2  and  J3»  J3  corresponds  to  j 

above  and  is  normal  to  the  plane  of  yaw,  whereas  J2  is  in  the  plane  of  yaw 

and  is  zero  in  the  elementary  theory.  J2  and  J3  are  output  in  milliradians 

as  a  function  of  time  but  they  approach  constant  values  as  time  increases. 

From  J  and  J  ,  swerve  components  in  the  units  of  P  (distance)  are  cal- 
2  3 

culated  from 


S2  =  10-3  J2P 

(348) 

S3  =  io-3j3p 

(349) 

These  components  are  transformed  into  an  £,q,t,  coordinate  system 
(Fig.  12)  by  means  of  the  equations  (from  Section  IV) 

S;  =  -  (S  cos  <|>  -  S3  sin  4>q)  sin  0q 

(350) 

=  (S2  cos  «|>o  -  S3  sin  4>q)  cos  0q 

(351) 

S;  =  S2  Sin  ^o  +  S3  C°S 

(352) 

where  0  is  the  elevation  angle  of  uq  above  the  horizontal  and  <j>Q  is  the 
angle  measured  about  uq  between  the  vertical  plane  and  the  initial  posi¬ 
tion  of  the  plane  of  yaw.  Components  of  q,  and  £,  are  given  by 

£  =  P  cos  0  +  S 

5  o  i 

(353) 
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(354) 


ri  =  P  sin  0  -  Q  +  S 

o  T| 


l  =  s. 


(355) 


The  £,r|,£  coordinate  system  is  shown  in  Fig.  12;  £  and  t,  are  horizontal, 
and  t|  is  vertical  and  is  measured  up.  The  system  is  right-handed,  so  t, 
is  measured  out  of  the  page.  P,  Q,  and  uq  are  in  the  | , tj  plane.  Expres¬ 
sions  for  0  and  d>  are  derived  in  Section  VI. 
o  To 

Values  of  and  used  in  20-mm,  M56  projectile  calculations  are 


J,  =  56 
2  o 


(356) 


J,  =  53  6 
3  o 

where  6  is  in  radians,  and  J_  and  J_  are  in  milliradians. 
o  2  3 


(357) 


3.4  Independence  of  Correction  Terms  -  The  yaw-drag  and  the  var¬ 
iable  air-density  corrections  were  derived  as  if  the  two  effects  aie  inde¬ 
pendent.  This  has  some  physical  justification  since  the  yaw-drag  effect 
is  important  at  the  beginning  of  the  trajectory, whereas  the  change  in  air 
density  is  greatest (e. g. ,  in  the  air-to-ground  case)  at  longer  ranges. 
There  is  a  difficulty,  however,  with  the  equation  containing  P.  The  ques¬ 
tion  arises  as  to  whether  Eq.  (314)  or  Eq.  (337)  should  be  used.  This  is 
resolved  by  including  both  effects  in  Eq.  (289).  The  following  relation  is 
obtained 


C  dU 
<r  ‘CJlUl 


—  e 


|/hPSine°tkoe-^]d 


This  is  a  good  approximation  since  in  practice  k  »  h.  (For  the  20-mm, 
M56  round,  k  «  0.  004,  whereas  h  0.  00003. )  It  follows  that  the  cor¬ 
rect  expression  to  be  used  must  be 
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-hP  sin  0 


-f  [S(u/V  -  S(“o/ao)]  -- 

!hP  sin  0  \  k 


«  P<1 


This  last  approximation  should  be  good  past  P  =*  1 , 000  ft. 


(358) 


3.  5  Summary  of  Results  -  Formulae  derived  in  this  section  are 
listed  below  for  convenience.  Q  and  t  denote  gravity  drop  and  time  of 
flight  as  obtained  from  the  basic  Siacci  theory,  whereas  and  t  denote 
corrected  values. 


k  f  hP  sin  0{ 

S(u/ao)  =  S(uo/aQ)  +  +  tT  P  }  1  1 


‘  =  -k;  [T<“/ao»  -  T<w] 


t  =  t  +  At  +  At, 
c  pi 


Q  =  Q  +  AQ„  +  AQ, 

C.  OX 


A(u/a  )  -  A(uo/aQ)  -  !(*/*<,>  P 


-  ■  fit 


aoh  sin  eo  |[H(u/ao)  '  H(uo/ao)] 

era  ) 
S(uo/ao)  -cr 1  | 


<r  a  1  +  k 

■  °  - - —  G(u  /a 

C  hu  o  c 
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AQ,  ■  -  TT  —z  W*o» 

ffao 


k 

I 

71 


oi_  ' 
I  2 
a 

A 


+  kP  -  2 


o-a  1  +  k  / 

- - SL  G(u  /a  )  J  1  - 

C  ku  '  o'  o'  ) 

°  ( 


-2.  -  kP 


H(U)  = 


V 

■f 

u 


S{U)  dU 

UGfuT 


AQ  = 
P 


— \  a  h  sin0 
(r  a  I  o  o 

°) 


Wi(u/ao)  -  Wl(uo/ao) 


'  2S(uo/ao)jA(u/ao)  -  A(uo/aQ)j 


-  | 


W2(u/ao)  -  W2(„o/ao) 


+  A(uo/ao)  +  I0(uo/a0) 


Wl(U)  =  j  du 


S{uo/ao)|  |S(u/ao)  -S(llo/ao,!j 


u 

V 


w2(u) 


f  S(U)  dU 

j  c.(ur 


u 


Exponential  terms  have  been  eliminated  from  the  correction  terms  for 
Atj  and  AQ^.  Additional  equations  (Sterne's  theory)  are 

k  =  2crc 


c  =  c '  + 


s  -  1 
o 
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c 1  = 


=fe-  !¥■  -.*« 


P  d 

1ST  Kd 


and  also 


A  N _ 

4B  p  d3u2  K.  , 
r  o  M 


s  -  1/2 

k  =  _2 - —  je  5^ 

o  s  -  1  D  2  o 
o  6 


h  =  3.  158  X  10'5  ft"1 


An  alternate  expression  for  c  derived  from  the  approximate  theory  is 


c 


] 


Equations  describing  the  windage  jump  are 


Alternate  windage-jump  equations  from  the  approximate  theory  (see 
Fig,  12)  are 
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^  c;  ■r^'^TST"  "•  '► 


s2  =  io-3j2p 

S3  =  10'3J3P 

=  -(S,  cos  <b  -  S„  sin  A  )  sin  9 
5  -  o  o  To  o 

ST1  =  {s2  COS  *0  ‘  s3  Sin  V  COS  °o 

S„  =  S.  sin  d>  +  S.  cos  A 
b  &  oi  To 

£  =  P  cos  G  +  S,. 

o  g 

ti  =  P  sin  0  -  Q  +  S 

O  t| 
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SECTION  VI 


AIRBORNE  FIRE  CONTROL  APPLICATIONS 


1 ,  General 

In  the  present  application,  the  gun-pointing  problem  is  complicated 
by  the  fact  that  both  the  gun  and  the  target  may  be  moving.  In  the  air-to- 
ground  application,  the  target  is  stationary  in  many  instances,  but  other¬ 
wise  its  motion  often  should  be  considered.  In  a  situation  where  the  pro¬ 
jectile  time  of  flight  is  1  sec,  for  example,  the  distance  moved  by  a  60-mph 
vehicle  is  88  ft.  The  problem  at  hand,  the  fire  control  problem,  involves 
predicting  the  motion  of  the  target  so  that  the  hit  position  may  be  obtained, 
and  involves  determination  of  the  correct  gun-pointing  direction  to  score 
a  hit  on  the  target.  Prediction  of  target  motion  is  known  as  kinematic 
prediction,  whereas  the  gun-pointing  problem  is  known  as  ballistic  pre¬ 
diction.  Kinematic  and  ballistic  prediction  are  discussed  briefly  in  Sub¬ 
section  2. 

In  Sections  IV  and  V,  equations  of  motion  were  developed  for  a 
projectile  in  flight,  but  initial  conditions  were  not  considered.  Initial 
conditions  will  be  developed  in  Subsection  4  for  the  equations  of  mo' ion 
developed  in  Sections  IV  and  V. 

Of  interest  in  current  applications  is  the  gatling  gun.  The  specifi¬ 
cation  of  initial  conditions  for  a  projectile  fired  from  such  a  weapon 
mounted  in  a  turret  in  a  moving  aircraft  Js  complicated.  The  combined 
motions  of  the  rotating  barrel  cluster,  the  turret,  and  the  aircraft  can 
cause  errors  in  calculations  if  they  are  not  accounted  for.  For  example, 
the  distance  from  the  aircraft  center  of  mass  to  the  gun  muzzle  may  be, 
say,  10  ft,  and  the  aircraft  angular  velocity  perpendicular  to  this  dis¬ 
tance  may  be,  say,  60  deg/sec  or  about  1  rad/sec.  The  component  of 
projectile  velocity  due  to  this  angular  rotation  is,  then,  10  ft  X  1  rad/sec 
=  10  ft/sec.  If  this  component  is  perpendicular  to  the  muzzle  velocity, 
which  is  about  3300  ft/sec  for  the  20-mm,  M56  round,  the  angular  error 
due  to  ignoring  this  angular  motion  is  10  ft/sec  divided  by  3300  ft/sec 
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and  this  equals  .003  rad  or  3  mr.  Slew  rates  of  the  gun  will  probably 
be  of  the  order  of  1  rad/ sec  and  the  distance  from  the  gimbals  to  the 
gun  muzzle  may  be  between  5  and  10  ft.  The  error  due  to  this  cause 
may  be  up  to  3  mr  also.  The  M6l  gatling  gun  fires  up  to  6000  rounds 
per  minute  from  six  barrels,  so  the  barrel  cluster  spins  at  1000  rota¬ 
tions  per  minute  or  2ir  X  1000/60  =  105  rad/sec.  The  distance  from  the 
barrel-cluster  axis  of  rotation  and  the  center  of  any  barrel  is  1.877  in. 

It  follows  that  the  rotating  barrel  cluster  imparts  a  velocity  of  105  rad/sec 
times  1.887/12  ft  equals  16.4  ft/sec  to  the  projectile.  This  amounts  to 
1000  X  16.4/3300  =  5  mr.  Evidently  these  effects  should  be  accounted 
for  in  initial  condition  calculations  and  Subsection  3  contains  the  neces¬ 
sary  coordinate  transformations. 

2.  Kinematic  and  Ballistic  Prediction 

The  fire  control  problem  may  be  divided  into  two  parts:  kinematic 
prediction  and  ballistic  prediction.  To  determine  the  future  target  path 
as  a  function  of  time  and  to  find  the  hit  position,  given  the  time  of  flight, 
is  the  kinematic  prediction  problem.  For  the  present  purposes,  it  will 
be  assumed  that  the  target  path  is  a  known  function  of  time.  Kinematic 
prediction  in  air-to-air  applications  is  the  subject  of  a  separate  study. 

The  determination  of  the  correct  gun-pointing  direction  to  score 
a  hit  on  the  target,  given  the  hit  position,  is  the  ballistic  prediction 
problem.  When  the  gun  and/or  the  target  is  moving,  hits  cannot  be 
scored  by  pointi  g  the  gun  directly  at  the  target  except  in  unusual  cir¬ 
cumstances.  The  correct  gun-pointing  direction  is  found  by  an  iterative 
trial-and-error  procedure.  A  first  guess  is  made  for  the  gun-pointing 
direction,  a  trajectory  is  calculated,  and  the  miss  'Ustance  is  used  to 
correct  the  gun-pointing  direction.  If  a  second  calculated  trajectory 
misses  the  target,  the  process  is  repeated.  The  iterative  procedure  is 
continued  until  the  gun  is  on  target. 

Solution  of  the  ballistic  prediction  problem,  as  described,  implies 
knowledge  of  the  hit  position,  which  in  turn  implies  knowledge  of  the 
v me  of  flight.  But  the  time  of  flight  to  ary  point  is  not  known  until  the 
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correct  gun-pointing  direction  to  score  a  hit  on  that  point  is  determined 
and  a  trajectory  Las  been  calculated.  In  other  words,  the  kinematic  and 
ballistic  problems  are  interdependent  and  cannot  be  solved  separately. 

An  iterative  solution  is  available:  a  first  estimate  to  the  time  of 
flight  is  chosen,  and,  given  the  target  path,  an  estimated  hit  position  is 
calculated.  Given  the  estimated  hit  position,  the  gun-pointing  direction 
to  score  a  hit  at  that  point  is  determined,  as  described  above,  and  the 
time  of  flight  is  calculated.  The  calculated  time  of  flight  and  the  first 
estimate  to  the  time  of  flight  are  used  to  correct  the  estimated  hit  posi¬ 
tion  and  the  process  is  repeated  until  the  correct  hit  position  is  obtained. 

The  iterative  solution  of  the  kinematic  and  ballistic  prediction  prob- 
’^ms,  as  described,  requires  the  calculation  of  several  trajectories. 

The  best  such  procedure  will  keep  the  required  number  of  trajectory  com¬ 
putations  to  a  minimum.  An  investigation  of  algorithms  for  onboard  kine¬ 
matic  and  ballistic  prediction  is  needed,  and  will  not  be  treated  here. 

The  purpose  of  this  subsection  is  to  point  out  the  existence  of  this  problem. 

3.  Coordinate  Systems  and  Transformations 

In  the  development  of  initial -condition  equations  for  a  projectile 
fired  from  a  turreted  gatling  gun,  five  right-handed,  rectilinear  coor¬ 
dinate  systems  will  be  used  and  they  are  as  follows:  (1)  an  earth-fixed, 
inertial  system,  S^;  (2)  a  system,  S^,  fixed  in  the  aircraft  with  its  origin 
at  the  center  of  mass,  with  the  x^  axis  along  the  body  longitudinal  axis 
and  directed  out  the  nose,  with  the  y^  axis  out  the  right  wing,  and  with 
the  axis  pointed  down  toward  the  aircraft  floor;  (3)  a  system,  S^,, 
fixed  in  the  aircraft  but  with  its  origin  on  the  axis  of  the  outer  turret 
gimbal;  (4)  a  system,  S^,,  attached  to  the  gun  frame  with  its  origin  on 
the  axis  of  the  inner  gimbal  and  which  rotates  about  the  two  gimbal  axes; 
and  (5)  a  system,  S^,  attached  to  the  barrel  cluster  with  its  origin  on 
the  gatling  gun  barrel-cluster  axis  of  rotation. 

Coordinate  transformations  between  these  systems  are  needed  so 

^  mm*  ■  A 

the  projectile  velocity,  uq,  angular  velocity,  wq,  position, Rq,  yaw,  6^, 
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and  precession  angle,  ^o,may  be  calculated  in  inertial  space,  Sj,  at 
firing  time.  To  calculate  6q  and  4>q,  the  direction  of  the  longitudinal 
axis  of  the  projectile,  i.e. ,  the  boreline  of  the  gun  barrel,  is  needed 
in  inertial  space. 

To  expedite  setting  up  these  transformations,  matrix  notation  will 
be  used.  The  appropriate  derivations  of  matrix  equations  are  given  in 
the  appendix.  A  summary  of  results  follows. 

3. 1  Matrix  Notation  -  The  results  of  the  appendix  may  be  sum¬ 
marized  as  follows:  The  transformation  of  the  coordinates  of  a  point  P 
as  observed  in  space  S'  into  the  coordinates  of  P  as  observed  in  space  S 
is  given  by  the  matrix  equation 

a  =  c  +  T  a' 

where  a'  represents  P  as  observed  in  S',  the  symbol  a  represents  P  in 
S,  c  is  the  position  of  the  origin  of  S'  as  observed  in  S,  and  T  is  a  3X3 
matrix  which  relates  the  angular  orientation  of  S'  to  that  of  S.  The 
transformation  of  the  velocity,  a',  of  P  as  observed  in  S'  into  the  velo¬ 
city,  a,  of  P  as  observed  in  S  is  given  by 

a  =  c  +  T(a'  +  ft' a') 

where  c  is  the  velocity  of  the  origin  of  S'  as  observed  in  S,  and  ft'  is 
the  3X3  matrix  representation  of  the  angular  velocity  of  S'  relative 
to  S.  ft'  may  be  written  as 


/  ° 


0 


ft 
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where  R  ,,  R  ,,  and  R  ,  are  the  components  of  angular  velocity  of  the 
x  y  z 

S'  system  with  respect  to  the  S  system  as  measured  in  the  S'  system. 

R'  may  be  transformed  into  the  S  system  by  the  relation 

R  =  TR'  T 

where  ¥  is  the  transpose  of  T.  The  column  vector  representation  of  R' 
is 


however,  and  the  transformation  from  S'  to  S  is 

u  =  Tu> 

The  column  vector  notation  is  obviously  easier  to  use  in  coordinate 
transformations . 

A  proof  that  angular  velocities  add  is  also  included  in  the  appendix. 
For  example,  if  oj'  is  the  angular  velocity  of  S'  with  respect  to  S  as 
measured  in  S',  and  w*  is  the  angular  velocity  of  a  solid  body  measured 
in  S',  then 

(jj  =  T(oj'  +  w*) 

is  the  angular  velocity  of  the  solid  body  measured  in  S. 

3.2  Subscript  Convention  -  It  will  be  convenient  to  employ  a  sub¬ 
script  notation  to  designate  coordinate  systems.  Thus,  the  coordinate 
space  J  is  designated  by  Sj  and  the  coordinate  space  K  is  S^.  The 
angular  orientation  of  the  coordinate  axes  of  Sj  with  respect  to  S^,  is 
given  by 


eJ  "  eK  TKJ 


or,  equivalently,  the  angular  orientation  of  the  coordinate  axes  of 
with  respect  to  Sj  is 

e  -  T  e 
K  KJ  J 

where  ej  and  e^,  are  row  matrices  composed  of  the  unit  vectors  directed 
along  the  axes  of  the  two  systems.  The  transformation  from  Sj  to  S^.  is 

bK  =  bKJ  +  TKJ  bJ 


where  bj  is  the  vector  (column  matrix)  representing  the  coordinates  of 
a  point  P  in  Sj  and  is  the  vector  representing  P  in  S^.  The  vector 
bKJ  rePresents  t^ie  distance,  measured  in  from  the  origin,  O^,  of 
to  the  origin,  Oj,  of  Sj. 

The  velocity  transformation  is  given  by 


VK  =  VKJ  +  TKJ  <VJ  +  “kjV 


where  VT  =  bT  and  V  =  b  represent  the  velocity  of  P  in  S  and  S  , 

J  J  J\  J  IV 

respectively.  V^j  =  is  the  velocity  of  Oj  with  respect  to  Oj^  as 

measured  in  S^,  and  is  the  angular  velocity  of  Sj  with  respect  to 

S..  as  measured  in  S  .  The  column  matrix  representation  of  is 
XV  J  IV  J 

“KJ* 

3.3  Coordinate  Transformations  -  Coordinate  systems  and  trans¬ 
formations  are  described,  starting  with  as  shown  m  Fig.  :  8.  There 
is  some  freedom  of  choice  in  the  way  these  systems  may  be  defined,  so 
definitions  were  made  as  convenient. 

In  Sg,  the  barrel-cluster  system,  the  x^  axis  is  along  the  barrel- 
cluster  rotation  axis,  and  the  gun  muzzle,  for  any  particular  barrel  at 
the  time  of  firing,  lies  on  the  y^  axis  a  distance  L  from  the  origin. 
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The  position  of  the  gun  muzzle  at  the  time  of  firing  is  given  by  the 
matrix 


and  the  direction  of  the  gunbore  axis  is  specified  by  the  matrix 


The  velocity  of  the  projectile  as  measured  in  Sg  is  parallel  to  Vg  and 
is  given  by 


UB  ~  VBVM 


(361) 


where  is  the  muzzle  velocity  (a  scalar),  and  the  angular  velocity  of 
the  projectile  as  measured  in  Sg  is 


Zn 

TkT  vmvb 


(362) 


where  n  is  the  distance  traveled  in  units  of  d  during  one  complete  rota¬ 
tion  of  the  projectile.  Units  of  n  are  calibers/turn.  It  is  our  purpose  to 
transform  these  four  vectors  into  inertial  space,  Sj. 

The  gun  system,  S^,,  is  shown  in  Fig.  19,  where  the  x^,  axis  is 
parallel  to  Xg  and  the  z^,  axis  lies  along  the  inner  gimbal  axis.  The 
distance  from  the  origin  of  S^,  to  the  origin  of  Sg  is  represented  by  the 
matrix 
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.3*»rrT-— »w« 


Fiizure  19 

Relation  of  Gun  Space,  S,,,  and  Barrel  Cluster  Space,  S 


(368) 


M 


“GB 


I  n 

|  u  I 
\  0  / 


The  projectile  velocity  in  S^.  is 


UG  =  TGB  +  nGBbB) 


(369) 


where  is  the  3X3  matrix  representation  of  uv~  • 

GB  GB 


(°  ■ 

0  , 

0 

rtGB  ' 

0  , 

0  , 

-<4i 

’ 

0  / 

(370) 


The  relation  of  the  gun  system,  S^.,  with  respect  to  the  turret 
system,  S,p,  is  shown  in  Fig.  20.  The  distance  from  the  origin  of  S^, 


to  that  of  S^.  is 


TG 


(371) 


The  transformation  matrix  relating  vectors  in  S^,  and  S,p  is 
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Inner  Gimbal  Axi 


G 


Figure  20 

Relation  of  Turret  Space,  S^,, 
0„  is  in  the  x^,  z 


and  Gun  Space 
,  Plane 


•**»*,  la 


cos  E1  cos  A1  ,  -sin  E’  co"  A1 ,  -sin  A' 

sin  E'  ,  cos  E'  ,  0 

cos  E1  sin  A'  ,  -sinE'  sin  A'  ,  cos  A' 


(372) 


and  the  position  of  the  gun  muzzle  as  measured  in  S^,  is 


bT  '  bTG  +  TTG  bG 


(373) 


The  direction  of  the  gunbore  as  measured  in  S,^  is 


VT  “  TTGVG 


(374) 


and  the  angular  velocity  of  the  projectile  in  S^,  is 


WT  =  TTG(wTG  +  wG} 


(375) 


where 


/-A* 


A'sin  E 


A'  cos  E’ 


(376) 


in  S^.  coordinates.  The  projectile  velocity  in  S^,  is  given  by 


UT  "  UTG  +  TTG  ^UG  +  ^TGbG^ 


(377) 


where 
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and 


(378) 


(379) 


Since  both  aircraft  space,  S^,  and  turret  space,  S,p,  are  fixed  in 
the  aircraft,  the  position  of  the  gun  muzzle  in  S.  is 


bA  =  bAT  +  TAT  bT 


(380) 


Whereas  it  is  possible  to  define  T^^  in  terms  of  orientation  angles 
relating  S^,  and  S^,  this  will  serve  no  particular  purpose  here  and 
will  be  omitted.  T^^  can  be  defined  in  any  particular  application  when 
the  need  arises.  The  vector  b^-j.  is  the  distance  from  the  airciaft  center 
of  gravity  to  the  origin  of  Srp  and  is  measured  in  S^.  The  direction  of 
the  gunbore  in  is 


VA  =  TATv1 


The  angular  velocity  of  the  projectile  in  S  is 

xx 


W, 


=  T 


(381) 


(382) 
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and  the  projectile  velocity  in  is 


UA  =  aATuT 


(383) 


The  position  of  the  gun  muzzle  in  inertial  space,  S^,  is 


bI  =  bIA  +  TLAbA 


(384) 


where  bTA  is  the  position  of  the  aircraft  center  of  mass  inST,  and  T 

I  IA 

is  the  matrix  which  relates  the  relative  orientations  of  ST  and  S.  .  Writ- 
ing  Tjj^  in  terms  of  orientation  angles  of  the  aircraft  would  serve  no 
purpose  here  and  will  be  omitted.  The  direction  of  the  gunbore  in  Sj  is 


VI  =  tia  va 


(385) 


and  the  angular  velocity  of  the  projectile  in  is 


“I  =  TIA("!A  +  “A> 


(386) 


where  is  the  angular  velocity  of  with  respect  to  as  measured  in 


and 


*1A 


(387) 


The  velocity  of  the  projectile  in  S  ^  is 


UI  =  UIA  +  llA  (uA  +  °IA  bA) 
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'  - - -  - -  - „■  SOW.' 


where 


0  -r  ,  i 

t 

T  0  ,  -p 

9 

■q  .  p  .  o 


(388) 


Parameters  p,  q,  and  r  are  the  usual  symbols  for  the  components  of 
aircraft  ang>  lar  velocity  in  the  x^,  y^,  and  z ^  directions,  respectively, 
and  are  ass.  Tied  to  be  known.  The  matrix  Uj^  is  equivalent  to  the  air¬ 
craft  velocity  vector  V^,  i.  e. , 


(389) 


Usually,  the  aircraft  position  matrix,  b^,  is  of  no  interest  (only  rela¬ 
tive  target  position  is  needed)  and  the  origin  of  Sj  is  taken  to  be  at  the 
instantaneous  position  of  the  aircraft  center  of  mass  at  the  time  of  fire. 
Thus 


bIA  =  0 


(390) 


and  bj  is  the  location  of  the  gun  muzzle  with  respect  to  the  aircraft  cen¬ 
ter  of  mass.  The  aircraft  altitude  is,  of  course,  assumed  to  be  known. 
The  symbol  used  in  the  next  subsection  for  the  initial  bullet  position  at 


the  time  of  fire  is  R  ,  so 

o 


(391) 


Al&v.  the  projectile  velocity  and  angular  velocities  at  time  of  fire  are, 
respectively, 


u  =  uT 
o  I 


(392? 


and 


o 


(393) 


and  the  direction  of  the  gunbore  or  of  the  shell  longitudinal  axis  is 


e 

o 


(394) 


4.  Initial  Conditions 

In  Subsection  3,  coordinate  systems  and  transformations  were 
introduced  for  the  calculation  of  projectile  velocity,  uq,  projectile 
angular  velocity,  u>q,  gun  muzzle  position,  Rq,  and  boreline  direction, 
eQ,  in  inertial  space  at  the  time  of  fire.  It  will  be  assumed  that  these 
vectors  are  available  (defined  in  inertial  space)  so  that  appropriate 
initial  conditions  may  be  derived  for  trajectory  computations  using  the 
equations  of  Sections  IV  and  V. 

4.1  Initial  Conditions  for  the  Matrix  Formulation  -  Required 

initial  conditions  for  use  with  the  matrix  formulation  of  the  equations 

of  motion.  Section  IV,  Subsection  2,  are  values  at  time  t  for  the  com- 
_  _  ° 
ponents  of  u,  u,  Rq,  and  the  direction  cosine  matrix 


By  definition, 


A 


x 


1 


e 

o 
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and  from  Eq.  (159)  it  is  seen  that  initial  values  ofi^,  an<*  ^3  are 
defined.  Also,  initial  position  Rq  is  given.  Other  initial  parameters 
depend  upon  the  choice  of  the  and  x^  directions  of  the  x^x^.x^ 
coordinate  system  and  a  convenient  initial  orientation  is  that  shown  in 
Fig.  21.  It  is  seen  that 


cos  p'  sina1, 
-sin(3!  sina’, 
-  cos  a'  , 


Comparison  of  this  equation  with  Eq.  (159)  shows  that 

£  j  =  cos  (3'  cos  a' 

£  2  =  cos  p'  sin  a' 

£  3  =  sin  p' 


and  it  follows  that 


cos  p' 


cos  a' 


sin  a' 


JT 


c 


2 

3 


£ 


2 


-  I 


2 

3 
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Relation  between  the  x,,x,,x,  Syste 


F igure  2 1 

vV-3 
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Given  a'  and  (3',  the  other  direction  cosines  may  be  calculated  and  so 

A  '-3  defined.  Components  of  u  and  w  in  the  x,,x^,x_  coordinate 

system  are  given  by  the  matrices  Auq  and  Auq,  respectively,  where  uq 

and  go  are  the  inertial  space  matrix  representations  of  H  and  w  ,  re- 
o  o  o 

spectively.  This  completes  the  set  of  initial  conditions  for  the  matrix 
equation  formulation  of  the  six-degree -of -freedom  equations. 

4.  2  Initial  Conditions  for  the  Euler  Angle  Formulation  -  Initial 
conditions  for  the  Euler  angle  formulation  of  Section  IV,  Subsection  3 
are  specified  by  values  for  the  following  parameters:  qq,  9q,  uq, 
o)po,  w30»  6q,  <{>0»  il0»  and  £  .  For  convenience,  the  origin  of  the 
ii  t  tial  space,  Sj,  is  taken  at  the  instantaneous  position  of  the  gun  muz- 
7,1  at  the  time  of  fire,  so  initial  position  coordinates  £  ,  q^,  and 
are  zero.  As  was  explained  in  Section  IV,  the  £,  q,  t,  axes  are  defined 
such  that  q  is  vertical,  positive  up,  and  the  £  axis  lies  along  the  pro¬ 
jection  of  uq  in  the  horizontal  plane.  If  the  Xj,  y^,  axes  are  defined 
as  shown  in  Fig.  22,  the  relation  between  the  two  systems  is  given  by 


In  this  equation,  an  arrow  over  a  symbol  denotes  a  unit  vector  in  the  direc¬ 
tion  associated  with  the  symbol.  If  u  ,  u  ,  and  u  are  the  components  of  u  in 

7  x  y  z  o 

the  Xj,  yj,  and  Zj  directions,  respectively, 


(396) 
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(397) 


where 


sin0 


cos  6  = 


(398) 


2  2 

u  +  u 
x  y 


(399) 


sinB1  = 


(400) 


cosB'  = 


(401) 


By  definition,  aQ  =  0.  If  the  components  of  eQ  in  the  Xj.yj,  Zj  system 

aTe  e  i  6  ,  and  e  , 
x  y  2 


u  •  e  , 

ros  s  =  _2 - 2  =  -  (ue  +ue  +ue) 

cos  °o  u  u  '  x  x  y  y  z  z 

o  o 


(402) 


sin  6  -I  1  -  cos  6 
on/  o 


(403) 
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An  expression  for  the  angle  4>q  can  be  obtained  with  the  aid  of  Fig.  23. 
In  the  figure,  tj*  is  in  the  vertical  plane  and  is  perpendicular  to  uq. 
Hence,  4>q  is  the  angle  between  t)'  and  the  projection,  e,  ,  of  eQ  in  the 
tj ' ,  C,  plane.  From  Eq.  (395),  the  components  of  eQ  in  the  £,ti,  £  coor¬ 
dinate  system  are 


Hence, 


e„  =  e  sin  B'+  e  cos  B' 
£  x  y 


(104) 


e  =  e  (405) 

q  z 


e.  =  e  cos  B'  -  e  sin  B' 
L  x  y 


(406) 


e  ,  =  -e.  sin  0  f  e  cos  8  (407) 

T)'  |  O  T|  o 


and  with 


e 


1 


(408) 


it  follows  chat 


sin  <J> 
ro 


>409) 


COS  d> 

Yo 


From  Fig.  22, 


(410) 


i  68 


Figure  23 


Geometiy  Defining  4> 
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(414) 


,  sin  6  cos  <i>  , 

o  To 

,  cos  6  cos  6  , 

o  To 

,  -sintb  , 
To 


sin  6  sin<b 
o  rol 


cos  6  sin  <b 
c  To 


cos  4> 


This  completes  the  set  of  initial  conditions  for  Euler  angle  formulation. 

4.  3  Initial  Conditions  for  the  Approximate  Equations  -  Initial  con¬ 
ditions  for  the  approximate  equations  of  Section  IV,  Subsection  4  are  de¬ 
veloped  in  the  same  manner  as  those  for  the  preceding  subsection.  New 
parameters  are 

P  =  0 

4 

P  =  u 

o 


Parameters  uq,  6q,  4>o>  9o,  and  B'  are  calculated  as  before. 

4.4  Siacci  Calculations  -  No  integration  is  required  when  the 
Siacci  equations  are  used.  However  uq,  0o>  B',  and  4>q  are  needed  and 
they  are  calculated  as  in  Subsection  4.2. 
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SECTION  VII 


FURTHER  COMMENTS 

Whereas  this  report  covers  most  of  the  aspects  of  exterior  bal¬ 
listics  which  are  of  interest  in  airborne  applications,  some  items  are 
excluded  because  they  are  adequately  covered  elsewhere.  Subjects 
omitted  for  this  reason  include  numerical  integration,  and  the  standard 
atmosphere. 

Methods  of  numerical  integration  are  required  since  it  is  believed 
to  be  impossible  to  obtain  a  closed-form  solution  to  the  projectile  equa¬ 
tions  of  motion.  Techniques  for  numerical  integration  are  well  known, 
and  can  be  found  in  many  of  the  standard  references  such  as  Ref.  68. 
Alternately,  specialized  methods  are  applicable  and  one  such  method 
may  be  found  in  Ref.  69. 

Since  all  trajectory  calculations  involve  the  air  density  and  the 
speed  of  sound  vs  altitude,  a  means  of  calculating  these  quantities  is 
needed.  This  need  is  met  by  Pef.  70  and  a  convenient  model  for  com¬ 
puter  application  can  be  found  in  Ref.  58. 

This  report  is  chiefly  concerned  with  the  application  of  more  or 
less  standard  methods  to  airborne  fire  control.  New  methods  applicable 
to  treating  windage  jump  may  be  found  in  Section  IV,  however,  and  in 
Refs.  71  and  72,  which  contain  a  simplified  set  of  approximate  equations 
adaptable  to  onboard  utilization.  Methods  developed  in  this  report  and 
in  Refs.  71  and  72  essentially  comprise  the  current  state  of  the  art  in 
exterior  ballistics  for  airborne  applications.  Models  developed  to  date 
are  applicable  to  the  20-mm,  M56  round.  Development  of  new  methods 
for  treating  windage  jump  may  be  needed  for  new  rounds,  such  as  the 
new  2  0-mm,  25 -mm,  and  30-mm  rounds  under  c evelopment. 
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APPENDIX 


MATRIX  NOTATION 


The  matrix  notation  explained  here  is  useful  in  defining  the 
coordinate  transformations  and  equations  which  are  necessary  for 
setting  up  initial  conditions  for  trajectory  computations  onboard  an 
aircraft. 

In  Fig.  24,  JetS  denote  a  right-handed,  rectangular  coordinate 
system  with  origin  O  and  coordinate  axes  x,y,z.  A  point  in  space,  as 
measured  in  S,  has  coordinates  x,y,  and  z  and  can  be  represented  as 
a  column  matrix 


or  as  a  vector 


a 


ex+ey+ez 
x  y  z 


where  e  ,  e  ,  and  e  are  unit  vectors  in  the  indicated  directions.  If 
x  y  z 


e 


) 


defines  a  row  matrix,  then  matrix  multiplication  provides  a  connection 
between  the  matrix  and  vector  notation  for  a  point  in  S:. 

a  =  ea  (415) 


173 


in  Fig.  24,  let  S'  denote  another  right-handed, rectangular  coordinate 


system  with  origin  O'  and  coordinate  axes  x',y',  z*.  If  the  matrix  c,  or 
the  vector  cT  =  ec,  denotes  the  position  of  O'  as  measured  inS,  the 
transformation  of  the  coordinates  of  a  point  P,  represented  by  a1  in  S-, 
into  the  coordinates  of  P,  represented  by  a  in  S,  is  given  by  the  relation 

a  =  c  +  T  a«  (416) 

T  is  a  3  X  3  matrix  defining  the  angular  orientation  of  the  S*  system  with 
respect  to  the  S  system.  Ii  e*  is  the  row  matrix  of  unit  vectors  defining 
the  x'.y'.z'  directions  inS1,  that  is,  if 


then 


e'  =  eT 


(417) 


and  also 

e  =  Te'  (418) 

where  the  tilde  (~)  over  a  symbol  denoting  a  matrix  represents  the  trans¬ 
pose  of  that  matrix. 

It  is  a  known  property  of  an  orthogonal  transformation  matrix,  T, 
that  the  inverse,  T  ^ ,  equals  the  transpose,  Y.  That  is 

W  _  ~  1 

or 

T  T  =  I  (419) 
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where  I  is  the  unit  matrix.  It  follows  from  Eq.(4l6)  that 

a*  -  T  (a  -  c) 

The  relations  above  define  transformations  of  position  coordinates 
between  S  andS1.  Equations  for  velocity  transformations  can  be  obtained 
by  differentiation  with  respect  to  time.  Differentiation  of  Eq.  (419),  where 
a  dot  denotes  time  diffe  entiation,  yields  the  result 

?T  +  f  T  =  o  (420) 


With 


12'  =  ?T 


(421) 


it  is  seen  that  £2'  =  ¥  T  and  so  Eq..  (420)  may  be  written  as 


£2'  +  £1'  =  0 


It  follows  that  the  principal  diagonal  of  £2*  is  zero  and  £2?^  =  -£2j^,  i.e., 
£2'  is  skew -symmetric,  and  can  be  written  in  the  form 


It  can  be  shown  that  the  components  of  the  column  matrix  £2* a'  equal  the 
corresponding  components  of  the  vector  £2'  X  a'.  That  is 

e'£2'  a'  =  £2'  X  a' 
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where 


fl'  =  eti*f  +  etn.  +  e,n, 

x*  x'  y*  y’  z'  z * 


Hence,!!’  is  identified  as  the  angular  velocity  of  S’  relative  to  S 
measured  in  S'.  It  follows  from  Eq.  (421)  that 


T  =  Tfi» 


(422) 


Differentiation  of  b  =  Tb'  yields 


b  =  Tb'  +  Tb' 


where  b  ~  eb  is  a  position  vector.  Hence,  from  Eq.  (422) 


b  =  Tb'  +  TQ'b1  =  T(b'+n’b') 


or,  since  e'  =  eT 


eb  =  e*  b*  +  e'  fi'b’ 


This  is  equivalent  to  the  familiar  vector  relation 


+  S2  X  b 


where  subsc.ipts  S  and  S'  denote  the  spaces  in  which  the  velocities  are 
observed.  If  e'b'  =  b'  represents  a  point  P  in  S',  then  e'b'  is  the 
velocity  of  P  as  observed  in  S',  e  b  is  the  velocity  of  P  as  observed  in  S, 
and  e’fl'b1  is  the  component  of  velocity  of  P  in  S  due  to  the  angular 
motion  of  S'. 
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Note  that  if  e  is  constant  (i.e. ,  if  S  is  inertial),  then  from  Eqs.  (417) 


and (422) 


e'  =  eT  =  eTQ'  =  e'Sl* 


If  the  column  vector  representation  of  Q*  in  S'  is 


u  =  1  OJ' 


(423) 


is  the  column  vector  representing  the  components  of  ft'  in  the  S  coordinate 
system.  This  follows  from  the  fact  that  there  are  vectors  SI  and  b  in  S 
such  that  0  X  b  =  «'  X  b*.  Hence, e£2b  =  e'n'b'  =  (eT)fl'('i'b)  = 
e(TQ*  T)b  and  so  SI  =  Tft'¥.  That  Eq,  (423)  follows  from  fit  =  TO' T  can 
be  shown  by  comparison  of  matrix  elements. 

If  U  is  the  matrix  representing  a  transformation  from  a  space  S" 


to  S',  that  is,  if 


a'  =  Ua" 


and  if  T  is  the  matrix  representing  a  transformation  from  S'  to  S,  that  is,  if 


a  =  T  a' 


then 
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a  -  TUa"  =  W  a” 


W  =  TU 


is  the  matrix  representing  a  transformation  from  S"  to  S.  Then 


W  =  TU  +  TU 


(424) 


FromEq.  (422),  we  may  write 


W  =  Wfi, 


T  =  Tfi, 


u  =  unT 


where  components  of  l and  £2^  are  written  with  respect  to  Sn  and 
components  of  T  are  written  in  S'  coordinates.  Substitution  of  these 
expressions  into  Eq.  (424)  yields 


Wfiw  =  T£2XU  +  T  UQy 


W«w  =  TUUQtU  +  TUfly 


Since  W  =  TU  and  W  =  U  T,  it  follows  that 


UfiwU  =  f2T  +  UQyU 
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If  and  cjy  are  the  vector  representations  of  QT,  and 

respectively,  then 

U<Jyy  =  Cdrp  +  U  Wy 

This  may  be  identified  with  the  vector  equations 

Uty  =  (JT  +  Uy 

This  development  is  a  proof  that  angular  velocities  add, which  is  a  result 
that  is  not  intuitively  obvious. 
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fire  control  system  design  are  documented  herein.  The  fundamentals  of  exterior 
ballistics  are  included,  along  with  a  description  of  currently  used  ballistic  and 
aerodynamic  notations,  a  discussion  of  the  limitations  of  the  semi-empirical 
aerodynamic  force  and  moment  system,  and  methods  of  preparing  aerodynamic 
data  for  use  in  trajectory  computation.  Tutorial  material  is  provided  to  give 
the  reader  an  understanding  of  windage  jump  caused  by  the  complicated  angular 
motion  of  a  spinning  projectile.  The  Siacci  method  is  described  and  means  for 
improving  its  accuracy  are  developed.  Six-degree-of-freedom  equations  are 
derived  in  several  different  formulations  for  exploratory  studies  and  digital 
computer  computations,  and  the  development  of  approximate  equations  for  rapid 
evaluation  of  trajectory  tables  is  included.  Methods  for  calculating  trajectory 
initial  conditions  are  provided  for  shell  fired  from  a  turreted,  gatling  gun  in  a 
maneuvering  aircraft,  and  the  problems  of  ballistic  and  kinematic  prediction 
are  discussed  briefly.  The  material  covered  herein  should  provide  personnel 
in  the  Air  Force  and  in  industry  with  sufficient  knowledge  of  exterior  ballistics 
for  advanced  fire  control  system  design. 
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